Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,901,405 C→T L272L (CTG→CTA bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,4050CT96.8% 334.2 / 0.7 95L272L (CTG→CTAbglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (53/39);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

GTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTGGC  >  NC_000913/3901261‑3901552
                                                                                                                                                |                                                                                                                                                   
gTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGcgcg                                                                                                                                               <  2:940532/151‑1 (MQ=255)
 tttGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGcgcg                                                                                                                                               >  2:199342/1‑150 (MQ=255)
    gCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGatat                                                                                                                                           <  1:110654/151‑1 (MQ=255)
    gCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGatat                                                                                                                                           >  2:451411/1‑151 (MQ=255)
      cGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATc                                                                                                                                          >  2:424323/1‑150 (MQ=255)
      cGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCa                                                                                                                                         <  2:1046583/151‑1 (MQ=255)
         gACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAAt                                                                                                                                       <  1:133355/150‑1 (MQ=255)
          aCATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAAttt                                                                                                                                     >  2:5102/1‑151 (MQ=255)
          aCATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAAttt                                                                                                                                     <  1:427965/151‑1 (MQ=255)
            aTAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTGTTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGAAATCTTCGCGAACTAGCGCGATATCAAtttt                                                                                                                                    >  1:477919/1‑150 (MQ=255)
                cGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCg                                                                                                                               >  2:514614/1‑151 (MQ=255)
                 gACCACTTTCCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGCTCGTAAACATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCgg                                                                                                                              <  1:398570/151‑1 (MQ=255)
                 gACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCgg                                                                                                                              <  1:540121/151‑1 (MQ=255)
                 gACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCgg                                                                                                                              >  2:775298/1‑151 (MQ=255)
                     aCTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCg                                                                                                                          <  2:1097987/151‑1 (MQ=255)
                         aCCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACt                                                                                                                       <  2:319547/150‑1 (MQ=255)
                          ccATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACttt                                                                                                                     >  1:234640/1‑151 (MQ=255)
                                 ggTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGttt                                                                                                              <  1:807108/151‑1 (MQ=255)
                                   tGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTc                                                                                                            >  2:821596/1‑151 (MQ=255)
                                    ggCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTcc                                                                                                           >  1:22838/1‑151 (MQ=255)
                                    ggCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTcc                                                                                                           >  2:616303/1‑151 (MQ=255)
                                       tttATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTc                                                                                                            >  1:576819/1‑147 (MQ=255)
                                            cccATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTa                                                                                                    >  1:1315383/1‑150 (MQ=255)
                                             ccATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTac                                                                                                   >  2:212065/1‑150 (MQ=255)
                                              cATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCGTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTacac                                                                                                 <  1:1326922/151‑1 (MQ=255)
                                               aTAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACAcc                                                                                                <  2:460716/151‑1 (MQ=255)
                                                 aaCGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACAccc                                                                                               <  1:517227/150‑1 (MQ=255)
                                                   cGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCt                                                                                            <  1:395685/151‑1 (MQ=255)
                                                   cGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCt                                                                                            >  2:584438/1‑151 (MQ=255)
                                                    gggATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTg                                                                                           >  2:1005766/1‑151 (MQ=255)
                                                    gggATATCTTTATAGCGGAAATCCATGGTAGTGGTATTAATCTGCTGTTTGTTTTGCATGATGAGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACTCGCTg                                                                                           >  1:1150396/1‑151 (MQ=255)
                                                      gATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCg                                                                                         <  2:700814/151‑1 (MQ=255)
                                                      gATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCg                                                                                         <  2:697819/151‑1 (MQ=255)
                                                        tatCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGc                                                                                       >  2:293145/1‑151 (MQ=255)
                                                         atCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGc                                                                                       >  1:1246386/1‑150 (MQ=255)
                                                         agcTTTATAGCGTAAATGATTGGTATTGGGATTTGTCGGCTGTTTGTGTTGCAGGCTGGGGTCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAGTCTTACGCGGACTGACTTTCCAGGTTTCCAGCCCTACACCCGCTGc                                                                                          >  1:82046/3‑147 (MQ=255)
                                                           cTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                                                    >  1:1138447/1‑151 (MQ=255)
                                                            tttATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                                                    >  2:932980/1‑150 (MQ=255)
                                                            tttATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                                                    >  2:414061/1‑150 (MQ=255)
                                                             ttATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                                                    <  1:1306773/149‑1 (MQ=255)
                                                               ataGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAg                                                                                  >  1:384329/1‑149 (MQ=255)
                                                                  gCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGc                                                                             >  2:1041254/1‑151 (MQ=255)
                                                                  gCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGGTTTGCAGGCTGCGCTCGTAATCCTCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCCGCATCAGTat                                                                               <  1:706936/149‑2 (MQ=255)
                                                                     tAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGCTCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGc                                                                          >  1:521769/1‑151 (MQ=255)
                                                                      aaaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCg                                                                         >  2:570497/1‑151 (MQ=255)
                                                                       aaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGt                                                                        >  1:788900/1‑151 (MQ=255)
                                                                        aTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGtt                                                                       <  1:390306/151‑1 (MQ=255)
                                                                         tCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGttt                                                                      >  1:79568/1‑151 (MQ=255)
                                                                         tCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGttt                                                                      >  2:341410/1‑151 (MQ=255)
                                                                          cAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGtttt                                                                     >  2:673437/1‑151 (MQ=255)
                                                                           aaTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTc                                                                    >  2:1082474/1‑151 (MQ=255)
                                                                           aaTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTc                                                                    >  2:69354/1‑151 (MQ=255)
                                                                           aaTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTGTACCCGGACCGACTTTCCAGTTTTCCAGGCCTTCACCCGCTGCGGCATCAGTGCGCTGCGTTTTc                                                                    <  1:1198049/151‑1 (MQ=255)
                                                                            aTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTcc                                                                   >  1:794846/1‑151 (MQ=255)
                                                                              ggtattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAg                                                                 <  2:1411456/151‑1 (MQ=255)
                                                                               gtattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACTCTGCGTTTTCCAGt                                                                <  1:417306/151‑1 (MQ=255)
                                                                               gtattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGt                                                                >  1:858637/1‑151 (MQ=255)
                                                                                tattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTc                                                               >  1:1024797/1‑151 (MQ=255)
                                                                                tattggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTc                                                               >  2:983146/1‑151 (MQ=255)
                                                                                 attggcattAATCTGCTGTTTGTTTTTCAGGCCGGGATCGTAATCATCGACATCTTCGCGAACTAGCGCGCAATCAATTTTACCCGGACCGAGTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGt                                                                >  1:731611/1‑149 (MQ=255)
                                                                                   tggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCaa                                                             <  1:344116/150‑1 (MQ=255)
                                                                                    ggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGc                                                           >  1:595438/1‑151 (MQ=255)
                                                                                       attAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATc                                                        >  1:1192239/1‑151 (MQ=255)
                                                                                         tAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCt                                                       <  1:858991/150‑1 (MQ=255)
                                                                                          aaTCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTgg                                                     >  1:800426/1‑151 (MQ=255)
                                                                                           aTCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTgg                                                     >  2:1350958/1‑150 (MQ=255)
                                                                                              tgctgTTTTTTTTGCAGGCTGCGACCGTATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGa                                                    <  2:357576/147‑1 (MQ=255)
                                                                                              tgctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAttt                                                 >  2:364499/1‑151 (MQ=255)
                                                                                              tgctgTTTGTTTTGCACGCTGCGATCGTAATCATCGACATCTTCGCGACCTAGCGCGATATCACTTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCCAGCATCTGGCttt                                                 >  1:766456/1‑151 (MQ=255)
                                                                                               gctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAttt                                                 <  2:991980/150‑1 (MQ=255)
                                                                                               gctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTc                                                <  2:20240/151‑1 (MQ=255)
                                                                                               actgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTc                                                >  1:648130/2‑151 (MQ=255)
                                                                                                       ttttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCg                                        >  2:741991/1‑151 (MQ=255)
                                                                                                         ttGCAGGCTGCGATCGTAATCATCGACAGCTTCGCGAACTAGCGCGATATCAATTTTCCCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTTCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCg                                        <  2:1221573/149‑1 (MQ=255)
                                                                                                           gCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgcgc                                    >  2:238333/1‑151 (MQ=255)
                                                                                                             aGGCTGCTATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCCGACCGACTTTCCAGTTTTCGAGACCTACTCCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgc                                      <  1:1197456/147‑1 (MQ=255)
                                                                                                             aGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCAc                                  <  1:1182451/151‑1 (MQ=255)
                                                                                                                cTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCAcc                                 >  2:977606/1‑149 (MQ=255)
                                                                                                                 tGCGATCGTAATCATCGACCGCTTCGCTAACTAGCGCGATATCAATTTGTCCCGGGGCGGCTTTCCAGTTTTCCAGACCCACCCCCGCTGCGGCCTTAGTACTCTGCGTTTTCCAGTCTAGCGTCTGGATTTCAATTTTCGGCGCCCCg                                <  1:251140/149‑1 (MQ=255)
                                                                                                                 tGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCgtg                              >  2:608854/1‑151 (MQ=255)
                                                                                                                      tCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTa                          <  1:1345875/150‑1 (MQ=255)
                                                                                                                          aaTCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTaccca                      <  2:452170/150‑1 (MQ=255)
                                                                                                                             catcGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTacccaccc                   <  1:921506/150‑1 (MQ=255)
                                                                                                                             catcGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTacccaccc                   >  2:362015/1‑150 (MQ=255)
                                                                                                                             catcGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTAcccaccca                  <  1:199333/151‑1 (MQ=255)
                                                                                                                                cgccATCTTCGCGAACTAGCGCGATCTCAATTTTACCCGGCCCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTTCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGa                <  1:615328/147‑1 (MQ=255)
                                                                                                                                  aCATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAAt             >  1:979019/1‑151 (MQ=255)
                                                                                                                                   cATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATc            <  2:750873/151‑1 (MQ=255)
                                                                                                                                      cTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTCACCCACCCAGAAATCAGc         >  2:323104/1‑151 (MQ=255)
                                                                                                                                        tCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGc         <  1:1022452/149‑1 (MQ=255)
                                                                                                                                         cgcgATCTACCGCGATCTCAATTTTACCCGGACCGACATCCCATTTTTCCAGACCGACTCCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCAACTGGATTTCCATTTTCGACGCACCGTGTTTACCCACCCAGAAATCAGcc        <  1:56204/149‑1 (MQ=255)
                                                                                                                                           cgAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTg    >  2:341979/1‑151 (MQ=255)
                                                                                                                                            gAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTgg   <  1:380661/151‑1 (MQ=255)
                                                                                                                                            gAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTgg   <  1:656474/151‑1 (MQ=255)
                                                                                                                                             aaCTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTGgc  >  2:800653/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTGGC  >  NC_000913/3901261‑3901552

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