Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,386,126 C→T 6.0% R1082H (CGT→CAT)  mscM ← miniconductance mechanosensitive channel MscM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,386,1260CT 6.0% 754.4 / 21.9 233R1082H (CGT→CAT) mscMminiconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (109/110);  new base T (7/7);  total (116/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 240 total aligned reads displayed.

TGGCTGTGGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGA  >  NC_000913/4385977‑4386275
                                                                                                                                                     |                                                                                                                                                     
        gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTgg                                                                                                                                               <  2:727813/150‑1 (MQ=255)
        gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTgg                                                                                                                                               <  1:425365/150‑1 (MQ=255)
         cAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGa                                                                                                                                              >  1:393611/1‑150 (MQ=255)
           gTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAg                                                                                                                                            >  1:794763/1‑150 (MQ=255)
           gTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAg                                                                                                                                            >  1:818070/1‑150 (MQ=255)
            ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTTCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaag                                                                                                                                           <  2:698531/150‑3 (MQ=255)
            ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAgg                                                                                                                                           <  1:458207/150‑1 (MQ=255)
              ttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaagcc                                                                                                                                         <  1:577670/150‑5 (MQ=255)
              ttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaagcc                                                                                                                                         <  1:1931188/150‑5 (MQ=255)
               tttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCg                                                                                                                                        <  1:1724135/150‑1 (MQ=255)
                ttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCgg                                                                                                                                       <  1:1764705/150‑1 (MQ=255)
                 tCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCggg                                                                                                                                      >  2:712336/1‑150 (MQ=255)
                  cTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                                                     >  2:1015334/1‑150 (MQ=255)
                   tGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGaa                                                                                                                                    >  2:796130/1‑150 (MQ=255)
                     cgcTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATg                                                                                                                                  <  1:1294029/150‑1 (MQ=255)
                      gcTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                                                                 <  2:1276408/150‑1 (MQ=255)
                      gcTATCAGTTTTGTTTGTGAGCCGGATTGGGTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                                                                 <  1:1734345/150‑1 (MQ=255)
                            aGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGCCGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                           >  1:1743610/1‑150 (MQ=255)
                            aGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                           <  1:788743/150‑1 (MQ=255)
                             gTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                          <  1:1527573/150‑1 (MQ=255)
                                   tttgtGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCATTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCg                                                                                                                    >  2:1711582/1‑150 (MQ=255)
                                     tgtgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                                  >  1:534527/1‑150 (MQ=255)
                                     tgtgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                                  >  2:187474/1‑150 (MQ=255)
                                     tgtgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                                  >  1:1003525/1‑150 (MQ=255)
                                      gtgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTgg                                                                                                                 >  1:125667/1‑150 (MQ=255)
                                       tgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTggg                                                                                                                <  1:278464/150‑1 (MQ=255)
                                        gAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGCGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGc                                                                                                               >  1:1221528/1‑150 (MQ=255)
                                             ggATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGa                                                                                                          <  1:712336/150‑1 (MQ=255)
                                               aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAag                                                                                                        <  2:1312376/150‑1 (MQ=255)
                                                  ggTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                     >  2:1485898/1‑150 (MQ=255)
                                                  ggTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                     >  2:222085/1‑150 (MQ=255)
                                                  ggTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                     >  2:1404170/1‑150 (MQ=255)
                                                    ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccagc                                                                                                   >  1:1106197/1‑150 (MQ=255)
                                                    ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccagc                                                                                                   >  2:270449/1‑150 (MQ=255)
                                                           ggCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATc                                                                                            <  2:1155420/150‑1 (MQ=255)
                                                           ggCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGGCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATc                                                                                            <  2:1342262/150‑1 (MQ=255)
                                                                 caaacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtgg                                                                                      >  2:639556/1‑150 (MQ=255)
                                                                   aacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggtg                                                                                    <  1:286889/150‑1 (MQ=255)
                                                                   aacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggtg                                                                                    <  2:402470/150‑1 (MQ=255)
                                                                     cttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGa                                                                                  >  1:206513/1‑150 (MQ=255)
                                                                     cttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGa                                                                                  <  1:1015334/150‑1 (MQ=255)
                                                                       tacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtc                                                                                >  1:1507681/1‑150 (MQ=255)
                                                                        acaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagc                                                                                                       >  2:1313970/1‑126 (MQ=255)
                                                                        acaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtct                                                                               >  2:246506/1‑150 (MQ=255)
                                                                         caaacttCCCGCCTGACGACCTTTGCCCGCAGCCGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                              >  2:763097/1‑150 (MQ=255)
                                                                         caaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtgg                                                                                      >  1:169979/1‑142 (MQ=255)
                                                                         caaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                              >  2:322643/1‑150 (MQ=255)
                                                                         caaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                              >  2:483244/1‑150 (MQ=255)
                                                                          aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                             <  1:1777420/150‑1 (MQ=255)
                                                                            acttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                           <  2:203802/150‑1 (MQ=255)
                                                                                cccGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGt                                                                       <  2:665151/150‑1 (MQ=255)
                                                                                     cTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGCGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCgg                                                                  <  2:1221528/150‑1 (MQ=255)
                                                                                      tGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGc                                                                 <  2:666740/150‑1 (MQ=255)
                                                                                          gacCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATAc                                                             <  2:1003525/150‑1 (MQ=255)
                                                                                            ccTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                                                           >  2:462710/1‑150 (MQ=255)
                                                                                            ccTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                                                           >  1:1193465/1‑150 (MQ=255)
                                                                                             cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGt                                                          >  1:1130497/1‑150 (MQ=255)
                                                                                                 gCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATTTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGAcc                                                      >  2:560075/1‑150 (MQ=255)
                                                                                                 gCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGAcc                                                      <  1:479756/150‑1 (MQ=255)
                                                                                                   ccGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCa                                                    >  1:1310572/1‑150 (MQ=255)
                                                                                                   ccGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCa                                                    <  2:1452895/150‑1 (MQ=255)
                                                                                                     gCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACGGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCAtc                                                  <  1:1451532/150‑1 (MQ=255)
                                                                                                        gACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCa                                               >  1:1804727/1‑150 (MQ=255)
                                                                                                            tCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcgg                                           <  2:306237/150‑1 (MQ=255)
                                                                                                               gCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCGGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGt                                        >  2:1250934/1‑150 (MQ=255)
                                                                                                                    ctcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATa                                   >  2:1888525/1‑150 (MQ=255)
                                                                                                                     tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGGATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGc                                                                         >  1:1606614/1‑111 (MQ=255)
                                                                                                                     tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATAc                                  <  2:1867940/150‑1 (MQ=255)
                                                                                                                      ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcgg                                           >  2:1766621/1‑140 (MQ=255)
                                                                                                                      ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACg                                 >  2:381182/1‑150 (MQ=255)
                                                                                                                      ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACg                                 >  2:1743339/1‑150 (MQ=255)
                                                                                                                       cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGc                                <  1:187474/150‑1 (MQ=255)
                                                                                                                       cccGTTAGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGCCCTAGCGGCATACGGTGACCCATCTCAGCCGCGTAAATAccc                                >  2:1650679/1‑148 (MQ=255)
                                                                                                                         cGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAg                              <  2:527952/150‑1 (MQ=255)
                                                                                                                             tgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCTGCGTATTTACGCCGCtga                          <  1:1485898/150‑3 (MQ=255)
                                                                                                                             tgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCg                          >  2:928509/1‑150 (MQ=255)
                                                                                                                              gtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGa                         >  2:978678/1‑150 (MQ=255)
                                                                                                                               tttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaa                        >  2:617201/1‑150 (MQ=255)
                                                                                                                                ttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCTGCGTATTTACGCCGCtgagat                       <  2:39747/150‑6 (MQ=255)
                                                                                                                                ttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaaa                       <  1:1404170/150‑1 (MQ=255)
                                                                                                                                  aCCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAAt                     >  2:1560131/1‑150 (MQ=255)
                                                                                                                                     gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGc                                                                         >  1:898652/1‑95 (MQ=255)
                                                                                                                                     gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGc                                                                         <  2:898652/95‑1 (MQ=255)
                                                                                                                                     gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTg                  <  1:270449/150‑1 (MQ=255)
                                                                                                                                     gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTg                  >  1:1864780/1‑150 (MQ=255)
                                                                                                                                     gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTg                  <  1:1538637/150‑1 (MQ=255)
                                                                                                                                     gTTGAGGCTTTCCAGACGCAGCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAACCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTg                  <  2:1594192/150‑1 (MQ=255)
                                                                                                                                       tGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGaa                <  1:918693/150‑1 (MQ=255)
                                                                                                                                       tGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGaa                >  1:65546/1‑150 (MQ=255)
                                                                                                                                          ggCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTg                  >  2:420346/1‑145 (MQ=255)
                                                                                                                                          ggCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcc             >  1:1325406/1‑150 (MQ=255)
                                                                                                                                           cccttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATccc            >  1:1982306/4‑150 (MQ=255)
                                                                                                                                            ccttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcccc           >  1:346380/3‑150 (MQ=255)
                                                                                                                                            cTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcccc           <  2:1130497/150‑1 (MQ=255)
                                                                                                                                             cttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCa                                                                                           >  1:404888/2‑69 (MQ=39)
                                                                                                                                                  aGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGc     <  2:1507681/150‑1 (MQ=255)
                                                                                                                                                    aCGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAg   <  2:1701134/150‑1 (MQ=255)
                                                                                                                                                    aCGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAg   <  2:1957505/150‑1 (MQ=255)
                                                                                                                                                     cGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGa  >  2:1146228/1‑150 (MQ=255)
                                                                                                                                                     cGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGa  >  1:1067944/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCTGTGGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGA  >  NC_000913/4385977‑4386275

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: