Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,950,559 (C)5→6 intergenic (‑37/‑273) aspS ← / → yecD aspartyl‑tRNA synthetase/isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,5541.C96.3% 280.6 / 2.0 81intergenic (‑32/‑278)aspS/yecDaspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (38/40);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

ATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950422‑1950683
                                                                                                                                      |                                                                                                                                 
aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccc                                                                                                                             >  2:1051437/1‑138 (MQ=255)
aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccc                                                                                                                             >  2:1244094/1‑138 (MQ=255)
aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccc                                                                                                                             >  1:573982/1‑138 (MQ=255)
  cAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           >  1:545148/1‑136 (MQ=255)
    ggCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctgct                                                                                                                         >  2:974751/1‑134 (MQ=255)
    ggCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctgct                                                                                                                         >  1:482479/1‑134 (MQ=255)
     gCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCctgctg                                                                                                                         >  1:1327090/1‑139 (MQ=255)
          cAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAg                                                                                                                    <  2:573982/139‑1 (MQ=255)
               tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGAc                                                                                                               >  1:731943/1‑139 (MQ=255)
               tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccc                                                                                                                             <  2:870529/124‑2 (MQ=255)
               tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccc                                                                                                                             >  1:870529/1‑123 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCctgctg                                                                                                                         <  2:528954/125‑1 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCctgctg                                                                                                                         >  1:528954/1‑125 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAg                                                                                                           >  1:583471/1‑139 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAg                                                                                                           >  1:91083/1‑139 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑ccc                                                                                                                                >  2:1783883/1‑117 (MQ=255)
                   gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑ccc                                                                                                                                <  1:1783883/117‑1 (MQ=255)
                     cgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         <  2:731943/139‑1 (MQ=255)
                     cgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                         <  1:1244094/139‑1 (MQ=255)
                       tcgGTTGACCCAACAACACAGAGCCACCTGCTGCCCCACGTGGGACAAACGCAGCTGTCCACAATATTCTGTACGCATGAGTTATCCCTTGACTTAGCTGCCGGCGGATCC‑CCCCCTCCTGCGCAGGTGACCAAGTCGc                                                                                                       >  1:1608405/2‑139 (MQ=255)
                       acgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       <  2:393275/139‑1 (MQ=255)
                        cgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCa                                                                                                      >  2:712745/1‑139 (MQ=255)
                               cccAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               <  2:1589160/139‑1 (MQ=255)
                                  aaCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:974751/139‑1 (MQ=255)
                                    ccACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                          >  2:875085/1‑139 (MQ=255)
                                            agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                  <  2:91083/139‑1 (MQ=255)
                                              tCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           >  2:1549383/1‑92 (MQ=255)
                                              tCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           <  1:1549383/95‑4 (MQ=255)
                                                      tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               <  1:918729/116‑1 (MQ=255)
                                                      tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               >  2:918729/1‑116 (MQ=255)
                                                        ccccACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTATATTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaaaa                                                                      >  1:825270/1‑139 (MQ=255)
                                                            aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  1:1051437/139‑1 (MQ=255)
                                                            aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:482479/139‑1 (MQ=255)
                                                             cGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggcg                                                                 >  2:1141503/1‑139 (MQ=255)
                                                                  gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                  >  2:1793070/1‑117 (MQ=255)
                                                                  gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                  <  1:1793070/117‑1 (MQ=255)
                                                                  gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  1:931594/133‑1 (MQ=255)
                                                                  gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  2:931594/1‑133 (MQ=255)
                                                                   aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  2:578087/1‑132 (MQ=255)
                                                                   aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  1:578087/132‑1 (MQ=255)
                                                                          gAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTgg                                                    >  1:209274/1‑139 (MQ=255)
                                                                                     aaTATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtgaa                                                                             <  1:330943/103‑1 (MQ=255)
                                                                                     aaTATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtgaa                                                                             >  2:330943/1‑103 (MQ=255)
                                                                                     aaTATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCg                                         <  1:1141503/139‑1 (MQ=255)
                                                                                     aaTATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCg                                         <  2:583471/139‑1 (MQ=255)
                                                                                          tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                    <  1:551606/75‑1 (MQ=255)
                                                                                          tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                    >  2:551606/1‑75 (MQ=255)
                                                                                            tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  <  2:1316832/139‑1 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  2:81511/104‑1 (MQ=255)
                                                                                               aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  1:81511/1‑104 (MQ=255)
                                                                                                 gCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               <  2:1129688/73‑1 (MQ=255)
                                                                                                 gCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               >  1:1129688/1‑73 (MQ=255)
                                                                                                 gCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGc                             >  2:125012/1‑139 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACTACTGAATGAAAAAAg                                                                    <  1:1200096/98‑1 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAg                                                                    >  2:1200096/1‑98 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                           >  2:795768/1‑139 (MQ=255)
                                                                                                   aTGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                           >  2:640456/1‑139 (MQ=255)
                                                                                                     gagaTATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           <  1:159132/40‑4 (MQ=255)
                                                                                                     gagaTATCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           >  2:159132/1‑37 (MQ=255)
                                                                                                       gaTATCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                       <  1:876784/139‑1 (MQ=255)
                                                                                                        atatCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTca                                                                                     <  1:559245/76‑1 (MQ=255)
                                                                                                        atatCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTca                                                                                     >  2:559245/1‑76 (MQ=255)
                                                                                                        atatCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           >  1:534527/1‑34 (MQ=255)
                                                                                                        atatCCCTTAACTTAGCTGCCGGCGGAT‑G‑cccccctg                                                                                                                           <  2:534527/37‑4 (MQ=255)
                                                                                                           tCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             <  2:151049/97‑1 (MQ=255)
                                                                                                           tCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             >  1:151049/1‑97 (MQ=255)
                                                                                                           tCCCTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCcgcg                   >  2:585111/1‑139 (MQ=255)
                                                                                                             ccTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       <  1:823514/53‑1 (MQ=38)
                                                                                                             ccTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                       >  2:823514/1‑53 (MQ=38)
                                                                                                             ccTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  >  1:247544/1‑58 (MQ=39)
                                                                                                             ccTTAACTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                  <  2:247544/58‑1 (MQ=39)
                                                                                                                 aaCTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAAc                                                                                 <  2:1790122/71‑1 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAAc                                                                                 >  1:1790122/1‑71 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  1:1126034/86‑1 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  2:1126034/1‑86 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCcgcg                   >  1:1349893/1‑133 (MQ=255)
                                                                                                                 aaCTTAGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCcgcg                   <  2:1349893/133‑1 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               <  1:104965/52‑1 (MQ=38)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                               >  2:104965/1‑52 (MQ=38)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             >  1:1620775/1‑86 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             <  2:1620775/86‑1 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGc                                           >  1:890781/1‑104 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGc                                           <  2:890781/104‑1 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  >  2:47891/1‑113 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  <  1:47891/113‑1 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTgg               >  1:1766934/1‑132 (MQ=255)
                                                                                                                      aGCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTgg               <  2:1766934/132‑1 (MQ=255)
                                                                                                                       gCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc       >  1:1050495/1‑139 (MQ=255)
                                                                                                                       gCTGCCGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc       >  2:775551/1‑139 (MQ=255)
                                                                                                                           ccGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  2:825270/139‑1 (MQ=255)
                                                                                                                           ccGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  2:1050495/139‑1 (MQ=255)
                                                                                                                           ccGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  1:775551/139‑1 (MQ=255)
                                                                                                                           ccGGCGGAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  1:712745/139‑1 (MQ=255)
                                                                                                                           ccGGCGGAT‑G‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  >  2:1754456/1‑75 (MQ=255)
                                                                                                                           ccGGCGGAT‑G‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                  <  1:1754456/75‑1 (MQ=255)
                                                                                                                            cggcggAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc  <  2:1327090/139‑1 (MQ=255)
                                                                                                                            cggcggAT‑GCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc  >  1:1649848/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                 
ATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGAT‑G‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950422‑1950683

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