Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,150 C→T 7.1% intergenic (+162/+28) sibC → / ← serA sRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1500CT 7.1% 267.9 / 8.8 85intergenic (+162/+28)sibC/serAsRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/48);  new base T (4/2);  total (35/50)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

CATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAG  >  NC_000913/3057002‑3057296
                                                                                                                                                    |                                                                                                                                                  
cATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTTTTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                                                                                    >  3:902410/1‑149 (MQ=255)
    ctgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                                                                                 <  2:568333/148‑1 (MQ=255)
     tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGcc                                                                                                                                                <  2:481483/148‑1 (MQ=255)
        cctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGt                                                                                                                                             >  1:219363/1‑148 (MQ=255)
        cctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGt                                                                                                                                             <  1:328313/148‑1 (MQ=255)
          tgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTg                                                                                                                                           >  2:734406/1‑148 (MQ=255)
              aTCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAttt                                                                                                                                       <  2:603713/148‑1 (MQ=255)
               tCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                                                                                      <  1:912695/148‑1 (MQ=255)
               tCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                                                                                      >  1:507139/1‑148 (MQ=255)
                cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                     <  2:108999/148‑1 (MQ=255)
                cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                                                                    <  4:842293/149‑1 (MQ=255)
                cGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGAAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                     >  1:641452/1‑148 (MQ=255)
                 gTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                                                                   >  3:398132/1‑149 (MQ=255)
                  tAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                                                                   <  1:396792/148‑1 (MQ=255)
                     aTCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGaa                                                                                                                               <  4:494295/149‑1 (MQ=255)
                     aTCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGaa                                                                                                                               <  3:111814/149‑1 (MQ=255)
                        cccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAggg                                                                                                                            <  3:362507/149‑1 (MQ=255)
                           aTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                     <  3:693925/137‑1 (MQ=255)
                           aTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                     >  4:693925/1‑137 (MQ=255)
                           aTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGa                                                                                                                          >  2:247549/1‑148 (MQ=255)
                           aTTAAATAAACACAACGCATTGATCTGACGTTGATTTATTTTCTGGAGCAGGCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                         >  3:404442/1‑149 (MQ=255)
                            ttAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                         <  2:587895/148‑1 (MQ=255)
                                 taaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTg                                                                                                                                                  <  2:767574/118‑1 (MQ=255)
                                 taaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTg                                                                                                                                                  >  1:767574/1‑118 (MQ=255)
                                   aaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                                    <  2:217384/146‑1 (MQ=255)
                                   aaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                                                                    >  1:217384/1‑146 (MQ=255)
                                       caACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTg                                                                                                                                           >  1:52109/1‑119 (MQ=255)
                                       caACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTg                                                                                                                                           <  2:52109/119‑1 (MQ=255)
                                       caACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagca                                                                                                              >  1:552694/1‑148 (MQ=255)
                                       caACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGGGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagca                                                                                                              >  2:732383/1‑148 (MQ=255)
                                         aCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGa                                                                                                            >  1:117848/1‑148 (MQ=255)
                                          cGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACg                                                                                                          <  4:398132/149‑1 (MQ=255)
                                             aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACggg                                                                                                        <  1:449791/148‑1 (MQ=255)
                                               tGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGCGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                                                                      >  2:447184/1‑148 (MQ=255)
                                                 aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                                                                                    >  1:240273/1‑148 (MQ=255)
                                                  tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                                                                  <  3:592165/149‑1 (MQ=255)
                                                  tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGa                                                                                                   >  1:884448/1‑148 (MQ=255)
                                                   cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                                                                  <  1:460463/148‑1 (MQ=255)
                                                   cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                                                                  >  1:270054/1‑148 (MQ=255)
                                                       cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGt                                                                                              >  1:125085/1‑148 (MQ=255)
                                                        tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                                                                             >  1:671440/1‑148 (MQ=255)
                                                         ttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                                                                            >  1:390325/1‑148 (MQ=255)
                                                          tGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAccc                                                                                          <  4:670834/149‑1 (MQ=255)
                                                            atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCg                                                                                         <  2:125085/148‑1 (MQ=255)
                                                             tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGa                                                                                       >  3:752812/1‑149 (MQ=255)
                                                              ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGaa                                                                                      <  3:45830/149‑1 (MQ=255)
                                                                atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATa                                                                                    <  3:819322/149‑1 (MQ=255)
                                                                         aGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAt                                                                            <  2:357638/148‑1 (MQ=255)
                                                                           cAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTg                                                                          <  2:671440/148‑1 (MQ=255)
                                                                            aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                                         <  1:247549/148‑1 (MQ=255)
                                                                            aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                                         <  2:853955/148‑1 (MQ=255)
                                                                                tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                    <  3:542548/149‑1 (MQ=255)
                                                                                 cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                             >  2:464978/1‑91 (MQ=255)
                                                                                 cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                             <  1:464978/91‑1 (MQ=255)
                                                                                 cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAg                                                                   <  4:552897/149‑1 (MQ=255)
                                                                                    aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc                                                                <  3:266392/149‑1 (MQ=255)
                                                                                          gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                                           <  1:400895/148‑1 (MQ=255)
                                                                                          gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTc                                                          >  3:296928/1‑149 (MQ=255)
                                                                                          gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTc                                                          <  4:503079/149‑1 (MQ=255)
                                                                                             tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGc                                                       <  3:590179/149‑1 (MQ=255)
                                                                                                 cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACg                                                   <  4:112867/149‑1 (MQ=255)
                                                                                                  aaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACg                                                   <  1:578561/148‑1 (MQ=255)
                                                                                                  aaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt                                                  <  4:296928/149‑1 (MQ=255)
                                                                                                   aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt                                                  >  1:84922/1‑148 (MQ=255)
                                                                                                   aCGCGGCAACGGTGTGGAGAAGGGACAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTc                                                 >  3:644538/1‑149 (MQ=255)
                                                                                                         cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgt                                            <  1:734406/148‑1 (MQ=255)
                                                                                                         cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgt                                            <  2:884448/148‑1 (MQ=255)
                                                                                                          aaCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtcg                                          <  4:269022/149‑1 (MQ=255)
                                                                                                            cGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCgg                                         <  2:240273/148‑1 (MQ=255)
                                                                                                             ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGc                                        <  2:390325/148‑1 (MQ=255)
                                                                                                                   gagaAGGGATAAAAAAACGGGTAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACTGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaaa                                  <  2:641452/148‑2 (MQ=255)
                                                                                                                   gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaata                                 <  3:332765/149‑1 (MQ=255)
                                                                                                                         ggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                    >  3:27194/1‑108 (MQ=255)
                                                                                                                         ggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                    <  4:27194/108‑1 (MQ=255)
                                                                                                                          gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATaa                          >  4:598989/1‑149 (MQ=255)
                                                                                                                           aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATaa                          >  1:226801/1‑148 (MQ=255)
                                                                                                                           aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATaa                          >  1:648735/1‑148 (MQ=255)
                                                                                                                                  aaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                  >  3:281559/1‑149 (MQ=255)
                                                                                                                                  aaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                  <  4:404442/149‑1 (MQ=255)
                                                                                                                                  aaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATa                   <  2:507139/148‑1 (MQ=255)
                                                                                                                                   aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAAc                 >  4:107213/1‑149 (MQ=255)
                                                                                                                                      ggCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCa              >  4:278019/1‑149 (MQ=255)
                                                                                                                                           gTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTg          >  4:4763/1‑148 (MQ=255)
                                                                                                                                           gTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTg          <  3:4763/148‑1 (MQ=255)
                                                                                                                                            tCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTgg         <  2:84922/148‑1 (MQ=255)
                                                                                                                                              aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCg      >  4:701589/1‑149 (MQ=255)
                                                                                                                                                   ccTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAg  >  2:114166/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAG  >  NC_000913/3057002‑3057296

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