Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,102,948 T→C 5.6% K83K (AAA→AAG) ‡ csgD ← DNA‑binding transcriptional dual regulator CsgD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,9480TC 5.6% 324.9 / 7.0 106K83K (AAA→AAG) ‡csgDDNA‑binding transcriptional dual regulator CsgD
Reads supporting (aligned to +/- strand):  ref base T (55/45);  new base C (2/4);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

Only 100 of 109 total aligned reads displayed.

TGCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAG  >  NC_000913/1102799‑1103097
                                                                                                                                                     |                                                                                                                                                     
tGCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATAtt                                                                                                                                                       >  1:913312/1‑150 (MQ=255)
      aaGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGata                                                                                                                                                 >  2:159843/1‑150 (MQ=255)
       aGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGatat                                                                                                                                                <  1:1426025/150‑1 (MQ=255)
        gTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAACCCCCGTTCAAAACACCCCAGTTTTCATTGTCCCTGAACGGTTAATATACCGGCGAATTTAGAAAAAAAATTTTTATAtt                                                                                                                                               >  1:430089/1‑150 (MQ=255)
         tAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACATCACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACTGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttg                                                                                                                                              <  2:1503477/150‑1 (MQ=255)
          aGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgt                                                                                                                                             <  2:346372/150‑1 (MQ=255)
          aGCAATCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATCGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgt                                                                                                                                             >  2:1176150/1‑150 (MQ=255)
           gCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgtt                                                                                                                                            >  2:890851/1‑150 (MQ=255)
            cATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgttg                                                                                                                                           <  2:90993/150‑1 (MQ=255)
            cATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCAGGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgttg                                                                                                                                           <  2:1131893/150‑1 (MQ=255)
             aTTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATtgttgt                                                                                                                                          >  1:1290020/1‑150 (MQ=255)
             aTTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgc                                                                                                                                          <  1:1479516/150‑2 (MQ=255)
              ttCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTgttgtt                                                                                                                                         >  2:139101/1‑150 (MQ=255)
              ttCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgct                                                                                                                                         <  1:1096914/150‑3 (MQ=255)
               tCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgcta                                                                                                                                        <  2:1330205/150‑4 (MQ=255)
               tCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgcta                                                                                                                                        <  1:1495096/150‑4 (MQ=255)
                cgccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGtttt                                                                                                                                       <  2:422830/150‑1 (MQ=255)
                cgccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGtttt                                                                                                                                       >  1:453821/1‑150 (MQ=255)
                cgccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGtttt                                                                                                                                       >  2:1365140/1‑150 (MQ=255)
                   cgcgCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTcc                                                                                                                                    <  1:281555/150‑1 (MQ=255)
                   cgcgCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTcc                                                                                                                                    <  2:1114824/150‑1 (MQ=255)
                    gcgcAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATTTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCt                                                                                                                                   <  2:212412/150‑1 (MQ=255)
                      gcAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGc                                                                                                                                 >  1:569858/1‑150 (MQ=255)
                       cAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCt                                                                                                                                <  1:930200/150‑1 (MQ=255)
                       cAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCt                                                                                                                                <  1:924359/150‑1 (MQ=255)
                       cAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCt                                                                                                                                >  2:1382191/1‑150 (MQ=255)
                         ggACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCa                                                                                                                              >  1:262306/1‑150 (MQ=255)
                           aCGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCaaa                                                                                                                            >  1:225056/1‑150 (MQ=255)
                             gCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGt                                                                                                                          <  2:633114/150‑1 (MQ=255)
                               cTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTAt                                                                                                                        >  2:1291169/1‑150 (MQ=255)
                               cTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTAt                                                                                                                        >  1:1180432/1‑150 (MQ=255)
                                 tGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATcc                                                                                                                      <  1:778291/150‑1 (MQ=255)
                                  gCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCt                                                                                                                     >  1:53273/1‑150 (MQ=255)
                                   cAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTg                                                                                                                    <  1:1291169/150‑1 (MQ=255)
                                     aCCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGcc                                                                                                                  <  1:233129/150‑1 (MQ=255)
                                     aCCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGcc                                                                                                                  >  1:678100/1‑150 (MQ=255)
                                     aCCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTATAAGTATCCTGcc                                                                                                                  <  2:290064/150‑1 (MQ=255)
                                        cATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCaat                                                                                                               >  1:162390/1‑150 (MQ=255)
                                           tGACAACACGTTCTTGATCCTCCATTGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCaataat                                                                                                            <  2:242855/150‑1 (MQ=255)
                                             acaacaCGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATgg                                                                                                          <  1:1159020/150‑1 (MQ=255)
                                               aacaCGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGAt                                                                                                        >  2:396013/1‑150 (MQ=255)
                                                acacGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTg                                                                                                                    >  2:697245/1‑137 (MQ=255)
                                                acacGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATa                                                                                                       >  2:1138716/1‑150 (MQ=255)
                                                acacGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATa                                                                                                       >  2:902662/1‑150 (MQ=255)
                                                   cGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAga                                                                                                    >  2:1213559/1‑149 (MQ=255)
                                                         tGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTa                                                                                              >  1:631797/1‑150 (MQ=255)
                                                          gATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTAt                                                                                             <  1:1138716/150‑1 (MQ=255)
                                                          gATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTAt                                                                                             >  2:1242467/1‑150 (MQ=255)
                                                            tcctccATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATcc                                                                                           >  1:234912/1‑150 (MQ=255)
                                                             cctccATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCg                                                                                          <  1:139101/150‑1 (MQ=255)
                                                             cctccATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATACTGCCAATAATGGATAAGCTTTTTATCCg                                                                                          >  1:250261/1‑150 (MQ=255)
                                                              ctccATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGc                                                                                         >  2:726924/1‑150 (MQ=255)
                                                               tccATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCt                                                                                        >  1:1518039/1‑150 (MQ=255)
                                                                 cATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTc                                                                                      >  2:1033345/1‑150 (MQ=255)
                                                                   tGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCca                                                                                    <  1:1382191/150‑1 (MQ=255)
                                                                    ggAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcat                                                                                   >  2:1411089/1‑150 (MQ=255)
                                                                       aTAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATACGCTTCcatcat                                                                                >  1:541118/1‑150 (MQ=255)
                                                                        tAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCata                                                                               >  2:1431562/1‑150 (MQ=255)
                                                                          aaaaaCGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATc                                                                             >  1:915950/1‑150 (MQ=255)
                                                                            aaaCGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCa                                                                           <  2:1186775/150‑1 (MQ=255)
                                                                               cGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaa                                                                        <  2:569858/150‑1 (MQ=255)
                                                                                  cGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAgaagaa                                                                     >  2:619703/1‑150 (MQ=255)
                                                                                   gTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaat                                                                    <  1:396013/150‑1 (MQ=255)
                                                                                    ttGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaata                                                                   <  2:1180432/150‑1 (MQ=255)
                                                                                    ttGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCAGCTTCCATCATATCCAGAAGaata                                                                   >  1:852452/1‑150 (MQ=255)
                                                                                      tatatGAGGACAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTAATTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTAAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaataat                                                                 <  2:678100/149‑1 (MQ=255)
                                                                                        tatGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATag                                                               <  1:902662/150‑1 (MQ=255)
                                                                                         atGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATaga                                                              >  2:1234021/1‑150 (MQ=255)
                                                                                          tGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATagag                                                             >  1:327886/1‑150 (MQ=255)
                                                                                            aGGCCAGTTTTCAATGTCGCGGTACGGTTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGcc                                                           <  1:1365140/150‑1 (MQ=255)
                                                                                             ggCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCt                                                          >  1:51603/1‑150 (MQ=255)
                                                                                               ccAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGa                                                        >  1:939237/1‑150 (MQ=255)
                                                                                               ccAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGa                                                        <  1:1333477/150‑1 (MQ=255)
                                                                                                  gTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAaga                                                     >  1:1132064/1‑150 (MQ=255)
                                                                                                  gTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAaga                                                     <  2:234912/150‑1 (MQ=255)
                                                                                                     ttCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGat                                                  <  2:162390/150‑1 (MQ=255)
                                                                                                       cgaTGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGatat                                                >  2:1454973/3‑150 (MQ=255)
                                                                                                           gTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTc                                            <  1:1108596/150‑1 (MQ=255)
                                                                                                            tCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATAATTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTcc                                           <  2:1110716/150‑1 (MQ=255)
                                                                                                             cgcgGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATACAGAAGAATAATAGAGCCTGAAGAGATATCGTCCa                                          >  2:204632/1‑150 (MQ=255)
                                                                                                               cgGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCaga                                        <  1:619703/150‑1 (MQ=255)
                                                                                                                ggTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTTATATTTTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCagag                                       <  1:1046972/150‑1 (MQ=255)
                                                                                                                      ggTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGtt                                 >  1:441040/1‑150 (MQ=255)
                                                                                                                          aTCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAat                             <  1:1411089/150‑1 (MQ=255)
                                                                                                                           tCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAata                            >  1:794311/1‑150 (MQ=255)
                                                                                                                                  gCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTa                     >  1:1480332/1‑150 (MQ=255)
                                                                                                                                     tATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAAtt                  >  2:1125060/1‑150 (MQ=255)
                                                                                                                                      aTTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAAttt                 >  2:1026026/1‑150 (MQ=255)
                                                                                                                                        aaaacaacaaTATCAAAATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTc               >  1:1017076/5‑150 (MQ=255)
                                                                                                                                         tAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTcc              >  2:1329231/1‑150 (MQ=255)
                                                                                                                                         tAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTcc              <  1:159843/150‑1 (MQ=255)
                                                                                                                                         aaacaacaaTATCAAAATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTcc              >  1:347229/4‑150 (MQ=255)
                                                                                                                                           gCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTg            <  1:1234021/150‑1 (MQ=255)
                                                                                                                                           gCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTg            <  1:1125060/150‑1 (MQ=255)
                                                                                                                                           gCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTg            <  2:1480332/150‑1 (MQ=255)
                                                                                                                                             aacaaTATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTa          >  1:1101138/1‑150 (MQ=255)
                                                                                                                                                 atatTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTg      >  1:1130111/1‑150 (MQ=255)
                                                                                                                                                   atTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGcc    <  2:852452/150‑1 (MQ=255)
                                                                                                                                                    ttttGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCa   >  2:89485/1‑150 (MQ=255)
                                                                                                                                                     tttGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAg  <  2:1061523/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAG  >  NC_000913/1102799‑1103097

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