Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,483,342 G→C 6.3% Q94H (CAG→CAC hrpA → ATP‑dependent RNA helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,483,3420GC 6.3% 537.9 / 19.5 158Q94H (CAG→CAChrpAATP‑dependent RNA helicase HrpA
Reads supporting (aligned to +/- strand):  ref base G (77/71);  new base C (7/3);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 168 total aligned reads displayed.

GCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGA  >  NC_000913/1483193‑1483491
                                                                                                                                                     |                                                                                                                                                     
gCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATAACCAg                                                                                                                                                       >  1:1589336/1‑150 (MQ=255)
gCCATTTTCCAGGAGATCGCGAAAGAGATTCACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTAAGTCAGAAAAAACAGGACATACTCCAAGCGATACGTGATCACCAg                                                                                                                                                       >  2:1925926/1‑150 (MQ=255)
  caatttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGAGATTCGTGATCACCAggt                                                                                                                                                     >  2:411651/4‑150 (MQ=255)
  cATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACTGGACATTCTCGAAGCGATTCGTGATCACCAggt                                                                                                                                                     <  1:94230/150‑1 (MQ=255)
     tttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt                                                                                                                                                  <  1:307760/150‑1 (MQ=255)
     tttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt                                                                                                                                                  <  2:362335/150‑1 (MQ=255)
       tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                                                >  2:1317337/1‑150 (MQ=255)
         cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                                                              >  2:1130584/1‑150 (MQ=255)
          aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcg                                                                                                                                             <  2:769206/150‑1 (MQ=255)
           ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAATTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCACCTGGTGAtcac                                                                                                                                            <  2:1098377/150‑3 (MQ=255)
           ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCACCTGGTGAtcac                                                                                                                                            <  1:1313192/150‑3 (MQ=255)
              gaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGc                                                                                                                                         >  1:1852974/1‑150 (MQ=255)
               aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGcc                                                                                                                                        >  2:1525582/1‑150 (MQ=255)
                         agaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCGGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACgg                                                                                                                              >  1:2316113/1‑150 (MQ=255)
                              gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCt                                                                                                                         >  1:493472/1‑150 (MQ=255)
                                cctggcggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                                                       >  1:1532137/4‑150 (MQ=255)
                                ccAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                                                       <  2:1286762/150‑1 (MQ=255)
                                ccAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                                                       <  1:1240899/150‑1 (MQ=255)
                                      ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAAc                                                                                                                 <  2:2207702/150‑1 (MQ=255)
                                      ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAAc                                                                                                                 <  1:1181019/150‑1 (MQ=255)
                                       gcAGGGAAAGTCCTGCTGCGTGACGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACg                                                                                                                >  1:1655018/1‑150 (MQ=255)
                                       gcAGGGAAAGTCCTGCTGCGTGAAGCTGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACATGACATTCTCGAAGCGATTCGTGATCAGCAGGTGGTGATCGTCGCCGGGGAAACGTGTTCTGGTAAAACg                                                                                                                <  2:712783/150‑1 (MQ=255)
                                       gcAGGGAAAGTCCTGCTGCGGGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACg                                                                                                                <  2:2221628/150‑1 (MQ=255)
                                        cAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGa                                                                                                               >  1:2175302/1‑150 (MQ=255)
                                         aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGAc                                                                                                              <  1:179339/150‑1 (MQ=255)
                                         aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGAc                                                                                                              <  2:13481/150‑1 (MQ=255)
                                          gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACt                                                                                                             >  2:1792442/1‑150 (MQ=255)
                                          gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACt                                                                                                             <  1:578200/150‑1 (MQ=255)
                                          gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACt                                                                                                             <  1:1011937/150‑1 (MQ=255)
                                           ggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTc                                                                                                            <  1:2462706/150‑1 (MQ=255)
                                           ggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTc                                                                                                            <  1:1623913/150‑1 (MQ=255)
                                            gAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGCCTCa                                                                                                           >  2:2470179/1‑150 (MQ=255)
                                            gAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCa                                                                                                           <  2:2251359/150‑1 (MQ=255)
                                              aaGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGt                                                                                                         >  2:1610287/1‑150 (MQ=255)
                                                gTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTa                                                                                                       >  1:1093405/1‑150 (MQ=255)
                                                  ccTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAcc                                                                                                     <  1:1284829/150‑1 (MQ=255)
                                                   ctgctgCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCg                                                                                                    >  1:2403193/1‑150 (MQ=255)
                                                     gctgcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaa                                                                                                  >  1:249667/1‑150 (MQ=255)
                                                              aGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                                                                         >  2:1493763/1‑150 (MQ=255)
                                                                     cGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTg                                                                                  >  2:753307/1‑150 (MQ=255)
                                                                       aCCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTggg                                                                                >  2:1482760/1‑150 (MQ=255)
                                                                       aCCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTggg                                                                                >  2:2009983/1‑150 (MQ=255)
                                                                          ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                                                             <  1:1471514/150‑1 (MQ=255)
                                                                          ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                                                             >  2:1150695/1‑150 (MQ=255)
                                                                          ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                                                             >  1:1619218/1‑150 (MQ=255)
                                                                          ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                                                             >  1:908484/1‑150 (MQ=255)
                                                                           gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGAGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                                                            >  1:2417082/1‑150 (MQ=255)
                                                                             aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCgg                                                                          <  1:1770762/150‑1 (MQ=255)
                                                                              aTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCggg                                                                         <  2:1852974/150‑1 (MQ=255)
                                                                                     aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                                                  <  1:2484283/150‑1 (MQ=255)
                                                                                       ccTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGGAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGa                                                                >  2:2484384/1‑150 (MQ=255)
                                                                                        cTGACAATTTACCGGTTAGTCAGAAAAAACAGGTCATTCTCGAAGCGATTCGTGATCACCAGTTGGTGATCGTCGCCGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTACCGATAATCTGTATGGAGCTGGGGCGCTGGATTAAAGGag                                                               <  2:1398144/150‑2 (MQ=255)
                                                                                        cTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc                                                               <  1:710789/150‑1 (MQ=255)
                                                                                        cTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc                                                               <  2:2301904/150‑1 (MQ=255)
                                                                                          gACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg                                                             <  1:1493763/150‑1 (MQ=255)
                                                                                                ttACCGGTTAGTCCGAAAAAACAGGACATTCTCGACGCGCTGCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGGTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGc                                                       <  2:1610651/150‑1 (MQ=255)
                                                                                                ttACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGc                                                       <  2:2454011/150‑1 (MQ=255)
                                                                                                 tACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGcc                                                      >  1:2108344/1‑150 (MQ=255)
                                                                                                     ggTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                                                  >  1:1105558/1‑150 (MQ=255)
                                                                                                     ggTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                                                  <  2:14229/150‑1 (MQ=255)
                                                                                                      gTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAcc                                                 <  2:908484/150‑1 (MQ=255)
                                                                                                         aGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAg                                              >  2:2121556/1‑150 (MQ=255)
                                                                                                         aGTCAGAAAAAACAGGAAATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAg                                              <  2:1024944/150‑1 (MQ=255)
                                                                                                          gTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGc                                             >  1:647767/1‑150 (MQ=255)
                                                                                                          gTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGc                                             >  2:383521/1‑150 (MQ=255)
                                                                                                           tCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGcc                                            >  2:1188069/1‑150 (MQ=255)
                                                                                                            cAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg                                           <  2:288561/150‑1 (MQ=255)
                                                                                                            cAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg                                           >  2:85810/1‑150 (MQ=255)
                                                                                                            cAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg                                           >  1:1445442/1‑150 (MQ=255)
                                                                                                              gAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgcg                                         >  1:857609/1‑150 (MQ=255)
                                                                                                              caaaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgcg                                         <  2:583874/149‑1 (MQ=255)
                                                                                                                 aaaaCAGGACATTCTCTAAGCGATTCGTGATCACCAGGTGGTGATCGGCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGTAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtcg                                      <  2:1369573/150‑1 (MQ=255)
                                                                                                                  aaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtcgt                                     >  2:2386119/1‑150 (MQ=255)
                                                                                                                       ggACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggc                                >  1:1127011/1‑150 (MQ=255)
                                                                                                                       ggACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggc                                <  2:267033/150‑1 (MQ=255)
                                                                                                                       ggACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggc                                >  1:1545028/1‑150 (MQ=255)
                                                                                                                       ggACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGAGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggc                                >  1:133145/1‑150 (MQ=255)
                                                                                                                          cATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggcggc                             <  1:746503/150‑1 (MQ=255)
                                                                                                                          cATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggcggc                             >  2:2212425/1‑150 (MQ=255)
                                                                                                                          cATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATAGTCCCCGGGGAAACGGGTTCAGGTAAAACCACTCAGTTACAGAAAATCTGTATGGAGCTCGGGCGCGGGATTAAAGGACTCATCGCCCATCCCCACCCCAGTCGACAggcggc                             >  2:128494/1‑150 (MQ=255)
                                                                                                                            ttCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGt                                                                                                         <  1:1075069/72‑1 (MQ=255)
                                                                                                                            ttCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGt                                                                                                         >  2:1075069/1‑72 (MQ=255)
                                                                                                                               tcGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGaa                        >  1:422742/1‑150 (MQ=255)
                                                                                                                                     cGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTgg                  <  1:1379259/150‑1 (MQ=255)
                                                                                                                                      gATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGc                 >  1:1297507/1‑150 (MQ=255)
                                                                                                                                         tCGTGATGACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCTGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGaa              <  2:1761788/150‑1 (MQ=255)
                                                                                                                                          gacgATCACCACCTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAAc             >  2:1458204/4‑150 (MQ=255)
                                                                                                                                           acgATCACCACCTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAAcc            >  1:536215/3‑150 (MQ=255)
                                                                                                                                            tGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCg           >  1:255633/1‑150 (MQ=255)
                                                                                                                                            tGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCg           <  2:1868897/150‑1 (MQ=255)
                                                                                                                                            cgATCACCACCTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCg           >  2:2060656/2‑150 (MQ=255)
                                                                                                                                             gATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGt          >  1:1774828/1‑150 (MQ=255)
                                                                                                                                             gATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGt          >  1:353114/1‑150 (MQ=255)
                                                                                                                                             gATCACCACCTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGt          >  1:631147/1‑150 (MQ=255)
                                                                                                                                               tCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTAt        >  2:1350402/1‑150 (MQ=255)
                                                                                                                                                cACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTAtt       >  2:212308/1‑150 (MQ=255)
                                                                                                                                                  ccAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGc     <  2:493472/150‑1 (MQ=255)
                                                                                                                                                   cAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCg    >  2:1719320/1‑150 (MQ=255)
                                                                                                                                                     ggtggtGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGa  <  2:2175302/150‑1 (MQ=255)
                                                                                                                                                     ggtggtGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGa  >  2:2276487/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGA  >  NC_000913/1483193‑1483491

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: