Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,191 C→G 6.9% intergenic (+25/+28) yicH → / ← yicI putative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,1910CG 6.9% 2203.2 / 110.0 607intergenic (+25/+28)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (299/265);  new base G (24/18);  total (324/283)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 615 total aligned reads displayed.

GTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAG  >  NC_000913/3832043‑3832339
                                                                                                                                                    |                                                                                                                                                    
 ttGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                                                                     <  3:1341537/149‑1 (MQ=255)
   gAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                   >  2:1459026/1‑149 (MQ=255)
   gAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                   >  2:33765/1‑149 (MQ=255)
        cTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                                                                              >  4:1842970/1‑149 (MQ=255)
         tACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcc                                                                                                                                             >  4:1896639/1‑149 (MQ=255)
          aCGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                                                                            <  2:1623183/149‑1 (MQ=255)
           cGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                                                                           <  4:1639864/149‑1 (MQ=255)
            ggTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                                                                          <  3:538997/149‑1 (MQ=255)
              ttAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                                                                   <  3:629270/138‑1 (MQ=255)
              ttAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGt                                                                                                                                                <  4:137383/141‑1 (MQ=255)
                  cGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                                                                               >  4:268493/1‑138 (MQ=255)
                        aTTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                                                                              >  3:738995/1‑149 (MQ=255)
                         ttGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                                                                  <  2:263692/128‑1 (MQ=255)
                         ttGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                             <  4:1607532/149‑1 (MQ=255)
                         ttGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                             >  4:721779/1‑149 (MQ=255)
                             aTGCCGTGAACGCCGCAAAGCAGCTATTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                                                                         >  2:1175874/1‑149 (MQ=255)
                             aTGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAccc                                                                                                                                            >  2:102800/1‑130 (MQ=255)
                                     aaCGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATAGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                                                                                 <  3:1617682/149‑1 (MQ=255)
                                          cgcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGt                                                                                                            >  1:8610/1‑149 (MQ=255)
                                           gcAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                                                                        <  2:1550913/120‑1 (MQ=255)
                                            cAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcc                                                                                                                                             >  3:955614/1‑114 (MQ=255)
                                             aaaGCAGCAAGTTAGTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAg                                                                                                         >  3:137731/1‑149 (MQ=255)
                                                  agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                                                                    <  3:1773064/149‑1 (MQ=255)
                                                         gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                                                             >  1:77997/1‑149 (MQ=255)
                                                          aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                            <  1:43941/149‑1 (MQ=255)
                                                           cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                                                                           >  1:1717886/1‑149 (MQ=255)
                                                            tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTATGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                                                                          >  2:1852543/1‑149 (MQ=255)
                                                            tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                                                                          >  1:139869/1‑149 (MQ=255)
                                                            tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                                                                          >  2:137918/1‑149 (MQ=255)
                                                              aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTACCTTg                                                                                             >  1:1003195/1‑144 (MQ=255)
                                                                aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTc                                                                                      >  4:105531/1‑149 (MQ=255)
                                                                  cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                                                                    >  1:1573766/1‑149 (MQ=255)
                                                                   aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCa                                                                                                          <  2:1314451/126‑1 (MQ=255)
                                                                   aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                            >  1:183744/1‑140 (MQ=255)
                                                                     gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcaccac                                                                                 >  3:1549313/1‑149 (MQ=255)
                                                                         gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGc                                                                             <  2:766603/149‑1 (MQ=255)
                                                                           tGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                                                                           <  3:332740/149‑1 (MQ=255)
                                                                           tGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                                                                           <  3:1173770/149‑1 (MQ=255)
                                                                            ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                                                                           <  1:295492/132‑1 (MQ=255)
                                                                             gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                            <  3:1273535/130‑1 (MQ=255)
                                                                               aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                            >  2:99977/1‑128 (MQ=255)
                                                                               aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                                                                       >  1:1745865/1‑149 (MQ=255)
                                                                                tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCa                                                                                                          >  3:1323027/1‑113 (MQ=255)
                                                                                tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCa                                                                                                          <  4:1323027/113‑1 (MQ=255)
                                                                                tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCa                                                                                                          >  3:799462/1‑113 (MQ=255)
                                                                                    ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                                                                  <  4:1806780/149‑1 (MQ=255)
                                                                                     tCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGATAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATCGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAcat                                                                 <  2:663495/149‑3 (MQ=255)
                                                                                        aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa                                                              <  2:432760/149‑1 (MQ=255)
                                                                                           ggTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                                                                            <  4:587950/116‑1 (MQ=255)
                                                                                            gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                                                                           >  4:127056/1‑116 (MQ=255)
                                                                                            gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCt                                                          >  3:1075059/1‑149 (MQ=255)
                                                                                             tGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg                                                         <  2:859075/149‑1 (MQ=255)
                                                                                                  ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                             >  4:746333/1‑76 (MQ=255)
                                                                                                  ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                                                                    <  1:1241140/85‑1 (MQ=255)
                                                                                                  ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc                                                    >  4:890854/1‑149 (MQ=255)
                                                                                                      aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAc                                                                                                                                              >  4:348028/1‑55 (MQ=37)
                                                                                                      aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGt                                                 >  2:468556/1‑148 (MQ=255)
                                                                                                      aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTc                                                >  3:863989/1‑149 (MQ=255)
                                                                                                        cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCtt                                              >  3:1211578/1‑149 (MQ=255)
                                                                                                          aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                                                                             >  1:1369165/1‑68 (MQ=255)
                                                                                                           cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGc                                                                             <  1:1473423/115‑1 (MQ=255)
                                                                                                           cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCa                                           <  4:499776/149‑1 (MQ=255)
                                                                                                           cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCa                                           <  3:417129/149‑1 (MQ=255)
                                                                                                            ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc                                                            <  3:1449816/131‑1 (MQ=255)
                                                                                                               gCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGAcc                                       <  3:1691228/149‑1 (MQ=255)
                                                                                                                    tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCTGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTa                                  <  2:210585/149‑1 (MQ=255)
                                                                                                                    tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTa                                  <  2:1745865/149‑1 (MQ=255)
                                                                                                                        gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCa                                                                                                          >  4:44492/1‑73 (MQ=255)
                                                                                                                        gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                                                                           <  3:474502/88‑1 (MQ=255)
                                                                                                                        gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACttt                              >  2:1793416/1‑149 (MQ=255)
                                                                                                                         tAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                                                                              >  1:602075/1‑84 (MQ=255)
                                                                                                                         tAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                                                                              <  2:575217/84‑1 (MQ=255)
                                                                                                                         tAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCtt                                              >  2:296312/1‑132 (MQ=255)
                                                                                                                            aCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTcaca                          >  1:1057166/1‑149 (MQ=255)
                                                                                                                               tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                                                                          >  1:222848/1‑98 (MQ=255)
                                                                                                                                 aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                                                                                 >  1:1789465/1‑57 (MQ=255)
                                                                                                                                 aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCa                                           <  3:709231/127‑1 (MQ=255)
                                                                                                                                 aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCa                                           >  4:1482965/1‑127 (MQ=255)
                                                                                                                                 aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACAtt                     >  1:701167/1‑149 (MQ=255)
                                                                                                                                  aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAtt                                                                                                              >  1:1126833/1‑59 (MQ=255)
                                                                                                                                   tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                                                                  >  1:1631304/1‑102 (MQ=255)
                                                                                                                                       gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACttt                              <  4:1206407/134‑1 (MQ=255)
                                                                                                                                        ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                                                                <  2:917648/99‑1 (MQ=255)
                                                                                                                                        ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                                                                >  1:917648/1‑99 (MQ=255)
                                                                                                                                         cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTc                                                <  2:1706242/114‑1 (MQ=255)
                                                                                                                                         cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTgc                   <  2:1070862/143‑1 (MQ=255)
                                                                                                                                         cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGc             >  1:325709/1‑149 (MQ=255)
                                                                                                                                         cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTg                    >  2:1220198/1‑142 (MQ=255)
                                                                                                                                         cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGc             >  4:695159/1‑149 (MQ=255)
                                                                                                                                          gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg                                                             <  3:1889992/100‑1 (MQ=255)
                                                                                                                                            gTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCacaa                         <  3:1839920/134‑1 (MQ=255)
                                                                                                                                            gTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGcaca          >  4:1159186/1‑149 (MQ=255)
                                                                                                                                            gTACGTCTGAGGCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCacaa                         >  2:958701/1‑134 (MQ=255)
                                                                                                                                            gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGcaca          <  3:234602/149‑1 (MQ=255)
                                                                                                                                                gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAgtgt      >  4:1704504/1‑149 (MQ=255)
                                                                                                                                                  cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTcc    <  1:47073/149‑1 (MQ=255)
                                                                                                                                                   tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCa   >  1:1879578/1‑149 (MQ=255)
                                                                                                                                                   tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCa   >  1:366684/1‑149 (MQ=255)
                                                                                                                                                   tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCa   >  1:966517/1‑149 (MQ=255)
                                                                                                                                                   tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCa   >  4:829600/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAG  >  NC_000913/3832043‑3832339

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: