Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,094,823 A→C V387G (GTC→GGC)  rhaA ← L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,8230AC98.8% 298.2 / ‑3.5 82V387G (GTC→GGC) rhaAL‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (45/36);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

TGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCATATTGCGTGTACCAATG  >  NC_000913/4094675‑4094966
                                                                                                                                                    |                                                                                                                                               
tGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGccc                                                                                                                                               >  2:796994/1‑151 (MQ=255)
 gCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTGGGCAATACATTTCCCAGCCCg                                                                                                                                              >  2:14235/1‑151 (MQ=255)
  ccGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGc                                                                                                                                             >  2:30437/1‑151 (MQ=255)
  ccGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGc                                                                                                                                             >  2:1163330/1‑151 (MQ=255)
   cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGcc                                                                                                                                            >  1:632618/1‑151 (MQ=255)
    gACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCt                                                                                                                                           <  1:1060920/151‑1 (MQ=255)
    gACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCt                                                                                                                                           <  2:747348/151‑1 (MQ=255)
      cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGc                                                                                                                                         >  2:535207/1‑151 (MQ=255)
        cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGcc                                                                                                                                        >  2:1191923/1‑150 (MQ=255)
        cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCa                                                                                                                                       <  1:842379/151‑1 (MQ=255)
         gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCAc                                                                                                                                      <  1:1079624/151‑1 (MQ=255)
                ggATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAAc                                                                                                                                >  2:145351/1‑150 (MQ=255)
                ggATAAGGCCCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAAc                                                                                                                                <  2:6214/150‑1 (MQ=255)
                  aTAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGAt                                                                                                                             <  2:651198/151‑1 (MQ=255)
                  aTAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGAt                                                                                                                             >  2:869920/1‑151 (MQ=255)
                    aaGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGAttt                                                                                                                           >  2:872030/1‑151 (MQ=255)
                        cgcTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGc                                                                                                                       >  2:326939/1‑151 (MQ=255)
                         gcTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGct                                                                                                                      >  2:262695/1‑151 (MQ=255)
                         gcTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGct                                                                                                                      >  1:825514/1‑151 (MQ=255)
                         gcTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGct                                                                                                                      >  1:1089975/1‑151 (MQ=255)
                         gcTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGct                                                                                                                      >  1:642989/1‑151 (MQ=255)
                           tCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGctct                                                                                                                    >  2:794007/1‑151 (MQ=255)
                              cgccgcATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTcc                                                                                                                 >  2:735417/1‑151 (MQ=255)
                                 cgcATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCcagc                                                                                                              >  1:1133192/1‑151 (MQ=255)
                                  gcATCCGGCAGTGTTTCCCCGCGGCGACTCAAAATTTCTTTCTCAAAAGCCCGCGCGCTCTCAAGCCATTCGCGACCTGCTGGCGTATCGGGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCcagca                                                                                                             <  2:367486/151‑1 (MQ=255)
                                      ccGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCCCGCTCTCCAGCCATTCGCTTCCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGc                                                                                                         <  1:43755/151‑1 (MQ=255)
                                      ccGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGc                                                                                                         >  2:590167/1‑151 (MQ=255)
                                         gCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAg                                                                                                      <  1:167527/151‑1 (MQ=255)
                                              gtTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGAcgcg                                                                                                 >  2:68870/1‑151 (MQ=255)
                                              gtTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGAcgc                                                                                                  <  1:578261/150‑1 (MQ=255)
                                                 tACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACgcgc                                                                                                >  2:622926/1‑149 (MQ=255)
                                                     cgcgGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTa                                                                                          <  1:4618/151‑1 (MQ=255)
                                                           gccTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTTCCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATcg                                                                                      <  2:332419/147‑1 (MQ=255)
                                                           gACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATcgcc                                                                                    >  2:208076/1‑151 (MQ=255)
                                                             cTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATcgcccg                                                                                  >  2:404169/1‑151 (MQ=255)
                                                              tCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATcgcccgc                                                                                 >  2:396025/1‑151 (MQ=255)
                                                                aaaaTTTCTTCCCCATAAGCCCGCACGCTCTCCAGCCATTCGCTCCCGGCTGGCTTATCGTGACGTTGGCAATACATTCCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCccgccg                                                                               <  1:321399/151‑1 (MQ=255)
                                                                 aaaTTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCcgccgc                                                                              >  2:1007916/1‑151 (MQ=255)
                                                                      tctttctCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTcc                                                                          <  1:269778/150‑1 (MQ=255)
                                                                      tctttctCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTcc                                                                          >  2:769118/1‑150 (MQ=255)
                                                                      tctttctCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTcc                                                                          >  2:387243/1‑150 (MQ=255)
                                                                         ttccCAAAAGCCCGCACGCTCTCCAGCCATTCGCTCCCGGCGGGCTTATCGTGCCGTGGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAg                                                                        <  2:1206459/149‑1 (MQ=255)
                                                                            tcATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTgc                                                                   >  2:1114097/1‑151 (MQ=255)
                                                                               tAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTgc                                                                   <  1:1015644/148‑1 (MQ=255)
                                                                                aaGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCa                                                                <  1:398638/150‑1 (MQ=255)
                                                                                 aGCCCGCACGCTCTCCAGCCATTCGCTCCCTGCTGGCGTATCGTGCCGTTGCCATTACATTCCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCa                                                                <  2:315570/149‑1 (MQ=255)
                                                                                  gCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCt                                                             >  2:1176374/1‑151 (MQ=255)
                                                                                    ccGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCAAGCACGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTACTCGCCCGCCGCTTCCAGCTTGCGCAGCTc                                                            >  1:681266/1‑150 (MQ=255)
                                                                                     cGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAg                                                          <  1:163145/151‑1 (MQ=255)
                                                                                      gCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGc                                                         >  2:205851/1‑151 (MQ=255)
                                                                                      gCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGc                                                         >  2:879828/1‑151 (MQ=255)
                                                                                        aCGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCgg                                                       >  1:784675/1‑151 (MQ=255)
                                                                                         cGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCgg                                                       >  1:659066/1‑150 (MQ=255)
                                                                                              tcCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGtt                                                 >  2:119154/1‑151 (MQ=255)
                                                                                               ccACCCATTCGCTACCTGCTGGCGTATCTTGACGTGGGCAATACATTCCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAg                                                    <  2:1022442/147‑1 (MQ=255)
                                                                                                cAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTcc                                               >  1:717268/1‑151 (MQ=255)
                                                                                                 aGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTCCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCa                                              <  1:30438/151‑1 (MQ=255)
                                                                                                   ccATTCGCTCCCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCa                                              <  1:481399/149‑1 (MQ=255)
                                                                                                   ccATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGAAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTcc                                               <  1:581899/148‑1 (MQ=255)
                                                                                                    cATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGc                                            >  1:711007/1‑150 (MQ=255)
                                                                                                    cATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGc                                            >  1:423534/1‑150 (MQ=255)
                                                                                                     aTTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCa                                           >  2:471517/1‑150 (MQ=255)
                                                                                                        cgttACCTGCGGGCGTATCGTGACGTTGGCAATACTTTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAacgc                                       <  2:432759/148‑1 (MQ=255)
                                                                                                        cGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAac                                         <  2:422192/149‑1 (MQ=255)
                                                                                                            cccggcgggggtttCGTGCCGTTGGCAATCCATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAacgcacgc                                   <  2:443826/138‑1 (MQ=255)
                                                                                                             ccTGCTGGCTTTTCGTGAGGTTGGCAATACATTCCCCACCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAacgcac                                     <  2:468599/148‑1 (MQ=255)
                                                                                                             ccTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAacgcacgc                                   <  1:1159421/150‑1 (MQ=255)
                                                                                                             ccTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAacgcacgc                                   <  1:1114456/150‑1 (MQ=255)
                                                                                                                gctgGCGTTCCGTGACTTTGGCATTCCTTTCCCCACCCCGCCTGCCACGGCAACGATTTCTGCTCTCCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAacgcac                                     <  2:891625/145‑1 (MQ=255)
                                                                                                                   ggCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGcagcag                              <  1:597642/149‑1 (MQ=255)
                                                                                                                   ggCGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAggg                            >  2:74577/1‑151 (MQ=255)
                                                                                                                   ggCGTAGTGAGACGTTGGCAATCCTTTTCCCACCCCGCCGGCCACGGCAAGGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAgg                             <  2:936596/150‑1 (MQ=255)
                                                                                                                     cGTATCGTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCt                          <  1:96460/151‑1 (MQ=255)
                                                                                                                       tgtcggcaCGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCGTCGATTTCTGCTCTACCAGCGGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCtt                         <  1:735453/143‑1 (MQ=255)
                                                                                                                          cgggACGTTGGCAATCCATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCttttt                      <  2:564552/147‑1 (MQ=255)
                                                                                                                           gTGACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCtttttt                     >  2:1122978/1‑150 (MQ=255)
                                                                                                                            tGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTc                    >  1:255911/1‑150 (MQ=255)
                                                                                                                             gACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCat                  >  1:1071407/1‑151 (MQ=255)
                                                                                                                             gACGTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCat                  >  1:54887/1‑151 (MQ=255)
                                                                                                                                gTTGGCAATACATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCATAtt               >  1:880861/1‑151 (MQ=255)
                                                                                                                                          cATTTCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCATATTGCGTGTAcc      <  1:626155/150‑1 (MQ=255)
                                                                                                                                             ttCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCATATTGCGTGTACCAATg  >  1:154319/1‑151 (MQ=255)
                                                                                                                                              tCCCAGCCCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCATATTGCGTGTACCAATg  >  2:356863/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTTCTGCTCTTCCAGCAGTGCCAGACGCGCAGTGTAATCGCCCGCCGCTTCCAGCTTGCGCAGCTCAGCGGTAGGTTCCAGCAACGCACGCAGCAGGGCTTTTTTCATATTGCGTGTACCAATG  >  NC_000913/4094675‑4094966

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