Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,140 A→G 17.3% intergenic (‑191/+38) ibsC ← / ← serA toxic membrane protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1400AG17.3% 485.1 / 91.8 223intergenic (‑191/+38)ibsC/serAtoxic membrane protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (83/99);  new base G (16/22);  total (100/123)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 227 total aligned reads displayed.

CAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTT  >  NC_000913/3057041‑3057239
                                                                                                   |                                                                                                   
    gCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCGACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCa                                                                                                 >  1:2192499/1‑100 (MQ=255)
     cATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAg                                                                                                <  1:2653763/100‑1 (MQ=255)
        tGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGa                                                                                             <  2:318782/100‑1 (MQ=255)
         gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTTGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                            <  2:1110520/100‑1 (MQ=255)
          aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGAct                                                                                           <  2:2143307/100‑2 (MQ=255)
          aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                           <  1:3743551/100‑1 (MQ=255)
          aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                           <  2:3381081/100‑1 (MQ=255)
          aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                           <  1:2153570/100‑1 (MQ=255)
           tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGActt                                                                                          >  2:922075/1‑98 (MQ=255)
            cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCATCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTg                                                                                         >  2:1788396/1‑100 (MQ=255)
                cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                     <  1:4016214/100‑1 (MQ=255)
                cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                     >  2:3089868/1‑100 (MQ=255)
                 tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAGCTGCCCg                                                                                     >  2:3802788/1‑99 (MQ=255)
                   tGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttt                                                                                  <  1:15613/100‑2 (MQ=38)
                   tGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttt                                                                                  <  1:1467198/100‑2 (MQ=38)
                    gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                 <  1:2062643/100‑3 (MQ=38)
                    gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                 <  1:2163569/100‑3 (MQ=38)
                     atttatttTCTGGAGCAGGCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt                                                                                >  1:1875591/1‑100 (MQ=255)
                     atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCTCGttttt                                                                                <  2:2190823/100‑4 (MQ=38)
                     atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                                <  1:386399/100‑4 (MQ=38)
                     atttatttTCTGGAGCAGAATCGCAAAGTAGAACGCTCAACGCGGAAACGATGTGGAGACGTGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                 <  2:3944820/99‑1 (MQ=255)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAATCAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                               <  1:3614927/100‑5 (MQ=25)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                               <  2:3458708/100‑5 (MQ=38)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                               <  2:904040/100‑5 (MQ=38)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                               <  1:922075/100‑5 (MQ=38)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                               >  2:1373800/1‑100 (MQ=255)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                               >  2:3589794/1‑100 (MQ=255)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGAAAGTCAGTGACCTGCCCGTTGAtt                                                                               >  2:1138009/1‑100 (MQ=255)
                       ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAttt                                                                              <  2:1519917/100‑1 (MQ=255)
                        tatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttttta                                                                             <  2:1964092/100‑7 (MQ=37)
                         atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttttat                                                                            <  1:1620329/100‑8 (MQ=37)
                         atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                            >  1:991209/1‑100 (MQ=255)
                         atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                            <  2:2349674/100‑1 (MQ=255)
                         atttTCTGGAGCAGACTCGCAAAGTAGAAAGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                            <  2:206171/100‑1 (MQ=255)
                           tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                          <  1:759625/100‑1 (MQ=255)
                           tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                          <  1:3705642/100‑1 (MQ=255)
                           tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                          <  2:497971/100‑1 (MQ=255)
                            ttCTGGAGCAGACTCGCAAAGTAGAATGCGCATCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                         <  1:1788400/100‑1 (MQ=255)
                                ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAATTTCAGAGAAg                                                                     >  1:649840/1‑100 (MQ=255)
                                 gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                    <  1:195026/100‑1 (MQ=255)
                                   ggaGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                  <  1:4073625/98‑1 (MQ=255)
                                   gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                  >  1:1533379/1‑100 (MQ=255)
                                   gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                  <  1:2781690/100‑1 (MQ=255)
                                     aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                    >  2:2576548/1‑96 (MQ=255)
                                       aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                              <  2:646639/100‑1 (MQ=255)
                                       aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                              <  1:967809/100‑1 (MQ=255)
                                         tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAg                                                            >  2:1203064/1‑100 (MQ=255)
                                          cGCAAAGTAGAATGCTCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                           <  1:4123917/100‑1 (MQ=255)
                                          cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                           >  1:3085638/1‑100 (MQ=255)
                                           gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGTAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                          <  1:3154070/100‑1 (MQ=255)
                                           gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                          >  2:1809132/1‑100 (MQ=255)
                                            cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                         >  2:4110590/1‑100 (MQ=255)
                                             aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAca                                                        <  1:1203065/100‑1 (MQ=255)
                                              aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                       <  1:125288/100‑1 (MQ=255)
                                               aGTAGACTGCCCCACGCGGCAGCGTTGTGGAGAAGGGATAAAATAACGGGCAAGTGAATGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                       <  2:3756602/99‑1 (MQ=255)
                                                gTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagca                                                     >  2:46156/1‑100 (MQ=255)
                                                gTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagca                                                     <  1:3466643/100‑1 (MQ=255)
                                                  aGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGa                                                   >  2:1751798/1‑100 (MQ=255)
                                                   gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagcaact                                                  >  1:3349941/1‑97 (MQ=255)
                                                    aaTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACg                                                 >  1:420205/1‑100 (MQ=255)
                                                       gcgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                              >  1:672662/1‑100 (MQ=255)
                                                        cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                             >  1:3901111/1‑100 (MQ=255)
                                                         gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcgc                                            >  1:3411811/1‑100 (MQ=255)
                                                         gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGAGCAAGTCAGTGACCTGCCCGTTAATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcgc                                            >  2:843589/1‑100 (MQ=255)
                                                              gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATg                                       <  1:3298379/100‑1 (MQ=255)
                                                                 gCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                    <  2:2110167/100‑1 (MQ=255)
                                                                  cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTAAGAGAAGGGGAATTAGTACCGCAGACCGGCGCGAATGGTAc                                   >  2:3212446/1‑100 (MQ=255)
                                                                      ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGTCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                               <  1:2694191/100‑1 (MQ=255)
                                                                         gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                            <  1:1751802/100‑1 (MQ=255)
                                                                           ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                          >  1:3784791/1‑100 (MQ=255)
                                                                           ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                          >  1:224012/1‑100 (MQ=255)
                                                                            gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                         >  2:749399/1‑100 (MQ=255)
                                                                            gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                         <  1:492003/100‑1 (MQ=255)
                                                                             agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCGCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCt                        >  1:2951239/1‑100 (MQ=255)
                                                                             agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCt                        <  1:1373802/100‑1 (MQ=255)
                                                                              gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCtt                       >  1:3381870/1‑100 (MQ=255)
                                                                                aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                     <  1:2880666/100‑1 (MQ=255)
                                                                                aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                     <  1:2084180/100‑1 (MQ=255)
                                                                                   tgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                  >  1:598684/9‑100 (MQ=37)
                                                                                   gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                  >  1:862819/1‑100 (MQ=255)
                                                                                    aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTg                 <  1:3554942/100‑1 (MQ=255)
                                                                                      aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc               >  1:1063231/1‑100 (MQ=255)
                                                                                      aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc               >  2:3250306/1‑100 (MQ=255)
                                                                                      aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc               <  1:408961/100‑1 (MQ=255)
                                                                                      aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc               >  1:384076/1‑100 (MQ=255)
                                                                                         aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc            <  2:306425/100‑1 (MQ=255)
                                                                                         aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc            >  1:529749/1‑100 (MQ=255)
                                                                                         aaaaCGGGCAAGGCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGTAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc            <  2:1579215/100‑1 (MQ=255)
                                                                                            aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc         >  1:1000886/1‑100 (MQ=39)
                                                                                             cGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg        >  2:2655099/1‑100 (MQ=255)
                                                                                             cGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg        >  1:3204605/1‑100 (MQ=255)
                                                                                              gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc       <  2:3994223/100‑1 (MQ=255)
                                                                                              gggCAAGTCAGTGACCTGCCCGTTGAGTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc       <  2:3642211/100‑1 (MQ=255)
                                                                                               ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt      >  2:3412722/1‑100 (MQ=255)
                                                                                               ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt      <  2:2707134/100‑1 (MQ=255)
                                                                                               ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt      <  1:3188911/100‑1 (MQ=255)
                                                                                                gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGCATGAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt     <  2:2943608/100‑1 (MQ=255)
                                                                                                gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt     <  2:2972791/100‑1 (MQ=255)
                                                                                                  aaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGCATTAGTACAGCAGACGGGCGCGAATCGTACCCGGAATAGCTTTCATTGCCTGCAgcgcgt     >  2:165063/1‑96 (MQ=255)
                                                                                                  aaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt   >  2:405558/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTT  >  NC_000913/3057041‑3057239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: