Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,165 T→C 6.3% intergenic (+24/+220) 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,1650TC 6.3% 322.5 / 8.3 95intergenic (+24/+220)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (44/45);  new base C (3/3);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

AGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCA  >  NC_000913/3639021‑3639259
                                                                                                                                                |                                                                                              
aGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAg                                                                                             <  2:36798/148‑1 (MQ=255)
    ggtggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCt                                                                                           >  1:321886/1‑146 (MQ=255)
    ggtggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCt                                                                                           <  2:321886/146‑1 (MQ=255)
                                  aGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCttt                                                           >  2:703598/1‑148 (MQ=255)
                                   gCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgtg                                                                         <  1:358141/133‑1 (MQ=255)
                                   gCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgtg                                                                         >  2:358141/1‑133 (MQ=255)
                                          gATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCacgc                                                   >  2:273696/1‑145 (MQ=255)
                                            tGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACAt                                                                     >  3:780883/1‑128 (MQ=255)
                                            tGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACAt                                                                     <  4:780883/128‑1 (MQ=255)
                                                 tGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCt                                                                                   >  3:167586/1‑109 (MQ=255)
                                                 tGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCt                                                                                   <  4:167586/109‑1 (MQ=255)
                                                 tGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTg                                                                 <  3:1051772/127‑1 (MQ=255)
                                                 tGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTg                                                                 >  4:1051772/1‑127 (MQ=255)
                                                 tGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgc                                            >  4:624748/1‑148 (MQ=255)
                                                   ggTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCa                                                                                              <  2:993620/96‑1 (MQ=255)
                                                   ggTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCa                                                                                              >  1:993620/1‑96 (MQ=255)
                                                      gtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGt                                       >  1:435967/1‑148 (MQ=255)
                                                       tttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCACGCCGCATCCGGCATg                                      <  1:273696/148‑1 (MQ=255)
                                                        ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          >  2:126087/1‑127 (MQ=255)
                                                        ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          <  1:126087/127‑1 (MQ=255)
                                                        ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      >  1:643520/1‑131 (MQ=255)
                                                        ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      <  2:643520/131‑1 (MQ=255)
                                                        ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGa                                     >  3:863512/1‑148 (MQ=255)
                                                           cccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          >  4:535441/1‑124 (MQ=255)
                                                           cccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          <  3:535441/124‑1 (MQ=255)
                                                            ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCt                                                             >  2:821743/1‑120 (MQ=255)
                                                            ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCt                                                             <  1:821743/120‑1 (MQ=255)
                                                            ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCt                                                             <  3:267706/120‑1 (MQ=255)
                                                            ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCt                                                             >  4:267706/1‑120 (MQ=255)
                                                            ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                 >  3:1133567/1‑148 (MQ=255)
                                                              ttCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                               <  1:703598/148‑1 (MQ=255)
                                                              ttCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                               <  3:624748/148‑1 (MQ=255)
                                                               tCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGt                                       >  1:364550/1‑139 (MQ=255)
                                                               tCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGt                                       <  2:364550/139‑1 (MQ=255)
                                                               tCCGGCTGCGGTACTGCTTGCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                              <  2:435967/148‑1 (MQ=255)
                                                                     tGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                        <  4:863512/148‑1 (MQ=255)
                                                                       cGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAACACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          <  2:119233/112‑1 (MQ=255)
                                                                       cGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          >  1:119233/1‑112 (MQ=255)
                                                                           aCTGCTTTCCGGTGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGc                                                                                    >  1:91171/1‑82 (MQ=255)
                                                                           aCTGCTTTCCGGTGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGc                                                                                    <  2:91171/82‑1 (MQ=255)
                                                                           aCTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtg                                                                           >  1:669463/1‑91 (MQ=255)
                                                                           aCTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtg                                                                           <  2:669463/91‑1 (MQ=255)
                                                                             tGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACAt                                                                     >  1:822146/1‑95 (MQ=255)
                                                                             tGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACAt                                                                     <  2:822146/95‑1 (MQ=255)
                                                                               cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTc                                                                                               >  4:1035949/1‑67 (MQ=255)
                                                                               cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTc                                                                                               <  3:1035949/67‑1 (MQ=255)
                                                                               cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggcg                                        <  3:448733/122‑1 (MQ=255)
                                                                               cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggcg                                        >  4:448733/1‑122 (MQ=255)
                                                                               cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggcg                                        <  4:982363/122‑1 (MQ=255)
                                                                               cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTAGTTCATGCCGGATgcggcg                                        >  3:982363/1‑122 (MQ=255)
                                                                                tttCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCgg                                                >  1:402130/1‑113 (MQ=255)
                                                                                tttCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCgg                                                <  2:402130/113‑1 (MQ=255)
                                                                                      gcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                 <  2:763389/122‑1 (MQ=255)
                                                                                      gcggGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                 >  1:763389/1‑122 (MQ=255)
                                                                                           cTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCa  <  4:1133567/148‑1 (MQ=255)
                                                                                             gTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      >  4:633973/1‑94 (MQ=255)
                                                                                             gTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      <  3:633973/94‑1 (MQ=255)
                                                                                                 tGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      <  1:1062823/90‑1 (MQ=255)
                                                                                                 tGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      >  2:1062823/1‑90 (MQ=255)
                                                                                                 tGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCACGCCGCATCCGGCATGa                                     >  4:805621/1‑107 (MQ=255)
                                                                                                 tGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCACGCCGCATCCGGCATGa                                     <  3:805621/107‑1 (MQ=255)
                                                                                                    ctctCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                    >  3:556695/1‑105 (MQ=255)
                                                                                                    ctctCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                    <  4:556695/105‑1 (MQ=255)
                                                                                                          ttGCAGTTCCTGTAATCGTACGCACCAAAACAAGCGGGCCAGCTGACCCGCtcg                                                                                 >  3:251415/1‑52 (MQ=25)
                                                                                                          ttGCAGTTCCTGTAATCGTACGCACCAAAACAAGCGGGCCAGCTGACCCGCtcg                                                                                 <  4:251415/54‑3 (MQ=21)
                                                                                                           tGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCACGCCGCATccggcat                                       <  4:531574/95‑3 (MQ=255)
                                                                                                           tGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCACGCCGCATccggcat                                       >  3:531574/1‑93 (MQ=255)
                                                                                                             cAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                               <  4:700711/101‑1 (MQ=255)
                                                                                                             cAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                               >  3:700711/1‑101 (MQ=255)
                                                                                                              aGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTc                                                                                 <  2:559266/50‑1 (MQ=255)
                                                                                                              aGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTc                                                                                 >  1:559266/1‑50 (MQ=255)
                                                                                                                 tCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          >  3:173054/1‑70 (MQ=255)
                                                                                                                 tCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTg                                                          <  4:173054/70‑1 (MQ=255)
                                                                                                                 tCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCgg                                                >  3:1065485/1‑80 (MQ=255)
                                                                                                                 tCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCgg                                                <  4:1065485/80‑1 (MQ=255)
                                                                                                                     gTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggc                                         >  2:586319/1‑83 (MQ=255)
                                                                                                                     gTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggc                                         <  1:586319/83‑1 (MQ=255)
                                                                                                                              cGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                            >  2:551546/1‑87 (MQ=255)
                                                                                                                              cGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                            <  1:551546/87‑1 (MQ=255)
                                                                                                                                   cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaaa              >  1:420112/1‑96 (MQ=255)
                                                                                                                                   cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaaa              <  2:420112/96‑1 (MQ=255)
                                                                                                                                    aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAAt     <  1:457263/104‑1 (MQ=255)
                                                                                                                                    aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAAt     >  2:1100504/1‑104 (MQ=255)
                                                                                                                                    aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAAt     >  2:457263/1‑104 (MQ=255)
                                                                                                                                    aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAAt     <  1:1100504/104‑1 (MQ=255)
                                                                                                                                     aaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                       >  4:215306/1‑85 (MQ=255)
                                                                                                                                     aaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                       <  3:215306/85‑1 (MQ=255)
                                                                                                                                      aaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCa  >  3:916554/1‑105 (MQ=255)
                                                                                                                                      aaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCa  <  4:916554/105‑1 (MQ=255)
                                                                                                                                       aCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa               >  3:751876/1‑91 (MQ=255)
                                                                                                                                       aCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa               <  4:751876/91‑1 (MQ=255)
                                                                                                                                       aCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAAt     <  1:676471/101‑1 (MQ=255)
                                                                                                                                       aCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAAt     >  2:676471/1‑101 (MQ=255)
                                                                                                                                              ggTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCCTACCGGATgcggcg                                        >  3:655532/1‑59 (MQ=38)
                                                                                                                                              ggTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATACCGGATgcggcg                                        <  4:655532/59‑1 (MQ=255)
                                                                                                                                                tCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      >  4:646268/1‑43 (MQ=255)
                                                                                                                                                tCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                      <  3:646268/43‑1 (MQ=255)
                                                                                                                                                |                                                                                              
AGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCA  >  NC_000913/3639021‑3639259

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: