Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,313,739 A→C E130D (GAA→GAC rpiB → ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,7390AC98.5% 234.2 / ‑4.3 66E130D (GAA→GACrpiBribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (24/41);  total (24/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTT  >  NC_000913/4313591‑4313878
                                                                                                                                                    |                                                                                                                                           
gAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACgg                                                                                                                                           >  1:686024/1‑151 (MQ=255)
 aaTTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGc                                                                                                                                          <  2:741166/151‑1 (MQ=255)
  aTTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCg                                                                                                                                         >  2:910368/1‑151 (MQ=255)
     cgcgcgGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGgtc                                                                                                                                      >  1:173830/1‑151 (MQ=255)
     cgcgcgGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGgtc                                                                                                                                      >  1:409741/1‑151 (MQ=255)
      gcgcgGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGgtcg                                                                                                                                     <  1:175877/151‑1 (MQ=255)
      gcgcgGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGgtcg                                                                                                                                     <  2:363890/151‑1 (MQ=255)
       cgcgGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGgtcgt                                                                                                                                    >  1:42515/1‑151 (MQ=255)
              gtcTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATcaac                                                                                                                             >  2:508255/1‑151 (MQ=255)
                 tGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCaacaac                                                                                                                          >  2:421252/1‑151 (MQ=255)
                  gTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATcaaca                                                                                                                            <  2:501770/148‑1 (MQ=255)
                      cGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTg                                                                                                                     >  1:805680/1‑151 (MQ=255)
                        aaCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTgg                                                                                                                    >  1:662144/1‑150 (MQ=255)
                         aCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAg                                                                                                                  <  1:899236/151‑1 (MQ=255)
                         aCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAg                                                                                                                  <  1:593428/151‑1 (MQ=255)
                              aTTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGAt                                                                                                             >  1:306857/1‑151 (MQ=255)
                                 ccGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTa                                                                                                           <  1:677902/150‑1 (MQ=255)
                                  cgcgCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACg                                                                                                         >  2:861371/1‑151 (MQ=255)
                                  cgcgCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACg                                                                                                         <  1:362362/151‑1 (MQ=255)
                                   gcgcAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACgg                                                                                                        <  2:577242/151‑1 (MQ=255)
                                        aCTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATa                                                                                                   <  2:554960/151‑1 (MQ=255)
                                        aCTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATa                                                                                                   <  2:377333/151‑1 (MQ=255)
                                           ttCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATag                                                                                                  <  2:176121/149‑1 (MQ=255)
                                             cgcgGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGagca                                                                                              <  1:706869/151‑1 (MQ=255)
                                              gcgGCAGCACAGCGACACCAACGTGCTGGCTTTTGGTTCACGAGGGTTTGCGCTCGAACGGGCAAAAATGATTGTGGATTCGTGGCTTGGCGCACAGTGCGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATGgaggag                                                                                             <  2:456926/151‑1 (MQ=255)
                                                 gcagcaTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCgg                                                                                          <  1:280958/151‑1 (MQ=255)
                                                   agcaTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGgag                                                                                        <  2:239727/151‑1 (MQ=255)
                                                        aaCGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAAt                                                                                    <  2:910616/150‑1 (MQ=255)
                                                         aCGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTg                                                                                  <  1:233934/151‑1 (MQ=255)
                                                           gACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTgag                                                                                <  2:455090/151‑1 (MQ=255)
                                                                    gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATcc                                                                        >  2:704221/1‑150 (MQ=255)
                                                                    gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCa                                                                       <  1:496673/151‑1 (MQ=255)
                                                                         ggCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCacta                                                                    <  1:194425/149‑1 (MQ=255)
                                                                          gCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACtt                                                                 <  1:547309/151‑1 (MQ=255)
                                                                             tttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGc                                                               <  1:786469/150‑1 (MQ=255)
                                                                              ttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCa                                                              <  2:188755/150‑1 (MQ=255)
                                                                                ggTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATg                                                            <  2:449153/150‑1 (MQ=255)
                                                                                ggTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATg                                                            <  2:827476/150‑1 (MQ=255)
                                                                                   tCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATg                                                        >  2:159350/1‑151 (MQ=255)
                                                                                      cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATg                                                        >  2:859533/1‑148 (MQ=255)
                                                                                       gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAg                                                      >  2:264864/1‑149 (MQ=255)
                                                                                       gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGt                                                     >  2:473764/1‑150 (MQ=255)
                                                                                       gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTa                                                    <  2:158509/151‑1 (MQ=255)
                                                                                         gTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAAt                                                  <  1:387325/151‑1 (MQ=255)
                                                                                           ggTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATg                                                 <  1:536206/150‑1 (MQ=255)
                                                                                           ggTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGa                                                <  2:528247/151‑1 (MQ=255)
                                                                                              tGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGAtt                                              >  2:635588/1‑150 (MQ=255)
                                                                                               ggCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTa                                             <  1:485712/150‑1 (MQ=255)
                                                                                               ggCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCCGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTa                                             <  1:196267/150‑1 (MQ=255)
                                                                                                       aCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATa                                    <  1:87930/151‑1 (MQ=255)
                                                                                                        cTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGAt                                     >  2:814378/1‑149 (MQ=255)
                                                                                                          ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATag                                   <  2:902612/149‑1 (MQ=255)
                                                                                                           gCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGt                                <  2:803940/151‑1 (MQ=255)
                                                                                                            cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGtt                               >  1:325747/1‑151 (MQ=255)
                                                                                                             aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGttt                              >  2:824363/1‑151 (MQ=255)
                                                                                                              aaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTc                             >  1:749911/1‑151 (MQ=255)
                                                                                                              aaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTc                             <  2:379313/151‑1 (MQ=255)
                                                                                                               aaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCt                            <  2:923563/151‑1 (MQ=255)
                                                                                                                      tgtgGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTg                       >  1:921754/1‑149 (MQ=255)
                                                                                                                         ggATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                  >  1:402393/1‑151 (MQ=255)
                                                                                                                         ggATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                  <  1:367770/151‑1 (MQ=255)
                                                                                                                           aTGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                <  1:271477/151‑1 (MQ=255)
                                                                                                                           aTGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                <  2:28129/151‑1 (MQ=255)
                                                                                                                                 ggCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATg          >  1:455600/1‑151 (MQ=255)
                                                                                                                                         gcACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCtt  <  1:605419/151‑1 (MQ=255)
                                                                                                                                         gcACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCtt  <  2:291738/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTT  >  NC_000913/4313591‑4313878

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