Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,541 G→C 10.3% P28R (CCG→CGG)  yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,5410GC10.3% 394.2 / 31.3 117P28R (CCG→CGG) yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (60/45);  new base C (5/7);  total (65/52)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 119 total aligned reads displayed.

CTGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCG  >  NC_000913/3101393‑3101687
                                                                                                                                                    |                                                                                                                                                  
cTGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCAcgcg                                                                                                                                                    >  3:394457/1‑149 (MQ=255)
       tCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGa                                                                                                                                              <  2:743874/148‑1 (MQ=255)
       tCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGAc                                                                                                                                             >  4:346133/1‑149 (MQ=255)
        ccTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGAc                                                                                                                                             >  1:725510/1‑148 (MQ=255)
         cTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                           >  3:307898/1‑149 (MQ=255)
          ttCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                          <  4:239695/149‑3 (MQ=255)
           tCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                          <  2:957203/148‑3 (MQ=255)
            cGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTg                                                                                                                                         <  1:949758/148‑1 (MQ=255)
            cGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACGATgt                                                                                                                                        <  3:158520/149‑2 (MQ=255)
             gCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCa                                                                                                                                       >  4:184800/1‑149 (MQ=255)
              ccTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACgatgtg                                                                                                                                       <  2:51034/148‑6 (MQ=255)
               cTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCAc                                                                                                                                      >  1:687531/1‑148 (MQ=255)
               cTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACt                                                                                                                                     >  4:269292/1‑149 (MQ=255)
                 gCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTg                                                                                                                                    <  2:891566/148‑1 (MQ=255)
                             cggcATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTg                                                                                                                        >  2:39512/1‑148 (MQ=255)
                                cATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCATCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGc                                                                                                                     <  1:169753/148‑1 (MQ=255)
                                 aTTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                                   >  4:909465/1‑149 (MQ=255)
                                   tCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGt                                                                                                                  >  2:874130/1‑148 (MQ=255)
                                   tCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGt                                                                                                                  >  1:943628/1‑148 (MQ=255)
                                      gggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACga                                                                                                                                           >  4:852085/1‑119 (MQ=255)
                                       ggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAtt                                                                                                              >  2:565347/1‑148 (MQ=255)
                                       ggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCAtt                                                                                                              >  1:311235/1‑148 (MQ=255)
                                       ggAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTg                                                                                                             >  4:371386/1‑149 (MQ=255)
                                         aaTATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                          <  4:805041/118‑3 (MQ=255)
                                         aaTATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACgat                                                                                                                                          >  3:805041/1‑116 (MQ=255)
                                           tatTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGAGTGCCGCTTTCCCCCATCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGTGAGTGGGCGTGAGCGGTCATTGCCGt                                                                                                         <  4:565039/149‑1 (MQ=255)
                                              tGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCAc                                                                                                      <  4:107629/149‑1 (MQ=255)
                                                  tACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCa                                                                                                                                       >  2:305458/1‑112 (MQ=255)
                                                  tACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGa                                                                                                   >  2:448610/1‑148 (MQ=255)
                                                   aCCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                 <  3:534758/149‑1 (MQ=255)
                                                   aCCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                 <  3:184800/149‑1 (MQ=255)
                                                   aCCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTCCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGTGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                 <  3:264603/149‑1 (MQ=255)
                                                     cgtaacaCATCACGTTGCCGTCTGGCTTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGATCGGTTTGCGGGCTGACAAGCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGt                                                                                               >  4:814492/7‑149 (MQ=255)
                                                          taCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgctg                                                                                           <  1:118494/148‑1 (MQ=255)
                                                           aCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCg                                                                                                                    <  3:36990/122‑1 (MQ=255)
                                                           aCATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCg                                                                                                                    >  4:36990/1‑122 (MQ=255)
                                                             aTCACGTTGCCGTCTGTCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgc                                                                                       >  3:913222/1‑149 (MQ=255)
                                                             aTCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccg                                                                                        >  2:260813/1‑148 (MQ=255)
                                                              tCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgcc                                                                                      >  3:857126/1‑149 (MQ=255)
                                                               cACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgcc                                                                                      >  1:438425/1‑148 (MQ=255)
                                                                aCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcag                                                                                    <  4:546724/149‑1 (MQ=255)
                                                                   ttGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagca                                                                                  >  2:795884/1‑148 (MQ=255)
                                                                       cGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTg                                                                                                                        <  1:846063/106‑1 (MQ=255)
                                                                       cGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTg                                                                                                                        >  2:846063/1‑106 (MQ=255)
                                                                        gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                            >  4:128752/1‑149 (MQ=255)
                                                                        gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                            >  3:674083/1‑149 (MQ=255)
                                                                        gTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAtt                                                                             <  1:39512/148‑1 (MQ=255)
                                                                         tCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                            <  1:260813/148‑1 (MQ=255)
                                                                               gTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcat                                                                     >  3:195900/1‑149 (MQ=255)
                                                                                   tCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCACCAGCAGCATTTTGCGCATCATa                                                                 <  4:548434/149‑1 (MQ=255)
                                                                                     aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAg                                                                >  1:129728/1‑148 (MQ=255)
                                                                                     aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGCCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                               >  3:654359/1‑149 (MQ=255)
                                                                                     aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGACGGCGCGAGCGGTCATTGCCGTCACTGACAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                               >  3:422784/1‑149 (MQ=255)
                                                                                      ccAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                               >  2:199506/1‑148 (MQ=255)
                                                                                       cAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTc                                                              >  1:686901/1‑148 (MQ=255)
                                                                                       cAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCt                                                             >  3:2482/1‑149 (MQ=255)
                                                                                           ttGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTcc                                                          >  2:886764/1‑148 (MQ=255)
                                                                                            tGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTccc                                                         >  2:408826/1‑148 (MQ=255)
                                                                                             gCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTc                                                       >  4:109473/1‑149 (MQ=255)
                                                                                              ccGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTc                                                       >  2:631510/1‑148 (MQ=255)
                                                                                               cGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTACGCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCa                                                      <  1:698073/148‑1 (MQ=255)
                                                                                                      gCCCCTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                              >  3:760606/1‑149 (MQ=255)
                                                                                                       cccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAt                                             >  4:294449/1‑149 (MQ=255)
                                                                                                       cccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAt                                             >  3:746265/1‑149 (MQ=255)
                                                                                                        ccTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAt                                             >  2:744071/1‑148 (MQ=255)
                                                                                                          ttCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCa                                           >  1:558046/1‑148 (MQ=255)
                                                                                                           tCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatcat                                                                  >  1:251451/1‑124 (MQ=255)
                                                                                                                  cACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTc                                                                                                                              >  3:907046/1‑57 (MQ=255)
                                                                                                                  cACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTc                                                                                                                              <  4:907046/57‑1 (MQ=255)
                                                                                                                           cGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGc                          <  2:765513/148‑1 (MQ=255)
                                                                                                                            gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                        <  3:294449/149‑1 (MQ=255)
                                                                                                                            gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcc                         <  2:275212/148‑1 (MQ=255)
                                                                                                                            gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcc                         <  1:886764/148‑1 (MQ=255)
                                                                                                                                tGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                    >  3:47255/1‑149 (MQ=255)
                                                                                                                                  aCAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                  >  4:720741/1‑149 (MQ=255)
                                                                                                                                  aCAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgatttt                  >  3:695732/1‑142 (MQ=255)
                                                                                                                                    aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                <  4:148068/149‑1 (MQ=255)
                                                                                                                                     aTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                <  1:734926/148‑1 (MQ=255)
                                                                                                                                     aTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCAACATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                <  2:725510/148‑1 (MQ=255)
                                                                                                                                      tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCa                                           >  1:471836/1‑120 (MQ=255)
                                                                                                                                      tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCa                                           <  2:471836/120‑1 (MQ=255)
                                                                                                                                      tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt               <  4:423856/148‑1 (MQ=255)
                                                                                                                                      tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt               >  3:423856/1‑148 (MQ=255)
                                                                                                                                      tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCAATGAAAGTGCTGACGCGAGAAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt               <  1:408826/148‑1 (MQ=255)
                                                                                                                                       cacaTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa             >  4:318041/1‑149 (MQ=255)
                                                                                                                                          agcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa          >  3:639358/3‑149 (MQ=255)
                                                                                                                                          agcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa          >  3:506978/3‑149 (MQ=255)
                                                                                                                                           tCGTCACGCGGCGTGATGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa         <  3:814492/149‑1 (MQ=255)
                                                                                                                                           tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa         <  3:640261/149‑1 (MQ=255)
                                                                                                                                           tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa         >  4:649472/1‑149 (MQ=255)
                                                                                                                                           tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa          >  2:323710/1‑148 (MQ=255)
                                                                                                                                            cGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                       <  1:139005/118‑1 (MQ=255)
                                                                                                                                            cGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                       >  2:139005/1‑118 (MQ=255)
                                                                                                                                            cGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa         <  2:943628/148‑1 (MQ=255)
                                                                                                                                            cGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                <  2:64540/141‑1 (MQ=255)
                                                                                                                                            cGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                >  1:64540/1‑141 (MQ=255)
                                                                                                                                              tCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt       <  1:795884/148‑1 (MQ=255)
                                                                                                                                              tCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa      >  3:33507/1‑149 (MQ=255)
                                                                                                                                               cACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa      >  2:213000/1‑148 (MQ=255)
                                                                                                                                                   cggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTGATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGGATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg  >  1:794413/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCG  >  NC_000913/3101393‑3101687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: