Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,374,052 C→T 5.2% E146E (GAG→GAA nanR ← sialic acid‑inducible nan operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,0520CT 5.2% 254.4 / 3.6 78E146E (GAG→GAAnanRsialic acid‑inducible nan operon repressor
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base T (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-02

GCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAA  >  NC_000913/3373904‑3374196
                                                                                                                                                    |                                                                                                                                                
gCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGc                                                                                                                                                  <  3:192495/149‑1 (MQ=255)
   gggCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCa                                                                                                                                               >  1:173999/1‑149 (MQ=255)
      cAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCg                                                                                                                                            <  2:79645/149‑1 (MQ=255)
            ggATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCa                                                                                                                                      >  3:135820/1‑149 (MQ=255)
               tCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGt                                                                                                                                   >  2:120967/1‑149 (MQ=255)
                  ccATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCa                                                                                                                                >  3:358529/1‑149 (MQ=255)
                     tGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTg                                                                                                                                         >  2:40094/1‑137 (MQ=255)
                     tGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTg                                                                                                                                         <  1:40094/137‑1 (MQ=255)
                     tGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGcc                                                                                                                             <  2:173999/149‑1 (MQ=255)
                         gATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGTTCATCGGTGGCATGTTCAGCCGCAt                                                                                                                         >  1:214733/1‑149 (MQ=255)
                         gATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGTTCATCGGTGGCATGTTCAGCCGCAt                                                                                                                         >  3:44274/1‑149 (MQ=255)
                         gATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCAt                                                                                                                         >  3:114506/1‑149 (MQ=255)
                            tGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAgc                                                                                                                      <  1:120967/149‑1 (MQ=255)
                                       ggATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTg                                                                                                           >  3:122947/1‑149 (MQ=255)
                                         aTCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGa                                                                                                         >  3:238740/1‑149 (MQ=255)
                                         aTCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGa                                                                                                         <  4:290693/149‑1 (MQ=255)
                                             ccGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCa                                                                                                     >  4:222013/1‑149 (MQ=255)
                                               gCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGTTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCaaa                                                                                                   <  2:214733/149‑1 (MQ=255)
                                               gCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGTTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCaaa                                                                                                   <  4:44274/149‑1 (MQ=255)
                                               gCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCaaa                                                                                                   <  4:87767/149‑1 (MQ=255)
                                                   gCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAAc                                                                                               <  2:99842/149‑1 (MQ=255)
                                                           tGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTaa                                                                                       <  3:468209/149‑1 (MQ=255)
                                                             gAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAAtt                                                                                     <  3:218911/149‑1 (MQ=255)
                                                               aaTCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGt                                                                                   <  4:68267/149‑1 (MQ=255)
                                                                 tCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTc                                                                                 <  1:281892/149‑1 (MQ=255)
                                                                 tCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTc                                                                                 <  4:358529/149‑1 (MQ=255)
                                                                  cAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCg                                                                                >  1:170851/1‑149 (MQ=255)
                                                                       tCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATg                                                                           >  2:194595/1‑149 (MQ=255)
                                                                         tGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATggg                                                                         >  4:110261/1‑149 (MQ=255)
                                                                         tGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATgg                                                                          <  2:477818/148‑1 (MQ=255)
                                                                         tGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATgg                                                                          >  1:477818/1‑148 (MQ=255)
                                                                           aaCGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGaaa                                                                             >  4:314873/1‑143 (MQ=255)
                                                                           aaCGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGaaa                                                                             <  3:314873/143‑1 (MQ=255)
                                                                              gaatgaatGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGa                                                                                                         >  4:251412/1‑112 (MQ=255)
                                                                              gaatgaatGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGa                                                                                                         <  3:251412/112‑1 (MQ=255)
                                                                                atgaatgCCGCCTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCGAAGAACAGACGTAATTGTTCGAAATGGGCAATccc                                                                  <  3:308816/149‑1 (MQ=255)
                                                                                   aatgCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCAcc                                                               >  3:161503/1‑149 (MQ=255)
                                                                                   aatgCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCAcc                                                               >  4:369313/1‑149 (MQ=255)
                                                                                    atgCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCg                                                              >  2:463849/1‑149 (MQ=255)
                                                                                     tgCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCgg                                                             >  3:242505/1‑149 (MQ=255)
                                                                                       ccGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGa                                                           <  1:463849/149‑1 (MQ=255)
                                                                                        cgcgTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGAt                                                          <  2:170851/149‑1 (MQ=255)
                                                                                         gcgTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATg                                                         <  4:122947/149‑1 (MQ=255)
                                                                                          cgTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATga                                                        <  2:247661/149‑1 (MQ=255)
                                                                                            ttgttATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAATAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATgaga                                                      >  3:405908/1‑149 (MQ=255)
                                                                                            ttgttATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATgaga                                                      >  4:135900/1‑149 (MQ=255)
                                                                                                   ccAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCaa                                                                      >  2:182558/1‑126 (MQ=255)
                                                                                                   ccAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCaa                                                                      <  1:182558/126‑1 (MQ=255)
                                                                                                     aGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGaa                                                                                                >  4:360423/1‑98 (MQ=255)
                                                                                                     aGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGaa                                                                                                <  3:360423/98‑1 (MQ=255)
                                                                                                         actgactgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGaaa                                                    <  4:129276/138‑1 (MQ=255)
                                                                                                         actgactgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGaaa                                                    >  3:129276/1‑138 (MQ=255)
                                                                                                         actgactgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTc                                         <  3:135900/149‑1 (MQ=255)
                                                                                                           tgactgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGc                                       >  3:236714/1‑149 (MQ=255)
                                                                                                            gactgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGcc                                      >  1:79857/1‑149 (MQ=255)
                                                                                                               tgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGAt                                                          >  3:246841/1‑126 (MQ=255)
                                                                                                               tgTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGAt                                                          <  4:246841/126‑1 (MQ=255)
                                                                                                                 ttaaTTTCCAGTGCTTTTCCCAGCAAATCGACTTAATCATCGTTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCTCAGAACGGACGTAATTATTCGAAATGCGCAATCCCACCGGGATGAGAAAGGCAATCTTTAGCCAAGcc                                 >  4:200265/4‑149 (MQ=255)
                                                                                                                 ttGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACa                                                                                              <  2:390487/88‑1 (MQ=255)
                                                                                                                 ttGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACa                                                                                              >  1:390487/1‑88 (MQ=255)
                                                                                                                 ttGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGcc                                 >  2:213048/1‑149 (MQ=255)
                                                                                                                         cAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCt                         >  3:214198/1‑149 (MQ=255)
                                                                                                                                  tGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGAGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCgatga                <  3:343561/149‑1 (MQ=255)
                                                                                                                                    ccAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  3:139400/149‑1 (MQ=255)
                                                                                                                                    ccAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  4:161503/149‑1 (MQ=255)
                                                                                                                                    ccAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  2:79857/149‑1 (MQ=255)
                                                                                                                                    ccAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  4:238740/149‑1 (MQ=255)
                                                                                                                                    ccAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGAAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  3:46183/149‑1 (MQ=255)
                                                                                                                                      aGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAgt            >  4:63080/1‑149 (MQ=255)
                                                                                                                                      aGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAgt            >  2:92589/1‑149 (MQ=255)
                                                                                                                                        cAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAgtgt          <  2:6738/149‑1 (MQ=255)
                                                                                                                                         aaaTCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTc         <  4:170768/149‑1 (MQ=255)
                                                                                                                                           aTCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              >  4:40109/1‑142 (MQ=255)
                                                                                                                                           aTCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  3:40109/142‑1 (MQ=255)
                                                                                                                                           aTCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              <  1:90601/142‑1 (MQ=255)
                                                                                                                                           aTCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATa              >  2:90601/1‑142 (MQ=255)
                                                                                                                                            tCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGc      >  3:251518/1‑149 (MQ=255)
                                                                                                                                                tttGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGaa  <  4:214198/149‑1 (MQ=255)
                                                                                                                                                tttGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGaa  <  4:135820/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCAGGGCAACGTGGATCGCCATGAAGATTGGGTTACCGGGGATCTCCGCCAGCACGCGGTGGAAATCAACGTCTGAACGAATGAATGCCGCGTTGTTATCCAGCGACTGACTGTTGATTTCCAGTGCTTTTGCCAGCAAATCGATTTGCTCATCGGTGGCATGTTCAGCCGCATAGCGCACCAGACTGGATTCAAAGAACAGACGTAATTGTTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAA  >  NC_000913/3373904‑3374196

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: