Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,400,330 C→G 5.5% L15V (CTG→GTG)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,400,3300CG 5.5% 236.1 / 5.0 85L15V (CTG→GTG) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (43/37);  new base G (3/2);  total (46/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 = 5.88e-01

TGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCG  >  NC_000913/4400181‑4400479
                                                                                                                                                     |                                                                                                                                                      
tgttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGATAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAc                                                                                                                                                        >  2:144428/1‑150 (MQ=255)
   tgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGTGGCAATCTTTACAAGATACGTTCCTGAACGCACTGc                                                                                                                                                     >  1:1034119/1‑150 (MQ=255)
   tgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGc                                                                                                                                                     >  2:1046781/1‑150 (MQ=255)
       gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGATAAGGGGCAATCTTTACAAGAACCGTTACTGAACGAACTGCG‑TCg                                                                                                                                                >  2:1044155/1‑150 (MQ=255)
        cTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGTAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCgg                                                                                                                                               >  1:213138/1‑150 (MQ=255)
        cTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCgg                                                                                                                                               >  2:631863/1‑150 (MQ=255)
        cTATCGCAGGCTGAATGTGTACAATTGACACGTATCGAGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATCTAAGGAAACGAGAGAATGGCTAAGGGGCAATCTTAACAAGATCCGTTCCTGAACGCACTGCG‑TCgg                                                                                                                                               <  1:52090/150‑1 (MQ=255)
         tATCGCAGGCTGAATGTTTAGAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTAAAACTACAAATAAGCATATAAGGAAAAGAGCGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCggg                                                                                                                                              <  2:92044/150‑1 (MQ=255)
         tATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCggg                                                                                                                                              >  2:780585/1‑150 (MQ=255)
          aTCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTCTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGa                                                                                                                                             >  2:66891/1‑150 (MQ=255)
             gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACg                                                                                                                                          <  2:565916/150‑1 (MQ=255)
             gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACg                                                                                                                                          >  1:1292768/1‑150 (MQ=255)
              cAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCCG‑ACGCAGTGCGTTCAGGAACgg                                                                                                                                         <  1:393861/150‑2 (MQ=255)
                ggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACgtgt                                                                                                                                       >  1:1227714/1‑150 (MQ=255)
                   tGAATGTGTACAATTGAGACGTAGCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAACGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTcc                                                                                                                                    <  1:1178667/150‑1 (MQ=255)
                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCt                                                                                                                             >  2:1297130/1‑150 (MQ=255)
                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCt                                                                                                                             >  1:778809/1‑150 (MQ=255)
                             cAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatt                                                                                                                          <  2:987010/150‑1 (MQ=255)
                                  gagaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTatttattt                                                                                                                     >  1:318842/1‑150 (MQ=255)
                                  gagaCCTATCGTGCGCAATTATTTCAGAATCGAAAGGGTCAAAGTACAAATAAGCATATAAGCAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTatttattt                                                                                                                     <  1:78731/150‑1 (MQ=255)
                                     aCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                                                                                                  <  2:577450/150‑1 (MQ=255)
                                      cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                                                                                                 >  1:1234744/1‑150 (MQ=255)
                                      cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                                                                                                 >  1:1229863/1‑150 (MQ=255)
                                      cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                                                                                                 <  1:35607/150‑1 (MQ=255)
                                           cGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTCCAAATAAGCATATAAGGAAAAGAGCGACTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGGGTTCCAGTTTCTATTTATTTGGTGAATgg                                                                                                            <  2:49643/150‑1 (MQ=255)
                                           cGTGCGCAATTTTTTCAGAATCGAAAGGATCAAAGTACAAATAAGCATATAAGGAAATGAGAGAATGGCTAAGGGGCAATCTTTACAACATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAttg                                                                                                            >  1:1023502/1‑148 (MQ=255)
                                             tGCGCAATTTTTTCAGAATCGAAAGGTTCAATGTACAAATACGCATATTAGGAAATGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAGAGTGATCCAGTTTCTATTTATTTGGGGTATGGTa                                                                                                          >  1:620741/1‑150 (MQ=255)
                                             tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                                                                                          >  1:438187/1‑150 (MQ=255)
                                               cgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAtt                                                                                                        >  1:766034/1‑150 (MQ=255)
                                                gcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                                                                                       >  2:213110/1‑150 (MQ=255)
                                                  aaTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCATGAACGAACTGCG‑TCGGGAACGTGTTCCAGATTCTATTTATTTGGTGAATGGTATTAAg                                                                                                     >  1:42400/1‑150 (MQ=255)
                                                     tttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                                                                                                  <  2:287378/150‑1 (MQ=255)
                                                       tttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                                                                                                <  2:238150/150‑1 (MQ=255)
                                                         tCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGATGCAAg                                                                                              >  1:358547/1‑150 (MQ=255)
                                                           aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                            >  2:842422/1‑150 (MQ=255)
                                                           aGAATCGAAAGGTTCAAAGTACAAATAAGCATAT‑AGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                           <  1:177910/150‑1 (MQ=255)
                                                              aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaa                                                                                         >  2:1077463/1‑150 (MQ=255)
                                                                 gAAAGGTTCAAAGTTCAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTTAATGGTATTAAGCTGCAAGGGCAAATc                                                                                      >  1:854383/1‑150 (MQ=255)
                                                                 gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGTAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc                                                                                      <  2:484960/150‑1 (MQ=255)
                                                                     ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCgtgt                                                                                  >  2:122913/1‑150 (MQ=255)
                                                                     ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGt                                                                                  <  2:809866/150‑1 (MQ=255)
                                                                         cAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCttt                                                                              <  1:1204357/150‑1 (MQ=255)
                                                                          aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGACCGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtttt                                                                             >  2:1251607/1‑150 (MQ=255)
                                                                          aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGAAAGGGCAAATCGAGTCtttt                                                                             <  1:313539/150‑1 (MQ=255)
                                                                           aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAAGGTGTTCCAGTTTCTAATTATTTGGTGAATGGGAGTAAGCTGCAAGGGCAAATCGAGTCTTTTg                                                                            >  2:182108/1‑150 (MQ=255)
                                                                            aGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAATTCGAGTCTTTTGa                                                                           >  1:544067/1‑150 (MQ=255)
                                                                            aGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCAGTTAATGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                           <  2:236915/150‑1 (MQ=255)
                                                                            aGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCGAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCG‑TCGGGAACCTGATCCAGTTTCTATATATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                           >  2:105684/1‑150 (MQ=255)
                                                                              tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCGCTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                         <  2:1186238/150‑1 (MQ=255)
                                                                                  aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGtt                                                                     >  2:836984/1‑150 (MQ=255)
                                                                                    tAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCg                                                                   <  2:1114951/150‑1 (MQ=255)
                                                                                      aGCATATAAGGAAAAGAGAGAATGGCTAAGGTGCAATCTTTACAAGAACCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTg                                                                 <  2:223517/150‑1 (MQ=255)
                                                                                          tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCNAGTTTNTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                             >  1:935522/1‑150 (MQ=255)
                                                                                               ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTg                                                        >  2:471506/1‑150 (MQ=255)
                                                                                                     gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                                  >  2:793234/1‑150 (MQ=255)
                                                                                                      agagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAaca                                                 <  1:563915/150‑1 (MQ=255)
                                                                                                       gagaATGGCGAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCAATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAacac                                                <  2:444665/150‑1 (MQ=255)
                                                                                                         gactGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGACCGCCCTGCG‑TCGGGAACGGGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                              <  2:543582/147‑1 (MQ=255)
                                                                                                         gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGACCGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                              <  2:1071334/150‑1 (MQ=255)
                                                                                                           aTGGCTATGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTc                                            <  2:916518/150‑1 (MQ=255)
                                                                                                              gCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGGCTTTTGATCAGTGCGTGATCCTGTTGAAAAACACGGTCAGc                                         >  2:475697/1‑150 (MQ=255)
                                                                                                                tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCa                                       <  1:1237007/150‑1 (MQ=255)
                                                                                                                 aaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGGCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAg                                      <  2:762407/150‑1 (MQ=255)
                                                                                                                   tgtgCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGAt                                    <  1:429313/147‑1 (MQ=255)
                                                                                                                   ggggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGAt                                    >  2:612883/1‑150 (MQ=255)
                                                                                                                    gggCAATCTTTACAAGATCCGTTCCTGAACGCCCTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                   >  2:271681/1‑150 (MQ=255)
                                                                                                                     ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                                  <  1:417239/150‑1 (MQ=255)
                                                                                                                      gCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCATTTTATATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGt                                 >  1:801658/1‑150 (MQ=255)
                                                                                                                        aaTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGttt                               <  2:61320/150‑1 (MQ=255)
                                                                                                                         aTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGTAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                              <  2:778809/150‑1 (MQ=255)
                                                                                                                         aTCTTTACAAGAACCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                              <  2:213138/150‑1 (MQ=255)
                                                                                                                           cTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                            <  1:1028565/150‑1 (MQ=255)
                                                                                                                            tttACAAGATCCGTTCCTGAACGCAGTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACaa                           <  2:548707/150‑1 (MQ=255)
                                                                                                                            tttACAAGATCCGTTCCTGAACGCACTGCG‑TCTGGAACGTGTTACAGTTTCTAATTATTTGGAGAAAGGTATTAAGCTCCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTCATCCTGTTTAAAAACACCGTCAGCCAGATGGTTTACaa                           >  2:959728/1‑150 (MQ=255)
                                                                                                                              tACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGCTCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGc                         >  2:623166/1‑150 (MQ=255)
                                                                                                                                 ggaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTCCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg                      >  1:498485/7‑150 (MQ=255)
                                                                                                                                  gaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgc                     >  1:397038/6‑150 (MQ=255)
                                                                                                                                    aTCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGCTCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                   >  1:1189452/1‑150 (MQ=255)
                                                                                                                                    aTCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                   >  1:272746/1‑150 (MQ=255)
                                                                                                                                     tCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAt                  <  1:1047532/150‑1 (MQ=255)
                                                                                                                                     cacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAt                  >  1:398397/3‑150 (MQ=255)
                                                                                                                                         ttCCTGAACGCACTGCG‑TCTGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCt              <  1:726152/150‑1 (MQ=255)
                                                                                                                                          tCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa             >  1:14146/1‑150 (MQ=255)
                                                                                                                                           ccTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAACGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc            <  1:217045/150‑1 (MQ=255)
                                                                                                                                             tGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCACCTCGAGGCTTTTGATCAGTTCGTGATCCTGTTGAAACACACCGTCAGCCAGATGGTTTACCAGCACGCGATTTCTACtg          >  2:802058/1‑150 (MQ=255)
                                                                                                                                              gAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgt         <  1:457886/150‑1 (MQ=255)
                                                                                                                                                     cTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCg  <  1:298849/150‑1 (MQ=255)
                                                                                                                                                     cTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCg  <  1:583240/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCG  >  NC_000913/4400181‑4400479

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: ‑