Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,719,806 A→T 6.5% intergenic (+38/‑243) yiaG → / → cspA HTH_CROC1 family putative transcriptional regulator/RNA chaperone and antiterminator, cold‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,719,8060AT 6.5% 289.1 / 6.3 104intergenic (+38/‑243)yiaG/cspAHTH_CROC1 family putative transcriptional regulator/RNA chaperone and antiterminator, cold‑inducible
Reads supporting (aligned to +/- strand):  ref base A (47/48);  new base T (4/3);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 110 total aligned reads displayed.

CGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAT  >  NC_000913/3719657‑3719955
                                                                                                                                                     |                                                                                                                                                     
cGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCtgatga                                                                                                                                                       >  2:385312/1‑150 (MQ=255)
 gCCATGGTAAAGGAATGGGAATCCAGACGCGTGCAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGAc                                                                                                                                                      >  2:1147489/1‑150 (MQ=255)
      ggTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                                                                                 >  2:552313/1‑150 (MQ=255)
       gTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGAGGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAcc                                                                                                                                                <  1:967331/150‑1 (MQ=255)
        tcaacgAATGGGAAGCCAGACGCGTGAAGCATTCACGGGCCGAACTAATATTGATGCGTTTGATTCACGCCAACCCGGCTTTAAGTATGCAGCTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGTCATCAGGACCg                                                                                                                                               <  2:57029/145‑1 (MQ=255)
        tAAAGGAAGGGGAATCCAGACGCGTCAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTTATTCAAGCCAACCCGGCATTAAGTAAGCAGTAGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGTCATCAGGACCg                                                                                                                                               <  2:102840/150‑1 (MQ=255)
        tAAACGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCTTATGATTCAAGCCAAGCCGGCAGTAAGTAAGCAGTTGATGGAATAGACTTTTAACCACTTTATTGCTGTTTACGGTCCTGTCATCAGGACCg                                                                                                                                               <  2:20020/150‑1 (MQ=255)
           aGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGttt                                                                                                                                            <  2:1008420/150‑1 (MQ=255)
             gAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGAGGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTc                                                                                                                                          <  1:433367/150‑1 (MQ=255)
              aaTGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTcc                                                                                                                                         <  2:364019/150‑1 (MQ=255)
                   gAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCg                                                                                                                                    <  2:351735/150‑1 (MQ=255)
                     aTCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                                                                                  <  2:826209/150‑1 (MQ=255)
                     aTCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                                                                                  >  2:604409/1‑150 (MQ=255)
                       ccAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTa                                                                                                                                >  2:875494/1‑150 (MQ=255)
                          gACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATc                                                                                                                             >  1:1220682/1‑150 (MQ=255)
                           aCGCGTGAAGACTTCAAGTGACGAACTAAAATTGATGCGTTTGATTCAAGCCAACCGGGCATTAAGTAAGCATTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCa                                                                                                                            <  1:149100/150‑1 (MQ=255)
                            cgcgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCAt                                                                                                                           >  1:911514/1‑150 (MQ=255)
                             gcgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATGCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCATGATGACAGGACCGTTTTCCAACCGATTAATCATa                                                                                                                          <  1:1189063/150‑1 (MQ=255)
                                tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCACCCGATTAATCATAAat                                                                                                                       >  2:238112/1‑150 (MQ=255)
                                 gAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAata                                                                                                                      <  1:430297/150‑1 (MQ=255)
                                             tGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGCCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAtt                                                                                                          <  1:707770/150‑1 (MQ=255)
                                             tGCCGAACTAAAATTGATGCGTTTGATTCAAGACAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAtt                                                                                                          <  1:1165114/150‑1 (MQ=255)
                                              gCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttg                                                                                                         <  2:410620/150‑1 (MQ=255)
                                              gCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttg                                                                                                         >  1:1041661/1‑150 (MQ=255)
                                              gCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttg                                                                                                         <  1:896706/150‑1 (MQ=255)
                                                cGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgtt                                                                                                       >  2:1215628/1‑150 (MQ=255)
                                                   aCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCa                                                                                                    <  1:110258/150‑1 (MQ=255)
                                                   aCTAAAATTGATGCGTTAGAATCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCa                                                                                                    <  2:90556/150‑1 (MQ=255)
                                                     tAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAATTAATCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAACGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATc                                                                                                  <  1:1214940/150‑1 (MQ=255)
                                                        aaTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAATTATGAAAAATAATTGTTGCATCAcc                                                                                               >  1:501173/1‑150 (MQ=255)
                                                        aaTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAcc                                                                                               <  1:888293/150‑1 (MQ=255)
                                                        aaTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAcc                                                                                               >  1:1091196/1‑150 (MQ=255)
                                                        aaTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAcc                                                                                               >  1:505695/1‑150 (MQ=255)
                                                         aTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAccc                                                                                              <  1:1204562/150‑1 (MQ=255)
                                                          ttGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCg                                                                                             >  1:1102530/1‑150 (MQ=255)
                                                            gATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGcc                                                                                           <  2:1256073/150‑1 (MQ=255)
                                                             aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGAGGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAGCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCa                                                                                          >  2:625758/1‑150 (MQ=255)
                                                               gCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAAt                                                                                        >  1:761014/1‑150 (MQ=255)
                                                                cGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGCCAGGACCGTTTTCCAACCGAGTTATCAGAAATATGAAAAATAATTGTTGCATCACCCGCCAATg                                                                                       >  2:39858/1‑150 (MQ=255)
                                                                cGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCCCCCGCCAATg                                                                                       >  2:1146950/1‑150 (MQ=255)
                                                                cGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATg                                                                                       >  1:579698/1‑150 (MQ=255)
                                                                 gTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCAAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGc                                                                                      >  1:736978/1‑150 (MQ=255)
                                                                    tGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTg                                                                                   >  1:545254/1‑150 (MQ=255)
                                                                    tGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTg                                                                                   >  1:522555/1‑150 (MQ=255)
                                                                      aTTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTggg                                                                                 >  2:828484/1‑149 (MQ=255)
                                                                      aTTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGc                                                                                 <  1:63184/150‑1 (MQ=255)
                                                                      aTTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGc                                                                                 <  2:281925/150‑1 (MQ=255)
                                                                       ttCAAGCCAACCCGGCATTAAGTAAGCTGTTGATGGAATAGGCTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAATATAATTGTTGCATCACCCGCCATTGCGTGGCt                                                                                >  1:185495/1‑150 (MQ=255)
                                                                           aGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAAt                                                                            >  2:1067666/1‑150 (MQ=255)
                                                                           aGCCAACCAGGCATTAACTAAGCAGTTGATGGAATAGACTTTTATCAACTTTATTGATGTTTACGGTCCTGATGACAGGACCGTTTTTCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGTCTTAAt                                                                            <  1:84852/150‑1 (MQ=255)
                                                                            gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATg                                                                           >  1:285356/1‑150 (MQ=255)
                                                                              cAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGca                                                                         >  2:753611/1‑150 (MQ=255)
                                                                                  ccGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCAGATCTCAGGCCCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATCCGTGGCTAAATCACATCa                                                                    >  1:584937/1‑150 (MQ=255)
                                                                                      cATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACg                                                                 <  2:237052/150‑1 (MQ=255)
                                                                                         taagtaagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGcc                                                                                           >  1:304950/1‑121 (MQ=255)
                                                                                         taagtaagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGcc                                                                                           <  2:304950/121‑1 (MQ=255)
                                                                                         taagtaagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGtt                                                              <  2:558238/150‑1 (MQ=255)
                                                                                          aagtaagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGttt                                                             <  2:176198/150‑1 (MQ=255)
                                                                                            gtaagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGa                                                           >  1:700182/1‑150 (MQ=255)
                                                                                                 ctgTTGATGTAATATACTTTTATCCACTTTTTTGCTGATTACGTACCTGATGACAGGACCGTTTTCCAACCTATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTAc                                                      <  2:193599/148‑1 (MQ=255)
                                                                                                  aGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACa                                                     <  2:477947/150‑1 (MQ=255)
                                                                                                   gTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAg                                                    <  2:319505/150‑1 (MQ=255)
                                                                                                      tattgAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGAcc                                                 <  1:579338/146‑1 (MQ=255)
                                                                                                         ggAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCAtt                                              <  2:687533/150‑1 (MQ=255)
                                                                                                           aaTAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTTTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACAATTaa                                            >  2:432345/1‑150 (MQ=255)
                                                                                                             tAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                          <  1:375971/150‑1 (MQ=255)
                                                                                                             tAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                          <  2:730905/150‑1 (MQ=255)
                                                                                                             tAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                          <  1:14624/150‑1 (MQ=255)
                                                                                                                aCTTTTATCCACTTTATTGCTGTTTACTGTCCTGATGACATGACCGTTTTCCAACCGATTAATCATAATTATTAAAAATAATTGTTGCATCACCCTCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAg                                       <  1:813279/150‑1 (MQ=255)
                                                                                                                 cTTTTATCCCCTTTATTGCTGTTTACGCTTCTGATGACAGAAACGTTTTCCAACCGATTAATCATAAATATGAAAAATAAATGTTGCATCACCCGCCAATGCGTGGCTTAATTCAAATCAACGGTTTGACGTACAGACCATAAAAGAAgt                                      >  2:528763/1‑150 (MQ=255)
                                                                                                                 cTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAgt                                      >  1:1122159/1‑150 (MQ=255)
                                                                                                                 cTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAgtg                                     >  2:1160943/1‑150 (MQ=255)
                                                                                                                  ttttATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAgtg                                     <  1:248991/150‑1 (MQ=255)
                                                                                                                   tttATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAgtgt                                    <  1:532208/150‑1 (MQ=255)
                                                                                                                    ttATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTgta                                   >  2:566354/1‑150 (MQ=255)
                                                                                                                     tataCACTTTATTGCTGTTTTCGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTgtag                                  <  1:282129/150‑1 (MQ=255)
                                                                                                                      aTCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTgtagt                                 <  1:875494/150‑1 (MQ=255)
                                                                                                                           cTTTTTTGCTGTTTACGGTCCTGATTACAGGACCTTTTTCCCAGCGATTAATCATAAATATGAAAAATAAGAGTTGCATCACCCGCCAATGCGTGGATTAATGCAGATCAACTGTGTGACGTACCGACCATTAAAGCAGTGTAGTAAGGc                            <  2:809814/150‑1 (MQ=255)
                                                                                                                             ttATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGAAGTGTAGTAAGGCaa                          <  2:531690/150‑1 (MQ=255)
                                                                                                                              tATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAATTATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAg                         >  1:836893/1‑150 (MQ=255)
                                                                                                                               aTTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGt                        <  2:1120426/150‑1 (MQ=255)
                                                                                                                                 tGCTGTTTACGGTCCTGATGACAGTACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCAAAACAACGGTTTGACTAACAGACCATTAAAGCAGTGTAGTAAGGCAAGTcc                      >  2:476101/1‑150 (MQ=255)
                                                                                                                                  gCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTccc                     >  2:1311808/1‑150 (MQ=255)
                                                                                                                                    ggaaaaCGGTCCTGTCATCAGGACCGTTTTCCAACCGATTTATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCtt                   >  2:446515/6‑150 (MQ=255)
                                                                                                                                     gTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCCCCCGCCAATGCGTGGCTTAATGCACATCAACGTTTTGACTTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTc                  >  2:77315/1‑150 (MQ=255)
                                                                                                                                     gTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTc                  >  1:152838/1‑150 (MQ=255)
                                                                                                                                     gTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTc                  >  1:766889/1‑150 (MQ=255)
                                                                                                                                       ttACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCaa                >  2:842426/1‑150 (MQ=255)
                                                                                                                                       aaaCGGTCCTGTCATCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCaa                >  2:459642/3‑150 (MQ=255)
                                                                                                                                        tACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAag               <  2:505602/150‑1 (MQ=255)
                                                                                                                                        tACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAag               >  1:219360/1‑150 (MQ=255)
                                                                                                                                        aaCGGTCCTGTCATCAGGGCCGTTTTCCAACAGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTATTACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAag               >  1:70216/2‑150 (MQ=255)
                                                                                                                                          cGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAagag             <  2:1220682/150‑1 (MQ=255)
                                                                                                                                             tCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTa          <  1:908559/150‑1 (MQ=255)
                                                                                                                                               cTGATGACAGGACCGTTTTCCACCCGATTAATAATAAATATGAAAAATAATTGTTGCATCACCCGACAATGCGTGGCTTAATGCACATCAAAGGTTTGACGTACAGAACATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATc        >  1:984333/1‑150 (MQ=255)
                                                                                                                                                 gatgaCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGt      <  2:522555/150‑1 (MQ=255)
                                                                                                                                                   tgaCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTg    <  2:1110937/150‑1 (MQ=255)
                                                                                                                                                   tgaCAGGACCGTTTTACAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTg    <  2:1283039/150‑1 (MQ=255)
                                                                                                                                                     aCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAt  >  2:109771/1‑150 (MQ=255)
                                                                                                                                                     aCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAt  >  2:1134676/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGTTGAT  >  NC_000913/3719657‑3719955

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