Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,102,949 T→A 5.7% K83I (AAA→ATA) ‡ 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,9490TA 5.7% 334.4 / 7.4 106K83I (AAA→ATA) ‡csgDDNA‑binding transcriptional dual regulator CsgD
Reads supporting (aligned to +/- strand):  ref base T (55/45);  new base A (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 = 6.36e-01

Only 100 of 108 total aligned reads displayed.

AAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAG  >  NC_000913/1102805‑1103097
                                                                                                                                                |                                                                                                                                                    
aaGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGata                                                                                                                                                 >  2:159843/1‑150 (MQ=255)
aaGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGata                                                                                                                                                 <  2:1055826/150‑1 (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)
       aTTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgc                                                                                                                                          <  1:1479516/150‑2 (MQ=255)
        ttCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTgttgtt                                                                                                                                         >  2:139101/1‑150 (MQ=255)
        ttCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgct                                                                                                                                         <  1:1096914/150‑3 (MQ=255)
         tCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgcta                                                                                                                                        <  1:1495096/150‑4 (MQ=255)
         tCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgcta                                                                                                                                        <  2:1330205/150‑4 (MQ=255)
          cgccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGtttt                                                                                                                                       >  1:1284507/1‑150 (MQ=255)
          cgccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGtttt                                                                                                                                       >  2:1365140/1‑150 (MQ=255)
             cgcgCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTcc                                                                                                                                    <  2:1114824/150‑1 (MQ=255)
             cgcgCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTcc                                                                                                                                    >  2:1240108/1‑150 (MQ=255)
             cgcgCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTcc                                                                                                                                    <  1:281555/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                                                                                                                        >  1:1180432/1‑150 (MQ=255)
                         cTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTAt                                                                                                                        >  2:1291169/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)
                               aCCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGcc                                                                                                                  >  1:359232/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                                                                                                                    <  1:697245/137‑1 (MQ=255)
                                          acacGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTg                                                                                                                    >  2:697245/1‑137 (MQ=255)
                                          acacGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATa                                                                                                       >  2:902662/1‑150 (MQ=255)
                                          acacGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATa                                                                                                       >  2:1138716/1‑150 (MQ=255)
                                             cGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAga                                                                                                    >  2:1213559/1‑149 (MQ=255)
                                                   tGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTa                                                                                              >  1:631797/1‑150 (MQ=255)
                                                    gATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTAt                                                                                             >  2:1242467/1‑150 (MQ=255)
                                                    gATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTAt                                                                                             <  1:1138716/150‑1 (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)
                                                                      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)
                                                                              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                                                             >  2:1499651/1‑150 (MQ=255)
                                                                                    tGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATagag                                                             >  1:327886/1‑150 (MQ=255)
                                                                                      agtccAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGTGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGcc                                                           <  2:430089/147‑1 (MQ=255)
                                                                                      aGGCCAGTTTTCAATGTCGCGGTACGGTTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGcc                                                           <  1:1365140/150‑1 (MQ=255)
                                                                                       ggCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCt                                                          >  1:51603/1‑150 (MQ=255)
                                                                                         ccAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGa                                                        <  1:1333477/150‑1 (MQ=255)
                                                                                         ccAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGa                                                        >  1:939237/1‑150 (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)
                                                                                                                                   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:1510987/1‑150 (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)
                                                                                                                                                |                                                                                                                                                    
AAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAG  >  NC_000913/1102805‑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: