Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 344,704 A→T 6.3% K177* (AAA→TAA) ‡ yahL → uncharacterized protein YahL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913344,7040AT 6.3% 227.7 / 5.5 80K177* (AAA→TAA) ‡yahLuncharacterized protein YahL
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base T (2/3);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

AGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATG  >  NC_000913/344555‑344853
                                                                                                                                                     |                                                                                                                                                     
aGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTa                                                                                                                                                       >  2:913808/1‑150 (MQ=255)
   atgatgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAg                                                                                                                                                    <  1:169613/150‑1 (MQ=255)
      atgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                                                                                 >  1:123834/1‑150 (MQ=255)
        gCGTTTCTTGTTATTAAAGCAAACGTGCAGTCTTGACGCTTTACCGGGCTTAAGTAATTAGTGCAGTTTATCGATTCAGTGCGAACATATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                               <  2:988110/150‑1 (MQ=255)
         ggTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCATTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                                                                              <  2:1150512/149‑1 (MQ=255)
         cGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                                                                              >  2:7442/1‑150 (MQ=255)
         cGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                                                                              <  1:913808/150‑1 (MQ=255)
            ttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTTTAGGGaaaattc                                                                                                                                           <  2:1425352/150‑5 (MQ=255)
             tCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgtg                                                                                                                                          <  2:281727/150‑1 (MQ=255)
              cTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTTTAGGGaaaattctc                                                                                                                                         <  2:860504/150‑7 (MQ=255)
              cTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgtgt                                                                                                                                         >  1:755126/1‑150 (MQ=255)
               ttGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                                        <  2:1302686/150‑1 (MQ=255)
               ttGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                                        <  1:209300/150‑1 (MQ=255)
               ttGATATTAAAGCAAACGTGCACGCGGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGATTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTTTAGGGaaaattctca                                                                                                                                        <  1:805846/150‑8 (MQ=255)
                   tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCa                                                                                                                                    <  1:18814/150‑1 (MQ=255)
                   tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCa                                                                                                                                    <  1:751609/150‑1 (MQ=255)
                   tatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCa                                                                                                                                    <  1:1319060/150‑1 (MQ=255)
                    atTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAt                                                                                                                                   >  2:1456963/1‑150 (MQ=255)
                     ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtg                                                                                                                                  >  1:761429/1‑150 (MQ=255)
                     ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtg                                                                                                                                  >  1:896181/1‑150 (MQ=255)
                     ttAAAGAAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtg                                                                                                                                  <  1:410094/150‑1 (MQ=255)
                         aGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGa                                                                                                                              <  1:607344/150‑1 (MQ=255)
                               cGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATg                                                                                                                        >  2:997144/1‑150 (MQ=255)
                                gTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATgg                                                                                                                       <  2:76738/150‑1 (MQ=255)
                                    aGGCTGGCCGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGAAAATGTATTCTTCTTGAATACAATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGGATGTGGAATAATGGCtta                                                                                                                   >  2:92838/1‑149 (MQ=255)
                                    aGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTc                                                                                                                   >  2:1304770/1‑150 (MQ=255)
                                         ggACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaa                                                                                                              >  2:421858/1‑150 (MQ=255)
                                         ggACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaa                                                                                                              >  1:1535828/1‑150 (MQ=255)
                                         ggACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGAGaa                                                                                                              >  2:719835/1‑150 (MQ=255)
                                              cTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATACATGTGGAAAAATGGCTTCGCGAAATgtc                                                                                                         >  2:536234/1‑150 (MQ=255)
                                                 tACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgtc                                                                                                      <  1:544349/150‑1 (MQ=255)
                                                      ggCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTAtta                                                                                                 >  2:430181/1‑149 (MQ=255)
                                                         ttAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCgg                                                                                              <  2:830787/150‑1 (MQ=255)
                                                          tAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                             <  1:471656/150‑1 (MQ=255)
                                                           aaGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgag                                                                                            >  2:673601/1‑150 (MQ=255)
                                                            aGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgaga                                                                                           >  1:311359/1‑150 (MQ=255)
                                                             gAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTTAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagac                                                                                          <  1:1308789/150‑1 (MQ=255)
                                                             gAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagac                                                                                          <  2:1424465/150‑1 (MQ=255)
                                                                        aGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                                                               <  1:117382/150‑1 (MQ=255)
                                                                           ttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTa                                                                            <  1:1332568/150‑1 (MQ=255)
                                                                            tATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGAGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTaa                                                                           >  2:1322778/1‑150 (MQ=255)
                                                                             aTCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAg                                                                          <  2:123834/150‑1 (MQ=255)
                                                                              tCGATTCAGTGCGAACCGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGCGACAGACGCGCCCATTAAGc                                                                         >  2:1128168/1‑150 (MQ=255)
                                                                              tCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGc                                                                         <  2:118311/150‑1 (MQ=255)
                                                                                    cAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgcgc                                                                   >  2:699440/1‑150 (MQ=255)
                                                                                     aGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGCCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCt                                                                  <  1:1394485/150‑1 (MQ=255)
                                                                                         cGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgcg                                                              <  1:903323/150‑1 (MQ=255)
                                                                                          gAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                                                                 <  1:43050/82‑1 (MQ=255)
                                                                                          gAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                                                                 >  2:43050/1‑82 (MQ=255)
                                                                                             cAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATACATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATg                                                          <  1:536234/150‑1 (MQ=255)
                                                                                               gATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgag                                                                                            >  2:1255521/1‑114 (MQ=255)
                                                                                               gATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgag                                                                                            <  1:1255521/114‑1 (MQ=255)
                                                                                                   gCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattat                                                    >  2:481475/1‑150 (MQ=255)
                                                                                                     aaaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTAATAAAGTGGAAAAATGGCTTCTCGAAATATCGTATAATCAAGAAACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatta                                                  >  2:1158705/1‑150 (MQ=255)
                                                                                                       aTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt                                                >  2:346788/1‑150 (MQ=255)
                                                                                                         gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTac                                              >  2:768935/1‑150 (MQ=255)
                                                                                                         gTATTCTTCTTGAATAACATATGCCTTTTCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTAGATAAATGTCTTCTATTAAGGATAAAGACGAGCCAATTAAGCAGGCGCTTGCGCATGGATTATTATTac                                              >  2:216167/1‑150 (MQ=255)
                                                                                                           aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTacac                                            <  1:768935/150‑1 (MQ=255)
                                                                                                            ttcttcTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACt                                           >  2:803918/1‑150 (MQ=255)
                                                                                                            ttcttcTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACt                                           <  2:148250/150‑1 (MQ=255)
                                                                                                              cttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                                                                       <  1:1068189/120‑1 (MQ=255)
                                                                                                              cttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                                                                       >  2:1068189/1‑120 (MQ=255)
                                                                                                                tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTaaa                                       <  1:1456963/150‑1 (MQ=255)
                                                                                                                 cttGAATACCATATGCCTTTTCTGAGAATTTTCAATAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATCTCGTATATTCCCGAGCCAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTaaaa                                      >  2:951362/1‑150 (MQ=255)
                                                                                                                      aTACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAatat                                 <  1:673601/150‑1 (MQ=255)
                                                                                                                       tACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAatata                                >  2:1150618/1‑149 (MQ=255)
                                                                                                                         ccATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTa                              >  2:115466/1‑150 (MQ=255)
                                                                                                                         ccATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTa                              >  1:9839/1‑150 (MQ=255)
                                                                                                                          cATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTAc                             <  1:7442/150‑1 (MQ=255)
                                                                                                                            tatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTAccc                           >  2:156099/1‑150 (MQ=255)
                                                                                                                              tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTg                         <  1:1522957/150‑1 (MQ=255)
                                                                                                                                 cTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGttt                      >  2:133461/1‑150 (MQ=255)
                                                                                                                                 cTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGttt                      >  2:599859/1‑150 (MQ=255)
                                                                                                                                     gCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTcc                  <  1:1444826/150‑1 (MQ=255)
                                                                                                                                      cTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATGTCCCCTGTTTTTCCt                 >  2:1321807/1‑150 (MQ=255)
                                                                                                                                          gaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGGTTTTCCTGAAt             <  1:1217870/150‑1 (MQ=255)
                                                                                                                                           ctctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATc            >  2:1119649/4‑150 (MQ=255)
                                                                                                                                           ctctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATc            >  2:1227946/4‑150 (MQ=255)
                                                                                                                                           aaTTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAACGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATc            <  1:409457/150‑1 (MQ=255)
                                                                                                                                                  ccTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTg     <  2:204388/150‑1 (MQ=255)
                                                                                                                                                   cTAAAGGGAAAGAGTGTATGCATGTGGAAAAATCGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGa    >  2:1234787/1‑150 (MQ=255)
                                                                                                                                                     aaaGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACCCTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATg  >  2:450698/1‑150 (MQ=255)
                                                                                                                                                     aaaGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATg  <  1:481475/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTGATG  >  NC_000913/344555‑344853

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: