Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,172 G→A 6.0% intergenic (+31/+213) pitA → / ← uspB phosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,1720GA 6.0% 742.9 / 24.2 218intergenic (+31/+213)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base G (90/115);  new base A (4/9);  total (94/124)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 220 total aligned reads displayed.

GTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAG  >  NC_000913/3639026‑3639304
                                                                                                                                                  |                                                                                                                                    
gtggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGcc                                                                                                                                    <  1:1147282/149‑1 (MQ=255)
 tggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGccc                                                                                                                                   >  4:702699/1‑149 (MQ=255)
  ggcggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCg                                                                                                                                  <  1:676919/149‑1 (MQ=255)
   gcggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGc                                                                                                                                 <  4:321565/149‑1 (MQ=255)
   gcggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGc                                                                                                                                 <  4:1480746/149‑1 (MQ=255)
   gcggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGc                                                                                                                                 <  2:1663995/149‑1 (MQ=255)
     ggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCtt                                                                                                                               <  3:1347570/149‑1 (MQ=255)
     ggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCtt                                                                                                                               >  4:1428022/1‑149 (MQ=255)
     ggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGCtc                                                                                                                               <  2:1064661/149‑2 (MQ=255)
       cTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCa                                                                                                                             >  1:397977/1‑149 (MQ=255)
         tACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGa                                                                                                                           <  3:1587358/149‑1 (MQ=255)
            aGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtg                                                                                                                        >  2:1799013/1‑149 (MQ=255)
            aGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtg                                                                                                                        >  3:433370/1‑149 (MQ=255)
             gCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCGTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgt                                                                                                                       >  1:1695614/1‑149 (MQ=255)
             gCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGGTTCAGATtgt                                                                                                                       <  3:1167505/149‑1 (MQ=255)
             gCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgt                                                                                                                       <  3:395081/149‑1 (MQ=255)
             gCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgt                                                                                                                       <  1:1308531/149‑1 (MQ=255)
             gCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgt                                                                                                                       <  1:1799013/149‑1 (MQ=255)
              cGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgtg                                                                                                                      >  1:591007/1‑149 (MQ=255)
               gTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGa                                                                                                                     <  1:212329/149‑1 (MQ=255)
                  aaaCCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACAt                                                                                                                  >  1:667651/1‑149 (MQ=255)
                   aaCCGTGACCAGCGTTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATa                                                                                                                 <  3:130173/149‑1 (MQ=255)
                     ccGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGcc                                                                                                                                    >  3:628282/1‑128 (MQ=255)
                     ccGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGcc                                                                                                                                    <  4:628282/128‑1 (MQ=255)
                       gTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGCt                                                                                                                                <  4:1684822/130‑1 (MQ=255)
                            cAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGccc                                                                                                                                   >  1:1400554/1‑122 (MQ=255)
                            cAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGccc                                                                                                                                   <  2:1400554/122‑1 (MQ=255)
                            cAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCttt                                                                                                        >  1:722275/1‑149 (MQ=255)
                             aGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                                                                       <  3:1011412/149‑1 (MQ=255)
                              gCATTCTGATGGCCTGGGTGTTTACCCTTACGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGt                                                                                                      <  2:471780/149‑1 (MQ=255)
                               cATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGtt                                                                                                     <  3:702699/149‑1 (MQ=255)
                               cATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGtt                                                                                                     <  1:1292164/149‑1 (MQ=255)
                                aTTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTc                                                                                                    <  4:906994/149‑1 (MQ=255)
                                    tGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGc                                                                                                >  2:117851/1‑149 (MQ=255)
                                        ggCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGAtt                                                                                                                         <  3:1754117/120‑1 (MQ=255)
                                           cTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgtg                                                                                                                      >  2:920282/1‑120 (MQ=255)
                                           cTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTc                                                                                                    >  3:1474465/1‑138 (MQ=255)
                                           cTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgc                                                                                         >  3:1567967/1‑149 (MQ=255)
                                               gtgtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGt                                                                                                      >  2:670872/1‑132 (MQ=255)
                                               gtgtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggcg                                                                                     <  1:1112825/149‑1 (MQ=255)
                                                tgtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGt                                                                                                      >  1:1032103/1‑131 (MQ=255)
                                                 gtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCttt                                                                                                        <  3:1411586/128‑1 (MQ=255)
                                                   ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                 <  4:433370/149‑1 (MQ=255)
                                                    tACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                <  2:374752/149‑1 (MQ=255)
                                                      cccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCTCCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCACGCCGCATCCGGCATGaa                                                                                 >  2:195609/1‑146 (MQ=255)
                                                      cccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                              <  2:1419434/149‑1 (MQ=255)
                                                      cccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                              <  2:943376/149‑1 (MQ=255)
                                                       ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGa                                                                                                                           >  1:1375478/1‑103 (MQ=255)
                                                       ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGa                                                                                                                           <  2:1375478/103‑1 (MQ=255)
                                                       ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                             <  2:722275/149‑1 (MQ=255)
                                                        cTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggc                                                                                      <  4:1004089/139‑1 (MQ=255)
                                                        cTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                            <  1:117851/149‑1 (MQ=255)
                                                         ttCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCACCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                           >  1:983115/1‑149 (MQ=255)
                                                           ccGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                         <  1:1335184/149‑1 (MQ=255)
                                                                    gTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCGGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                <  1:461045/149‑1 (MQ=255)
                                                                    gTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                >  1:1578849/1‑149 (MQ=255)
                                                                     tACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                   <  2:545707/145‑1 (MQ=255)
                                                                     tACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                   >  1:545707/1‑145 (MQ=255)
                                                                        tGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGa                                                                                                                     <  4:107493/92‑1 (MQ=255)
                                                                        tGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGa                                                                                                                     >  3:107493/1‑92 (MQ=255)
                                                                            ttCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACAt                                                        <  2:1578849/149‑1 (MQ=255)
                                                                              ccGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGt                                                                                                               >  1:639762/1‑92 (MQ=255)
                                                                              ccGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGc                                                      <  4:177165/149‑1 (MQ=255)
                                                                               cggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTgc                                                                                                             <  1:1013112/93‑1 (MQ=255)
                                                                               cggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTgc                                                                                                             >  2:1013112/1‑93 (MQ=255)
                                                                               cggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                         <  1:905956/129‑1 (MQ=255)
                                                                                ggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgt                                                                                                                       <  4:297480/82‑1 (MQ=255)
                                                                                ggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgt                                                                                                                       >  3:297480/1‑82 (MQ=255)
                                                                                ggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCt                                                                                                          <  2:636689/95‑1 (MQ=255)
                                                                                ggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                    >  3:915061/1‑133 (MQ=255)
                                                                                ggcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                    <  2:648727/133‑1 (MQ=255)
                                                                                          aCTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggcg                                                                                     <  2:696271/106‑1 (MQ=255)
                                                                                           cTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAt                                         <  3:1722095/149‑1 (MQ=255)
                                                                                           cTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAt                                         <  1:34799/149‑1 (MQ=255)
                                                                                            tGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAtt                                        >  2:1409439/1‑149 (MQ=255)
                                                                                             ggCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTg                                       <  1:926601/149‑1 (MQ=255)
                                                                                               ctctCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCa                                     <  4:1567967/149‑1 (MQ=255)
                                                                                               ctctCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCa                                     <  2:4415/149‑1 (MQ=255)
                                                                                                 ctcCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCaga                                   <  1:253190/149‑1 (MQ=255)
                                                                                                     ttGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAtt                               >  4:1227705/1‑149 (MQ=255)
                                                                                                       gCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                    >  2:1349826/1‑110 (MQ=255)
                                                                                                          gTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                        >  4:1376112/1‑103 (MQ=255)
                                                                                                             ccTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGc                       <  2:983115/149‑1 (MQ=255)
                                                                                                               tGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                 <  3:1419345/89‑1 (MQ=255)
                                                                                                                 tAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGa                   <  2:591007/149‑1 (MQ=255)
                                                                                                                   aTCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCaa                                              <  3:1250965/120‑1 (MQ=255)
                                                                                                                   aTCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCa                                     <  4:1418660/129‑1 (MQ=255)
                                                                                                                              cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                >  2:1593514/1‑75 (MQ=255)
                                                                                                                              cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                <  1:1593514/75‑1 (MQ=255)
                                                                                                                              cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaaa                                                           <  4:1048927/96‑1 (MQ=255)
                                                                                                                              cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaaa                                                           >  3:1048927/1‑96 (MQ=255)
                                                                                                                              cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGc                       >  2:1146067/1‑132 (MQ=255)
                                                                                                                               aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                >  1:1439052/1‑74 (MQ=255)
                                                                                                                                aaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAgaga                                 <  3:900461/120‑1 (MQ=255)
                                                                                                                                aaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAgaga                                 >  4:900461/1‑120 (MQ=255)
                                                                                                                                   gaaGCGGGCCAGCTGACCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATa                                            >  2:787044/3‑106 (MQ=255)
                                                                                                                                    aaGCGGGCCAGCTGACCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACa                                                         <  1:627058/91‑1 (MQ=255)
                                                                                                                                         ggTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAg               <  4:1263354/129‑1 (MQ=255)
                                                                                                                                                  gCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                         >  2:767849/1‑110 (MQ=255)
                                                                                                                                                  gCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                         <  1:767849/110‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                    
GTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAG  >  NC_000913/3639026‑3639304

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: