Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,868,477 A→G L294P (CTG→CCG)  b3691 ← D‑galactonate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,868,4770AG98.9% 591.1 / ‑3.5 185L294P (CTG→CCG) b3691D‑galactonate transport
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (102/81);  total (104/81)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 196 total aligned reads displayed.

TCAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCAGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGTT  >  NC_000913/3868388‑3868565
                                                                                         |                                                                                        
tCAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCgggc                                                                                       >  1:4012443/1‑92 (MQ=255)
tCAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGa                                                                                       <  1:322551/93‑1 (MQ=255)
tCAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGa                                                                                       <  2:627144/93‑1 (MQ=255)
 cAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGAc                                                                                      >  1:3863857/1‑93 (MQ=255)
 cAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGAc                                                                                      <  2:7323/93‑1 (MQ=255)
     aaGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCg                                                                                  >  2:1693983/1‑93 (MQ=255)
      aGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGa                                                                                 <  2:774686/93‑1 (MQ=255)
       gtcGCAGATAATAGGCGTTTTACGCGCAAAGCCCAGTGCAAAGCCCTTAAGTACCAGCAGATCCACTACCTAGCCAGAGAGCGGgacgacgac                                                                                >  2:3298368/3‑93 (MQ=255)
        ccGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACa                                                                               <  1:4685614/93‑1 (MQ=255)
         cGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACaa                                                                              <  1:1008313/93‑1 (MQ=255)
         cGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACaa                                                                              >  1:4059705/1‑93 (MQ=255)
          gCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACaaa                                                                             >  2:4177948/1‑93 (MQ=255)
           cAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTGCCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAc                                                                            <  2:2295376/93‑1 (MQ=255)
             gATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgc                                                                          >  1:4354790/1‑93 (MQ=255)
             gATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgc                                                                          >  2:2261552/1‑93 (MQ=255)
               tAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgccg                                                                        >  2:4682523/1‑93 (MQ=255)
                aaTCGGCGTTTTACGCGCAAAGCCCAGTGCAAAACCCTTACGTACCAGCACGTCCGCTACACAACCAAAGAGCGGGACGCCGACAGAcgccgc                                                                       >  2:3022150/1‑93 (MQ=255)
                aaTCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgccgc                                                                       >  2:3494823/1‑93 (MQ=255)
                aaTCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgccgc                                                                       >  1:4280622/1‑93 (MQ=255)
                 aTCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCcgcg                                                                      >  1:4502498/1‑93 (MQ=255)
                    ggCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAgg                                                                   >  2:553270/1‑93 (MQ=255)
                     gCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGa                                                                  <  2:4439766/93‑1 (MQ=255)
                      cGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGaa                                                                 >  2:4475657/1‑93 (MQ=255)
                      cGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGaa                                                                 <  2:4111008/93‑1 (MQ=255)
                       gTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCATCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAAt                                                                <  1:1680755/93‑1 (MQ=255)
                        ttttACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATg                                                               <  1:899198/93‑1 (MQ=255)
                                gcAAAGCCCAGTGATAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTg                                                       >  1:336978/1‑93 (MQ=255)
                                gcAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTg                                                       <  1:2878330/93‑1 (MQ=255)
                                gcAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTg                                                       <  1:1356223/93‑1 (MQ=255)
                                 cAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTgg                                                      >  2:1088175/1‑93 (MQ=255)
                                   aaGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgccgc                                                                       <  2:4719297/74‑1 (MQ=255)
                                   aaGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAAcgccgc                                                                       >  1:4719297/1‑74 (MQ=255)
                                   aaGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTc                                                    >  2:1653011/1‑93 (MQ=255)
                                    aGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCa                                                   >  2:1597418/1‑93 (MQ=255)
                                    aGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCa                                                   <  2:279230/93‑1 (MQ=255)
                                       ccAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCAtga                                                >  1:3733268/1‑91 (MQ=255)
                                       ccAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATaa                                                >  1:1048754/1‑93 (MQ=255)
                                       ccAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATaa                                                <  2:2174655/93‑1 (MQ=255)
                                       ccAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATaa                                                >  1:279027/1‑93 (MQ=255)
                                       ccAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATaa                                                >  2:839118/1‑93 (MQ=255)
                                        cAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATaaa                                               <  1:788240/93‑1 (MQ=255)
                                         aGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAg                                              >  2:3178388/1‑93 (MQ=255)
                                         aGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAg                                              >  2:4190923/1‑93 (MQ=255)
                                           tGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGcc                                            <  2:1206183/93‑1 (MQ=255)
                                             aaaaGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTg                                          >  1:2004180/1‑93 (MQ=255)
                                              aaaGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGc                                         >  1:2129601/1‑93 (MQ=255)
                                              aaaGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAACCCTGc                                         >  2:3810300/1‑93 (MQ=255)
                                               aaGCCCTTCCGTCCCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCt                                        <  2:923676/93‑1 (MQ=255)
                                               aaGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCt                                        >  1:2067325/1‑93 (MQ=255)
                                                aGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCtt                                       <  1:2385106/93‑1 (MQ=255)
                                                 gCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCttt                                      >  2:2702862/1‑93 (MQ=255)
                                                  cccTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTc                                     >  2:3065246/1‑93 (MQ=255)
                                                  cccTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTc                                     >  1:119728/1‑93 (MQ=255)
                                                   ccTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCa                                    <  2:947596/93‑1 (MQ=255)
                                                     ttACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAgc                                  >  1:4021638/1‑93 (MQ=255)
                                                     ttACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAgc                                  >  1:2243201/1‑93 (MQ=255)
                                                     ttACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAgc                                  >  2:2431138/1‑93 (MQ=255)
                                                      tACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAgcg                                 >  1:3434492/1‑93 (MQ=255)
                                                         gTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCg                              <  1:3541212/93‑1 (MQ=255)
                                                         gTACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCg                              >  1:41800/1‑93 (MQ=255)
                                                          tACCAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGt                             <  2:3353455/93‑1 (MQ=255)
                                                          gaCCCGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGt                             <  1:2404183/92‑1 (MQ=255)
                                                            ccAGCAGATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGa                           <  2:544509/93‑1 (MQ=255)
                                                             cagcagATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGAt                          <  2:277588/93‑1 (MQ=255)
                                                             cagcagATCCGCTACACAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCTCCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGAt                          <  1:565481/93‑1 (MQ=255)
                                                                cagATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTcc                       <  1:549015/93‑1 (MQ=255)
                                                                 agATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCTTGCTTTCAGCGCCGTGATTCCt                      <  2:2991050/93‑1 (MQ=255)
                                                                 agATCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCt                      >  2:1453102/1‑93 (MQ=255)
                                                                   aTCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCttt                    <  2:1417221/93‑1 (MQ=255)
                                                                   aTCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCttt                    >  1:1542555/1‑93 (MQ=255)
                                                                   aTCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCttt                    >  1:1770662/1‑93 (MQ=255)
                                                                    tCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCtttt                   <  1:2474996/93‑1 (MQ=255)
                                                                    tCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCtttt                   >  1:734734/1‑93 (MQ=255)
                                                                    tCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCtttt                   <  2:177161/93‑1 (MQ=255)
                                                                    tCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCtttt                   >  2:1885463/1‑93 (MQ=255)
                                                                    tCCGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCAGTGATTCCtttt                   >  2:2423563/1‑93 (MQ=255)
                                                                     ccGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCttttt                  <  1:4205459/93‑1 (MQ=255)
                                                                      cGCTACCCAGCCAGTGAGCGGGGCGCCGACAAACGCCGCGGGGAATTGCACCGTGGTCATAAAGCGTGCTTTCCGCGCCGGGATTCCTTTTTc                 >  2:262920/1‑93 (MQ=255)
                                                                      cGCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTc                 >  1:998115/1‑93 (MQ=255)
                                                                       gCTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTcc                >  2:1996841/1‑93 (MQ=255)
                                                                        cTACCCAGCCATAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAATCCTGCTTTCAGCGCCGTGATTCCTTTTTCCt               <  1:3336517/93‑1 (MQ=255)
                                                                        cTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCt               >  1:3076419/1‑93 (MQ=255)
                                                                        cTACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGCAATAAAGCCTGCTTTCAGCGCCGTGATTCTTTTTTCCt               >  1:1352970/1‑93 (MQ=255)
                                                                         tACCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTg              <  1:2675812/93‑1 (MQ=255)
                                                                          aCCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTgg             <  2:2115628/93‑1 (MQ=255)
                                                                          aCCCAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTgg             <  2:485983/93‑1 (MQ=255)
                                                                           cccAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTggg            <  2:4280622/93‑1 (MQ=255)
                                                                           cccAGCCAGAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTggg            >  2:902404/1‑93 (MQ=255)
                                                                            ccagccagAGAGCGGGACGCCGACAAAAGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGGGATTCCTTTTTCCTGGGt           >  2:1967664/1‑93 (MQ=255)
                                                                             cagccagAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGtt          <  1:1256092/93‑1 (MQ=255)
                                                                               gccagAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTaa        >  1:2398388/1‑93 (MQ=255)
                                                                               accagAGAGCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTaa        <  1:3509881/92‑1 (MQ=255)
                                                                                  agagagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATa     <  1:528283/93‑1 (MQ=255)
                                                                                   gagagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAg    >  1:4443706/1‑93 (MQ=255)
                                                                                   gagagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAg    >  2:1016603/1‑93 (MQ=255)
                                                                                    agagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGt   >  1:4149957/1‑93 (MQ=255)
                                                                                    agagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGt   >  1:70678/1‑93 (MQ=255)
                                                                                    agagCGGGACACCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGt   <  2:2905579/93‑1 (MQ=255)
                                                                                     gagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGtt  >  1:2624786/1‑93 (MQ=255)
                                                                                     gagCGGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGtt  >  1:2229948/1‑93 (MQ=255)
                                                                                         |                                                                                        
TCAGCAAGCCGCAGATAATCGGCGTTTTACGCGCAAAGCCCAGTGAAAAGCCCTTACGTACCAGCAGATCCGCTACCCAGCCAGAGAGCAGGACGCCGACAAACGCCGCGAGGAATGGCACCGTGGTCATAAAGCCTGCTTTCAGCGCCGTGATTCCTTTTTCCTGGGTTAAATAGTT  >  NC_000913/3868388‑3868565

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