Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,432,994 (C)6→7 coding (567/1125 nt) smf ← DNA recombination‑mediator A family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,9881.C98.5% 234.2 / ‑3.6 67coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (32/34);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑C‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTAC  >  NC_000913/3432850‑3433110
                                                                                                                                           |                                                                                                                           
tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                              >  1:224603/1‑139 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑c‑ccccc                                                                                                                      <  1:257516/139‑1 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑c‑ccccc                                                                                                                      >  2:52677/1‑139 (MQ=255)
               aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAg                                                                                                              >  1:1320385/1‑139 (MQ=255)
               aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑c‑cccccctgtt                                                                                                                 <  2:1138429/135‑6 (MQ=255)
               aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑c‑cccccctgtt                                                                                                                 >  1:1138429/1‑130 (MQ=255)
                           aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGc                                                                                                  >  1:542532/1‑139 (MQ=255)
                           aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGc                                                                                                  <  2:542532/139‑1 (MQ=255)
                                    caccaccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGc                                                                                         <  1:589641/139‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                <  2:224603/139‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                >  1:726688/1‑139 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                <  2:1294120/139‑1 (MQ=255)
                                               agtaaaggtgtAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGa                                                                                                          >  2:731988/11‑111 (MQ=255)
                                               agtaaaggtgtAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGa                                                                                                          <  1:731988/101‑1 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                                 >  2:900833/1‑139 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGATGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             >  2:835014/1‑123 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                                  <  2:1179987/118‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                                  >  1:1179987/1‑118 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             >  2:624312/1‑123 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             <  1:624312/123‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             <  1:835014/123‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg‑c                                                                                                                            >  2:1182781/1‑75 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg‑c                                                                                                                            <  1:1182781/75‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             >  1:25698/1‑122 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             <  2:25698/122‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                            <  1:52677/139‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                            >  2:1345155/1‑139 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:1045800/139‑1 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:726688/139‑1 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  1:900833/139‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCa                                                                                                            <  2:610810/84‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCa                                                                                                            >  1:610810/1‑84 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       <  1:1415244/121‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       >  2:1415244/1‑121 (MQ=255)
                                                                             ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                                                >  2:807500/1‑139 (MQ=255)
                                                                                  cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg‑c                                                                                                                            <  1:1257673/58‑1 (MQ=255)
                                                                                  cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg‑c                                                                                                                            >  2:1257673/1‑58 (MQ=255)
                                                                                      ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGc                                       >  1:557475/1‑139 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       >  2:1021282/1‑103 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       <  1:1021282/103‑1 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     >  2:975466/1‑121 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     <  1:975466/121‑1 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    >  1:905926/1‑138 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    <  2:905926/138‑1 (MQ=255)
                                                                                              tAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:1157687/112‑1 (MQ=255)
                                                                                              tAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          >  1:1157687/1‑112 (MQ=255)
                                                                                                aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       <  2:625221/97‑1 (MQ=255)
                                                                                                aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       >  1:625221/1‑97 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGAc                                                                                                         <  1:772881/58‑1 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGAc                                                                                                         >  2:772881/1‑58 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                <  1:548762/83‑1 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                >  2:548762/1‑83 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCt                        >  2:634458/1‑139 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCt                        >  1:545508/1‑139 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                            <  2:666340/134‑1 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                            >  1:666340/1‑134 (MQ=255)
                                                                                                          ggAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAg                   <  1:634458/139‑1 (MQ=255)
                                                                                                                tCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          >  2:690409/1‑94 (MQ=255)
                                                                                                                tCGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  1:690409/94‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                      <  2:502022/81‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                      >  1:502022/1‑81 (MQ=255)
                                                                                                                  gAGGGGAAATTCCGGGACGAGTGCG‑CACCCCCTGTTCAAACCGACTGGCAGGCAGTAGGGCATCAGGGCGGGGATCCATGGTACTAGGTCCATTCCCCAATACAGCAATGCGGACGCCAGTTAGGTGTTAGGCTGCttt           >  2:880685/1‑139 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                            <  2:1103740/57‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                            >  1:1103740/1‑57 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     <  2:893804/96‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     >  1:893804/1‑96 (MQ=255)
                                                                                                                     gggAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg        <  1:807500/139‑1 (MQ=255)
                                                                                                                      ggAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                  <  1:152122/112‑1 (MQ=255)
                                                                                                                      ggAAATTCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                  >  2:152122/1‑112 (MQ=255)
                                                                                                                           ttCCGAGACGAGAGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTac  <  2:557475/139‑1 (MQ=255)
                                                                                                                              cgagacgagaGCGCC‑CCCCCTGTTCAAGCAGACTGGc                                                                                                     <  2:652981/37‑1 (MQ=255)
                                                                                                                              cgagacgagaGCGCC‑CCCCCTGTTCAAGCAGACTGGc                                                                                                     >  1:652981/1‑37 (MQ=255)
                                                                                                                                agacgagaGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGACGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    <  1:171646/100‑1 (MQ=255)
                                                                                                                                agacgagaGCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGACGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    >  2:171646/1‑100 (MQ=255)
                                                                                                                                      gagCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          >  1:1351370/1‑72 (MQ=255)
                                                                                                                                      gagCGCC‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:1351370/72‑1 (MQ=255)
                                                                                                                                           |                                                                                                                           
TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑C‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTAC  >  NC_000913/3432850‑3433110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑