Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,992,083 G→C 8.1% noncoding (35/76 nt) glyW ← tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,992,0830GC 8.1% 249.9 / 13.7 86noncoding (35/76 nt)glyWtRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (29/50);  new base C (4/3);  total (33/53)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

GAACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTTCACC  >  NC_000913/1991934‑1992232
                                                                                                                                                     |                                                                                                                                                     
gAACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGATACGAGACTCGAACTCGCGACCCCGACCTTgg                                                                                                                                                       >  1:479526/1‑150 (MQ=255)
     ggACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAgg                                                                                                                                                  >  1:790453/1‑150 (MQ=255)
        cTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCg                                                                                                                                               <  2:381015/150‑1 (MQ=255)
         tAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGt                                                                                                                                              >  2:543035/1‑150 (MQ=255)
          aGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                                                                                                                             >  2:984570/1‑150 (MQ=255)
          aGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                                                                                                                             >  2:297249/1‑150 (MQ=255)
           gACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCgggg                                                                                                                                            <  2:713385/150‑4 (MQ=255)
            aCGGATTTGCAAGCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCATGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggggt                                                                                                                                           <  2:628637/150‑5 (MQ=255)
            aCGGATATGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggggt                                                                                                                                           <  2:344494/150‑5 (MQ=255)
                aTTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAc                                                                                                                                       <  2:1202679/150‑1 (MQ=255)
                 tttGCAATCCGCTACATAACCGCTTAGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAcc                                                                                                                                      <  1:231833/150‑1 (MQ=255)
                   tGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGCCCCCGACCTTGGCAAGGTCGTGCGCTACCaa                                                                                                                                    >  2:1110612/1‑150 (MQ=255)
                    gCAATCCGCTACATAACCGCTTTCTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAAc                                                                                                                                   <  2:232287/150‑1 (MQ=255)
                      aaTCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTg                                                                                                                                 >  2:407723/1‑150 (MQ=255)
                           gCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTa                                                                                                                            >  1:8444/1‑150 (MQ=255)
                               cATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACTTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTcc                                                                                                                        >  2:1247566/1‑150 (MQ=255)
                               cATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTcc                                                                                                                        <  1:40213/150‑1 (MQ=255)
                                   aCCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCa                                                                                                                    >  1:465280/1‑150 (MQ=255)
                                     cGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCACGCAtt                                                                                                                  >  2:369013/1‑150 (MQ=255)
                                      gCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCGTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtc                                                                                                                 <  2:712818/150‑1 (MQ=255)
                                       cTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtca                                                                                                                <  1:850273/150‑1 (MQ=255)
                                       cTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtca                                                                                                                <  1:297249/150‑1 (MQ=255)
                                          tGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAGACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtcatca                                                                                                             <  1:637873/150‑1 (MQ=255)
                                          tGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtcatca                                                                                                             <  2:995990/150‑1 (MQ=255)
                                          tGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCACGCATtcatca                                                                                                             >  1:1247682/1‑150 (MQ=255)
                                            ttAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTACACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAg                                                                                                           <  2:343784/150‑1 (MQ=255)
                                              aaCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCa                                                                                                         <  1:291863/150‑1 (MQ=255)
                                                   gcCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAg                                                                                                    <  2:47799/150‑1 (MQ=255)
                                                    ccAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGt                                                                                                   <  1:549022/150‑1 (MQ=255)
                                                           cttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACTAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTTTTCCCGCATTCATCAAGCAATCAGTTAATGAc                                                                                            >  2:805370/1‑150 (MQ=255)
                                                           cttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCAc                                                                                            <  1:661633/150‑1 (MQ=255)
                                                            ttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTTAACTGGAGCCGGAAACGAGACTCGAACTCGCGACCACGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCACTCAGTTAATCACt                                                                                           <  1:633857/150‑1 (MQ=255)
                                                              cAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCAGCATTCATCAAGCAATCAGTTAATCACTTg                                                                                         <  2:338036/150‑1 (MQ=255)
                                                                ggCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAt                                                                                       >  2:666634/1‑150 (MQ=255)
                                                                        aGCCAGACATCCGCTTGACGACGATGTCTTTTAAACTGTAGCGGCAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGTTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTACTCACTTGATGTtattat                                                                               <  2:69935/150‑1 (MQ=25)
                                                                         gCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATc                                                                              >  1:318277/1‑150 (MQ=255)
                                                                         gCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGATGAGCTCTACCAACTGAGCTAATCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATc                                                                              <  1:866785/150‑1 (MQ=255)
                                                                         aaaaGACATCGGCGTCAAGCGGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCGATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTTTc                                                                              >  1:1297514/4‑150 (MQ=21)
                                                                          ccAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCg                                                                             <  2:790453/150‑1 (MQ=255)
                                                                              aCATCCGCTTGCCGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTg                                                                         <  1:212711/150‑1 (MQ=255)
                                                                               cATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTgg                                                                        <  2:596706/150‑1 (MQ=255)
                                                                                   cGCTTGACGCCGATGTCTTTTACACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTATACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCaaa                                                                    >  2:919537/1‑149 (MQ=255)
                                                                                    gCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATc                                                                   <  1:766347/150‑1 (MQ=255)
                                                                                           gCCGATGTCTTTTAAACTGGAGCGTGAAACGAGACTCGAACTCGCGACCCCGTCCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCTAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgccg                                                            <  2:1232045/150‑1 (MQ=255)
                                                                                            ccGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgccgc                                                           <  2:847574/150‑1 (MQ=255)
                                                                                               aTGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGt                                                        <  2:591391/150‑1 (MQ=255)
                                                                                                    ttttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTTTTCCCGCATTCTTCAAGCAAGCGGTTTATCTCTTGATTGTATTATCGTCTGGCCATCAGTGCCGCCGTTCGAt                                                   >  2:935569/1‑150 (MQ=21)
                                                                                                     tttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATg                                                  <  2:609759/150‑1 (MQ=255)
                                                                                                      ttATACTGGAGCGGGATATGAGACTCGTTCTCTCGAGCACGACCTTGGCAAGGTCGTGCTCTTCCAAATGAGCTATTCGCGCATTCATCTAGCAATCAGTTAATCTCTTGATTTTATTATGGTCTGGCAATCAGTGCCGCCGTTCGATGc                                                 <  1:867418/150‑1 (MQ=25)
                                                                                                         aaCTGGAGCGGGAAACGAGACGGGAACTCGCGACCCCGACCTTCGCAAGGTCGTGCTCTACCAACTGTGCTATTACCGCATTCATCATGCAATGAGTTAATCACTTGTTTTTATTATCGTCTGGGAATCAGTGCCGCCGTTCGATGCGtt                                              <  1:1147747/150‑1 (MQ=25)
                                                                                                          aCTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTg                                             <  2:586958/150‑1 (MQ=255)
                                                                                                               aGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTc                                        >  1:728571/1‑150 (MQ=255)
                                                                                                                gCGGGAAACGAGACTCGCACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCCTCAAGCAATCAGTTAATCACTTGATTTTATTCTCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCt                                       <  1:201577/150‑1 (MQ=255)
                                                                                                                gCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCt                                       <  2:1285118/150‑1 (MQ=255)
                                                                                                                  gggAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCTCCGTTCGCTGCGTTGCATTCTAc                                     >  1:1286341/1‑150 (MQ=255)
                                                                                                                  gggAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTAc                                     <  1:369013/150‑1 (MQ=255)
                                                                                                                   ggAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCTCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACt                                    >  2:192831/1‑150 (MQ=255)
                                                                                                                   ggAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACt                                    >  2:575499/1‑150 (MQ=255)
                                                                                                                      aaCGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCAATGCGTTGCATTCTACTTAc                                 <  2:1056413/150‑1 (MQ=255)
                                                                                                                        cGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCt                               >  1:334593/1‑150 (MQ=255)
                                                                                                                         gagaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTg                              <  1:806479/150‑1 (MQ=255)
                                                                                                                         gagaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTg                              >  1:497623/1‑150 (MQ=255)
                                                                                                                            aCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGgcg                           <  1:492011/150‑1 (MQ=255)
                                                                                                                             cTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGgcgc                          >  2:554835/1‑150 (MQ=255)
                                                                                                                                 aaCTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATg                      <  2:795097/150‑1 (MQ=255)
                                                                                                                                     cgcgACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGCGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTc                  >  2:726198/1‑150 (MQ=255)
                                                                                                                                     cgcgACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTc                  <  1:1291593/150‑1 (MQ=255)
                                                                                                                                     cgcgACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTc                  >  1:20474/1‑150 (MQ=255)
                                                                                                                                     cgcgACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTc                  <  1:728742/150‑1 (MQ=255)
                                                                                                                                       cgACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCaa                >  2:886808/1‑150 (MQ=255)
                                                                                                                                        gACGCCGACCTTGGCAAGGTCGTGCTGTACCAACTGAGCTATTCCCGAAGTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAAc               <  2:806631/150‑1 (MQ=255)
                                                                                                                                        gACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAAc               <  2:787180/150‑1 (MQ=255)
                                                                                                                                         agcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACg              >  2:538030/5‑150 (MQ=255)
                                                                                                                                         aCCCCGACCTTGGCAAGGTCGTGCTCTACCATCTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACg              <  1:717155/150‑1 (MQ=255)
                                                                                                                                         aCCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACg              <  2:1133153/150‑1 (MQ=255)
                                                                                                                                          ccccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTTCATTCTACTTACCTGGCGCGATGAGTCAACGa             <  1:782215/150‑1 (MQ=255)
                                                                                                                                          ccccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGa             >  1:899795/1‑150 (MQ=255)
                                                                                                                                          ccccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGa             >  2:892059/1‑150 (MQ=255)
                                                                                                                                           cccGACCTTGGCAAGGTCGTGCTCTACCACCTGAGCTATACCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGat            <  1:815194/150‑1 (MQ=255)
                                                                                                                                           cacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGat            >  2:588624/3‑150 (MQ=255)
                                                                                                                                            acgcccTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGata           >  2:86022/5‑150 (MQ=255)
                                                                                                                                            acGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGata           >  1:1059005/2‑150 (MQ=255)
                                                                                                                                                  ttGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAGTTTATTATCGTCTGGCACTCAGTGCCGCCGTGCGATGCGTTGCATTCTACTTACCGGGCGCGATGAGTCAACGATATTTTTc     <  1:512877/150‑1 (MQ=255)
                                                                                                                                                   tGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTTca    <  2:479526/150‑1 (MQ=255)
                                                                                                                                                   tGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTTca    <  1:666634/150‑1 (MQ=255)
                                                                                                                                                   tGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTTca    <  1:1247566/150‑1 (MQ=255)
                                                                                                                                                   tGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTAGCGCGATGAGTCAACGATATTTTTca    <  2:1179423/150‑1 (MQ=255)
                                                                                                                                                     gCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTTcacc  <  2:465280/150‑1 (MQ=255)
                                                                                                                                                     gCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCACTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGATAGATGCGTTGCATTCTACTTACCTGGCGCGATGCGTCCACGATAATGTTcacc  >  1:101891/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTTCACC  >  NC_000913/1991934‑1992232

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