Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,055,191 A→G 5.2% L85P (CTT→CCT)  ygfB ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,055,1910AG 5.2% 294.9 / 5.3 96L85P (CTT→CCT) ygfBUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base A (54/37);  new base G (3/2);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GATAGCTTCACCGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGTTACCGCCACA  >  NC_000913/3055043‑3055333
                                                                                                                                                    |                                                                                                                                              
gATAGCTTCACCGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCGAAGACCAAGCAGGAAGTGATTGACCCCACCTGCCAATGCATCAGCCCGATCGAAAACGCTGCCATCATCGCCATCAGGCAGAGAAAGCTGATAAAg                                                                                                                                               >  1:667806/1‑150 (MQ=255)
gATAGCTTCACCGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAg                                                                                                                                               >  2:692186/1‑150 (MQ=255)
     cTTCACCGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGc                                                                                                                                          <  1:5857/150‑1 (MQ=255)
      ttCACCGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGcc                                                                                                                                         >  2:14707/1‑150 (MQ=255)
          ccGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCa                                                                                                                                     >  2:1022316/1‑150 (MQ=255)
            ggTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATc                                                                                                                                   >  2:658758/1‑150 (MQ=255)
             gTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATcc                                                                                                                                  <  1:1094891/150‑1 (MQ=255)
                 cGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCa                                                                                                                              <  1:420925/150‑1 (MQ=255)
                 cGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTCCa                                                                                                                              >  1:223426/1‑150 (MQ=255)
                  gCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCCACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAg                                                                                                                             >  2:282673/1‑150 (MQ=255)
                   ccGGTCACTTTATCCAGCTTCGGTTTCGTAACGCCAAGAGCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAgg                                                                                                                            <  2:1263873/150‑1 (MQ=255)
                   ccGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATTAAGCTGAAAAAGGAAGCCGTCATCCTGCAgg                                                                                                                            >  2:460969/1‑150 (MQ=255)
                         aCTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCg                                                                                                                      >  1:303675/1‑150 (MQ=255)
                          cTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAAAGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGc                                                                                                                     >  1:285295/1‑150 (MQ=255)
                             tATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAACCGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtgg                                                                                                                  >  1:1077739/1‑150 (MQ=255)
                             tATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCCGATAAAGCTGAAAAGGGAAGCCGTCATCCTGCAGGGCATCGCtgg                                                                                                                  >  1:908866/1‑150 (MQ=255)
                             tATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAGGGAAGCCGTCATCCTGCAGGGCATCGCtgg                                                                                                                  >  1:684677/1‑150 (MQ=255)
                             tATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtgg                                                                                                                  <  2:477342/150‑1 (MQ=255)
                             tATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtgg                                                                                                                  >  1:436363/1‑150 (MQ=255)
                               tCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCACCCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtggtg                                                                                                                <  1:66101/150‑1 (MQ=255)
                               tCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtggtg                                                                                                                <  1:5510/150‑1 (MQ=255)
                               tCCAGCTTCGGGTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtggtg                                                                                                                <  1:1101560/150‑1 (MQ=255)
                                ccAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCtggtgg                                                                                                               <  2:1180200/150‑1 (MQ=255)
                                 cAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATTCATCAGCCCGATAGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAATGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGc                                                                                                              >  2:317070/1‑150 (MQ=255)
                                   gCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGACAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCag                                                                                                            <  1:128607/150‑1 (MQ=255)
                                   gCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCag                                                                                                            <  2:1216759/150‑1 (MQ=255)
                                   gCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCag                                                                                                            >  1:432216/1‑150 (MQ=255)
                                        ggTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGCCATCATCGCCAGCAGGCAGATAAAGGTGAAAAAGGAAGCCGTCATCCTGCAGGGCATGGCTGGTGGCAGAGTGc                                                                                                       >  1:254578/1‑150 (MQ=255)
                                            gCGTAACGCCAAGACCAAGCAGGAAGTGATTGATCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCAttt                                                                                                   <  2:683152/150‑1 (MQ=255)
                                            gCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCAttt                                                                                                   <  1:986609/150‑1 (MQ=255)
                                            gCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCAttt                                                                                                   <  1:1271999/150‑1 (MQ=255)
                                            gCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCGGGTGGCAGAGTGCAttt                                                                                                   <  1:594897/150‑1 (MQ=255)
                                                 aCGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGc                                                                                              <  1:407873/150‑1 (MQ=255)
                                                    ccAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAGGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGt                                                                                           >  1:422208/1‑150 (MQ=255)
                                                    ccAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGt                                                                                           >  1:24654/1‑150 (MQ=255)
                                                    ccAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGt                                                                                           >  2:953391/1‑150 (MQ=255)
                                                         aCCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCtg                                                                                      <  1:976026/150‑1 (MQ=255)
                                                          ccAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCtgt                                                                                     >  2:310856/1‑150 (MQ=255)
                                                          ccAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCtgt                                                                                     >  2:755729/1‑150 (MQ=255)
                                                           cAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCtgtg                                                                                    >  1:420485/1‑150 (MQ=255)
                                                             aGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGcc                                                                                  <  1:366199/150‑1 (MQ=255)
                                                               cAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAACACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCGGTGCCAg                                                                                >  1:1123545/1‑150 (MQ=255)
                                                                     gTGATTGACCCATCCTGCGAATGCATCAGCCCGATCGAAAACGCTGACATCTTCGCCTTCAGGCAGATATAGCTGAAAAAGGAAGCCGTCATCCTGGAGGGCATCGCTGGTGTCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATg                                                                          <  1:1281223/150‑1 (MQ=255)
                                                                      tGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGa                                                                         >  1:827077/1‑150 (MQ=255)
                                                                       gATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAAGGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGAc                                                                        <  1:181263/150‑1 (MQ=255)
                                                                           gACCCAACCTGCGAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATAAGGCAGCTAAAGCTGAACAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCCGCTCATGCCCGaa                                                                    >  1:830236/1‑150 (MQ=255)
                                                                            aCCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGaaa                                                                   >  2:823479/1‑150 (MQ=255)
                                                                             cccACCCTGCCAAGGCATCAGCCCGATCGAAAACGCTGCCATCATCGCCATCAGGCAGATAACGCGGAAAAAGGAAGCCGTCCTCCTGCAGGGCATCGCTGGTGGCAGGGTGCCTGTTACGCAGTGCCTGTGCCAGCTCCTGACCGacac                                                                  >  1:737315/1‑147 (MQ=255)
                                                                             cccAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAg                                                                  >  1:581735/1‑150 (MQ=255)
                                                                                  ccTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCCTCCTGCAGGGCATCGCTGGTGGAAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATg                                                             >  2:1023776/1‑150 (MQ=255)
                                                                                  ccTGCCAATGCAGCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATg                                                             <  1:843628/150‑1 (MQ=255)
                                                                                    tGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGcc                                                           <  1:1280717/150‑1 (MQ=255)
                                                                                     gCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCt                                                          >  1:952785/1‑150 (MQ=255)
                                                                                      ccAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACAGAAAGCCATGCCtt                                                         >  1:1019246/1‑150 (MQ=255)
                                                                                      ccAATGCATCAGACCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGCCCGATAGCCATGCCtt                                                         >  1:1026353/1‑150 (MQ=255)
                                                                                       cAATGCATCAGACCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGCAAGCCGTCATCCTGCAGGGCATCGCTGCTGGCAGAGTGCACTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttc                                                        <  2:1186189/150‑1 (MQ=255)
                                                                                         aTGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAGGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttcgt                                                      <  2:684677/150‑1 (MQ=255)
                                                                                         aTGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAGGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttcgt                                                      <  2:908866/150‑1 (MQ=255)
                                                                                         aTGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttcgt                                                      <  2:1020892/150‑1 (MQ=255)
                                                                                         aTGCATCAGCCCGATCGAAAACGCTGACATAATCTCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttcgt                                                      <  2:94491/150‑1 (MQ=255)
                                                                                          tGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGTGCATCGCTGGTGGCAGAGTGAATTTTACGAAGTGCCTGTGCCAGCTCATGACCGAAAGCAATGCCttcgtt                                                     >  2:269743/1‑150 (MQ=255)
                                                                                          tGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttcgtt                                                     >  1:1285524/1‑150 (MQ=255)
                                                                                            cATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCttcgttcg                                                   <  1:777613/150‑1 (MQ=255)
                                                                                                   ccGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTc                                            >  2:329419/1‑150 (MQ=255)
                                                                                                     gATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGt                                          >  2:807022/1‑150 (MQ=255)
                                                                                                       tCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCGTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCTGCTCATGACCGAAAGCCAAGCCTTCGTTCGTCAGGTCGTGa                                        <  2:351137/150‑1 (MQ=255)
                                                                                                       tCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGa                                        <  1:1022316/150‑1 (MQ=255)
                                                                                                        cGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGaa                                       <  1:800551/150‑1 (MQ=255)
                                                                                                         gAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAg                                      >  2:440212/1‑150 (MQ=255)
                                                                                                          aaaaCGCTGACATCATCGCCATCAGGCAGATAAAGCTGAATAAGTAATCCGTCATCCTGCAGGGCATCCCTGGTGTAATAGTGAATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCGATGCCATCGTTCTTCAGGTCGTGAAga                                     >  2:130012/1‑149 (MQ=255)
                                                                                                               gCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCgg                                >  1:1026513/1‑150 (MQ=255)
                                                                                                               gCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCgg                                <  1:451978/150‑1 (MQ=255)
                                                                                                                  gACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGATGCCGTCATCCGGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAg                             >  2:50589/1‑150 (MQ=255)
                                                                                                                  gACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAg                             >  2:1066201/1‑150 (MQ=255)
                                                                                                                  gACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAg                             >  2:1353405/1‑150 (MQ=255)
                                                                                                                  gACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAg                             >  2:673797/1‑150 (MQ=255)
                                                                                                                   aCATCATCGCCATCAGGCTGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGTCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGc                            <  1:18356/150‑1 (MQ=255)
                                                                                                                   aCATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGc                            <  2:1221854/150‑1 (MQ=255)
                                                                                                                    catcatCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGcc                           >  2:22146/1‑150 (MQ=255)
                                                                                                                     atcatcGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTAATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTATCGGTAGCCa                          >  2:572110/1‑150 (MQ=255)
                                                                                                                       tatcGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATg                        <  1:1113982/149‑1 (MQ=255)
                                                                                                                       catcGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGTATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATg                        >  2:17318/1‑150 (MQ=255)
                                                                                                                       catcGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGTATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATg                        >  2:17319/1‑150 (MQ=255)
                                                                                                                       catcGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCTTTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATg                        <  1:591918/150‑1 (MQ=255)
                                                                                                                       catcGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATg                        >  1:487936/1‑150 (MQ=255)
                                                                                                                            ccATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGCTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGATCATGACCGAAAGCCATGCCTTCGTTAGTCAGGTCGTCAAGTAGCGGTAGCCATGAGCTg                   >  2:1201166/1‑150 (MQ=255)
                                                                                                                             cATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGt                  <  1:692186/150‑1 (MQ=255)
                                                                                                                                    cAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGtt           >  1:103288/1‑150 (MQ=255)
                                                                                                                                     aGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGAGCCTGTGCCAGCTCATGACCGAAAGCCATGCCGTGGTACGTCAGGTCGTCAAGTAGCGGTAGCCATGAGCTGTCAGCGTTa          >  1:808884/1‑150 (MQ=255)
                                                                                                                                     aGATAAAGCTGAAAAAGGAAGCCGACATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCCTTa          >  2:1315641/1‑150 (MQ=255)
                                                                                                                                       atatAGCTGTAAAAGGAAGCCGTCATCCTGCAGGGCAGCGCTGGTGGCATAGTGCTTTTTACGCAGTGCCTGTTCCAGCTCATGACCGATAGCCAAGCCTTCGTTCTTCAGGTCGTGAAGTAGCGGTAGCCATGTGCTGTCATCGTTAcc        <  1:1225565/150‑1 (MQ=255)
                                                                                                                                        tAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGTTACCg       >  2:1022738/1‑150 (MQ=255)
                                                                                                                                          aaGCTGAAAAAGGAAGCCGTCCTCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCGTCGTTCGTCAGGTCGTGAAGTAGCGGTTGCCATGCGCTGTCATCGTTACCGcc     >  1:146237/1‑150 (MQ=255)
                                                                                                                                           aGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCCGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGTTACCGccc    >  1:1200636/1‑149 (MQ=255)
                                                                                                                                            gCTGAATAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGTTACCGCcac   <  1:523286/150‑1 (MQ=255)
                                                                                                                                             cTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGTTACCGCcaca  >  2:1006375/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GATAGCTTCACCGGTTTCGCCGGTCACTTTATCCAGCTTCGGTTGCGTAACGCCAAGACCAAGCAGGAAGTGATTGACCCAACCTGCCAATGCATCAGCCCGATCGAAAACGCTGACATCATCGCCATCAGGCAGATAAAGCTGAAAAAGGAAGCCGTCATCCTGCAGGGCATCGCTGGTGGCAGAGTGCATTTTACGCAGTGCCTGTGCCAGCTCATGACCGAAAGCCATGCCTTCGTTCGTCAGGTCGTGAAGTAGCGGTAGCCATGAGCTGTCATCGTTACCGCCACA  >  NC_000913/3055043‑3055333

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: