Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,070,777 G→C 6.1% A60G (GCA→GGA)  rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,7770GC 6.1% 223.9 / 4.2 66A60G (GCA→GGA) rutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (24/38);  new base C (2/2);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ACCAACTGCGGGCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGT  >  NC_000913/1070628‑1070922
                                                                                                                                                     |                                                                                                                                                 
aCCAACTGCGGGCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTg                                                                                                                                                   >  2:577381/1‑150 (MQ=255)
 ccAACTGCGGGCATCACCGTGGGGAAACAGTAGTGGTATCCGTTCATAAATTTCTCTGTCCCTGGGGACGTTTGCGGGGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGc                                                                                                                                                  <  2:691217/150‑1 (MQ=255)
  cAACTGCGGGCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGcc                                                                                                                                                 <  2:1233909/150‑1 (MQ=255)
  cAACTGCGGGCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGcc                                                                                                                                                 <  2:73811/150‑1 (MQ=255)
      tGCGGGCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGc                                                                                                                                             >  1:968143/1‑150 (MQ=255)
           gCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTa                                                                                                                                        <  1:107991/150‑1 (MQ=255)
             aTCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTaaa                                                                                                                                      >  2:1000011/1‑150 (MQ=255)
             aTCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTaaa                                                                                                                                      <  2:429382/150‑1 (MQ=255)
              tCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAac                                                                                                                                     <  1:577381/150‑1 (MQ=255)
                aCCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAacac                                                                                                                                   <  1:298415/150‑1 (MQ=255)
                   gTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGTTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGAt                                                                                                                                <  2:264373/150‑1 (MQ=255)
                    tGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATc                                                                                                                               <  2:1302836/150‑1 (MQ=255)
                      gggAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCg                                                                                                                             <  1:1345272/150‑1 (MQ=255)
                       ggAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGt                                                                                                                            >  1:1105419/1‑150 (MQ=255)
                       ggAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGt                                                                                                                            >  2:19719/1‑150 (MQ=255)
                        gAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTCCCGTGCGGGTAAACACGATCCGTg                                                                                                                           >  2:235273/1‑150 (MQ=255)
                        gAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTCCCGTGCGGGTAAACACGATCCGTg                                                                                                                           >  2:230159/1‑150 (MQ=255)
                          aaCAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGcc                                                                                                                         <  2:867825/150‑1 (MQ=255)
                          aaCAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGcc                                                                                                                         <  2:75339/150‑1 (MQ=255)
                           aCAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCg                                                                                                                        >  1:523693/1‑150 (MQ=255)
                            cAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCgg                                                                                                                       >  2:7079/1‑150 (MQ=255)
                              gtagtGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGgag                                                                                                                     <  1:109167/150‑1 (MQ=255)
                               tagtGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGgaga                                                                                                                    <  2:892169/150‑1 (MQ=255)
                                agtGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGaa                                                                                                                   >  2:1238524/1‑150 (MQ=255)
                                 gtGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAAc                                                                                                                  <  1:1253495/150‑1 (MQ=255)
                                    gTAACCGTTCATAAAATTCACTCTCCCAGCCGACGAAAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGt                                                                                                               >  2:1113978/1‑150 (MQ=255)
                                        ccGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGc                                                                                                           <  2:523693/150‑1 (MQ=255)
                                         cGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCt                                                                                                          <  1:1386063/150‑1 (MQ=255)
                                          gTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTg                                                                                                         <  2:968143/150‑1 (MQ=255)
                                           ttCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGa                                                                                                        <  1:1178294/150‑1 (MQ=255)
                                              aTAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGt                                                                                                     >  1:553844/1‑150 (MQ=255)
                                                aaaaTTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCg                                                                                                   <  1:588941/150‑1 (MQ=255)
                                                 aaaTTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCgg                                                                                                  <  2:1344875/150‑1 (MQ=255)
                                                    ttCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGccc                                                                                               <  2:139524/150‑1 (MQ=255)
                                                           tCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTtc                                                                                        <  2:808291/150‑1 (MQ=255)
                                                             ccAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTtcat                                                                                      <  1:435289/150‑1 (MQ=255)
                                                                ggCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAg                                                                                   <  2:1174051/150‑1 (MQ=255)
                                                                  cgacgaTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAggg                                                                                 <  2:1383126/150‑1 (MQ=255)
                                                                         aGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAAt                                                                          <  2:857810/150‑1 (MQ=255)
                                                                              tGACGGGCGCGGTCAGGGTTTTTTGCAGATTGTCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTccc                                                                     <  2:1037969/150‑1 (MQ=255)
                                                                                   ggCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAAc                                                                <  1:1000011/150‑1 (MQ=255)
                                                                                    gcgcGGTCAGGGGTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGAGCAGACGTTCTTTTCCTTCTGCCGTGCTTGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTAGTACCCCATTTCATCTGTGCATAAATTTCCCGTaaag                                                               <  2:728983/150‑3 (MQ=255)
                                                                                      gcgGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTCCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGtc                                                             <  1:230159/150‑1 (MQ=255)
                                                                                      gcgGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTCCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGtc                                                             <  1:235273/150‑1 (MQ=255)
                                                                                           cAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATc                                                        >  2:1074756/1‑150 (MQ=255)
                                                                                             gggTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGc                                                      >  1:1351623/1‑150 (MQ=255)
                                                                                                   tttGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAac                                                                                                                                     <  2:249171/65‑1 (MQ=255)
                                                                                                   tttGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAac                                                                                                                                     >  1:249171/1‑65 (MQ=255)
                                                                                                   tttGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACgg                                                >  1:108856/1‑150 (MQ=255)
                                                                                                   tttGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACgg                                                >  1:657417/1‑150 (MQ=255)
                                                                                                      gCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGgtg                                             <  1:1113978/150‑1 (MQ=255)
                                                                                                       cAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGgtgt                                            >  2:1238341/1‑150 (MQ=255)
                                                                                                           ttGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTcc                                        <  1:1238524/150‑1 (MQ=255)
                                                                                                              ccGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGc                                     <  1:7079/150‑1 (MQ=255)
                                                                                                               cGCTGGAAAGTGCCGGGCGCAGATGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGcc                                    >  1:581478/1‑150 (MQ=255)
                                                                                                                      aaGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAGcgatg                             <  2:1105419/150‑5 (MQ=255)
                                                                                                                         tGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCTGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTtgtg                          >  1:383919/1‑150 (MQ=255)
                                                                                                                         tGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTtgtg                          >  2:835919/1‑150 (MQ=255)
                                                                                                                         ggCCGGGCGCAGACGTGCTTTTCCTTCTGCCGTGCGGGTAAACTCGATGCGTGCCGGAGATCAGTTAGCTGATGGCGGCCCCCATTTGATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGGtgtg                          <  2:1232204/149‑1 (MQ=255)
                                                                                                                                  cAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGc                 >  1:1195565/1‑150 (MQ=255)
                                                                                                                                  cAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGc                 >  2:947462/1‑150 (MQ=255)
                                                                                                                                   aGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCa                <  2:168774/150‑1 (MQ=255)
                                                                                                                                            tttCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAAc       >  1:767005/1‑150 (MQ=255)
                                                                                                                                             ttCCTTCTGCCGTGCGGGTAAACACGATCCATGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACa      >  1:1305089/1‑150 (MQ=255)
                                                                                                                                               ccTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAgg    >  2:1424913/1‑150 (MQ=255)
                                                                                                                                               ccTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAgg    >  2:382578/1‑150 (MQ=255)
                                                                                                                                                 ttCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCTGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGt  <  2:383919/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCAACTGCGGGCATCACCGTGGGGAAACAGTAGTGGTAACCGTTCATAAAATTCACTGTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGT  >  NC_000913/1070628‑1070922

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: