Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,344,036 G→C 7.3% T104S (ACT→AGT)  yciT ← global regulator of transcription; DeoR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,344,0360GC 7.3% 250.8 / 10.3 93T104S (ACT→AGT) yciTglobal regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base C (5/2);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAAC  >  NC_000913/1343892‑1344184
                                                                                                                                                |                                                                                                                                                    
cGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCg                                                                                                                                                  <  2:1004587/149‑1 (MQ=255)
  tgtCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCggg                                                                                                                                                <  1:1382082/149‑1 (MQ=255)
   gtCTAAGTGCCAACCATAGTATCGCTTTTTTTCTGTTACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGc                                                                                                                                               <  2:917317/149‑1 (MQ=255)
   gtCAAAGGGCCAACCATAGTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGAGGTAGCTGCTGCCCGTGATGATAGTGCCATTTTTCTTCTGCTCGCCCAGAGTCCGGGc                                                                                                                                               >  2:796526/1‑149 (MQ=255)
     cAAATGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGAGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                                                                             <  1:402915/149‑1 (MQ=255)
     cAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTGACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                                                                             <  1:973362/149‑1 (MQ=255)
     cAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGCCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGATTCCGGGCGa                                                                                                                                             >  1:1021574/1‑149 (MQ=255)
     cAAAGGGCCAACCATACATTCGCTTTTTATCTGATACACGCCACCGAGCAAAATAACTTCACAACGCGCGTCTTTCAGCAAAAGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGa                                                                                                                                             <  2:696155/149‑1 (MQ=255)
      aaaGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGAAATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAg                                                                                                                                            >  2:1476224/1‑149 (MQ=255)
         gggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAGCAg                                                                                                                                         <  2:1270320/149‑1 (MQ=255)
          ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                        >  2:976301/1‑149 (MQ=255)
          ggCAAACCATACTTTCGCTTTTTTTCTGATACACGCCACGGAGCACAATAACTTCACAAGGCGCGTCTTTCAGAAAATGCGCTATGTAGCTGCTGACCTTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                        <  2:291960/149‑1 (MQ=255)
           gCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCAAGAGTCCGGGCGAGTAggg                                                                                                                                       >  1:856541/1‑149 (MQ=255)
            ccAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCCCAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGATCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                                                                      >  1:165994/1‑149 (MQ=255)
            ccAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                                                                      >  1:438693/1‑149 (MQ=255)
            ccAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                                                                      >  1:155913/1‑149 (MQ=255)
             cAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGCGCa                                                                                                                                     >  2:239948/1‑149 (MQ=255)
                  aTACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGAGATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtggt                                                                                                                                >  1:721137/1‑147 (MQ=255)
                  aTACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCGTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgat                                                                                                                                >  2:1426850/1‑147 (MQ=255)
                  aTACTTTCGATTTTTTTCTGATACACGCCCCCGAGCAAAATAACTTCCCAAGGCGCGTCTTTCAGCAAATGCGCGAGGTAGCTGCTGACCGTGATGCTAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgct                                                                                                                                <  2:405842/149‑1 (MQ=255)
                     atttCGCTTCTTTTCTGATACACTCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGAATTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgct                                                                                                                             <  2:134640/148‑1 (MQ=255)
                      tttCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                                                                            >  2:7330/1‑149 (MQ=255)
                          gCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAt                                                                                                                        >  1:191410/1‑149 (MQ=255)
                           cTTTTTTTCTGATACACGCCACCGAGCAAATTACCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtt                                                                                                                       >  2:1295661/1‑149 (MQ=255)
                            tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                                                                      >  2:1265809/1‑149 (MQ=255)
                            tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCCTGATGATAGTGACATGTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                                                                      >  2:722136/1‑149 (MQ=255)
                             ttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtttt                                                                                                                     >  2:368336/1‑149 (MQ=255)
                              tttttCTGATACACGCCACCGAGCAAAATACCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATTATCGTGAGATTTTTCTTCTGCTCGCCGAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                                                                    <  2:89238/149‑1 (MQ=255)
                              tttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                                                                    <  2:405381/149‑1 (MQ=255)
                                   gtGATACTCGCCACGGAGCAAATTAACTTGACATGGCGCGTCTTTCAGCAAATGCGCGATGTTGCTGCTGAGCGTGATGATTGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaa                                                                                                               <  2:321061/148‑1 (MQ=255)
                                    tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaaa                                                                                                              >  2:1384763/1‑149 (MQ=255)
                                    tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCAGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaaa                                                                                                              >  1:58597/1‑149 (MQ=255)
                                     gATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAg                                                                                                             <  1:765351/149‑1 (MQ=255)
                                        acacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATg                                                                                                          <  1:471583/149‑1 (MQ=255)
                                          acGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                        >  2:1531758/1‑149 (MQ=255)
                                           cGCCACCGAGCAAAATAACTGCAGAAGGCGCGTCTTTCAGCAAATGCGCGAGGTAGCTGCGGCCCGGGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATATAGATGGtt                                                                                                       <  2:1382709/149‑1 (MQ=255)
                                            gCCACCGAGCAAAATAACTTCTCAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGttt                                                                                                      <  2:1482014/149‑1 (MQ=255)
                                             ccACCGATCATAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                     <  2:165994/149‑1 (MQ=255)
                                             ccACCGAGCAAAATAAGTTCACACTGCGCGTCTTTCAGCAAATGCGCGATGTAGCTTCGGACCGTGAGGATAGTGACATTTTTCTACTGCTCGCCCAGAGTCCCTGCGAGTAGGGCATTGCTGCTTCCATTTTCGATAAAGATGGTTTc                                                                                                     <  2:1037618/149‑1 (MQ=255)
                                             ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                     <  1:1292967/149‑1 (MQ=255)
                                                ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGcc                                                                                                  >  1:359556/1‑149 (MQ=255)
                                                  gAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCGAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCg                                                                                                <  1:580693/149‑1 (MQ=255)
                                                  gAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCg                                                                                                <  2:1471562/149‑1 (MQ=255)
                                                   aGCAAAATACCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCgg                                                                                               >  1:107128/1‑149 (MQ=255)
                                                    gCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCTCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGt                                                                                              <  2:1308207/149‑1 (MQ=255)
                                                     ctaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGtt                                                                                             <  2:568125/147‑1 (MQ=255)
                                                        aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTCCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGCTTGaa                                                                                          >  1:697077/1‑149 (MQ=255)
                                                          tAATTTCACATGGCGCGTTTTTCTGCAAATGCGCGATGTAGCTGCTGCGGGTGATGATAGTGACTTTTTTCTTCTTATCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCAATTTTCGATAAAGATGGTTTCGCCCGGTTGAAcc                                                                                        <  1:1528457/149‑1 (MQ=255)
                                                                 acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAcg                                                                                 >  2:1116420/1‑149 (MQ=255)
                                                                     ggCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGc                                                                             <  2:15193/149‑1 (MQ=255)
                                                                      gcgcGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCa                                                                            >  1:396050/1‑149 (MQ=255)
                                                                            cTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCg                                                                      >  1:961202/1‑149 (MQ=255)
                                                                            cTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGACATTTTCGATAAACATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCg                                                                      >  1:357675/1‑149 (MQ=255)
                                                                              ttCAGCAAATGCGCGATGTAGCTGCTGAGCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGc                                                                    <  2:1523923/149‑1 (MQ=255)
                                                                              ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGGGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGc                                                                    <  1:251598/149‑1 (MQ=255)
                                                                                cAGCTACTGCGCCATGTAGCTGCTGACCGTGATGATCGGGACCATTTTCTGCTGCGCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCTTTTTCGATAAAGATGGTTTCGCCCGGTTGAACGAGCGACGCGGCAAACTCGGCGa                                                                  <  2:864395/149‑1 (MQ=255)
                                                                                 aGCAAATGCGCTATGTAGCGTCTGAGCGTGATGATAGTGACATTTTTCTTCTGCTCTCCCAGATTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTGCGCCCGGTTTAACCAGCGACGCGGCAAACTCGGCGAg                                                                 <  2:1158459/149‑1 (MQ=255)
                                                                                          gcgATGTAGCTCCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCGACCATTTTCGATAAAGATGGTTACGACCGCTTGCACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                        >  2:1208673/1‑149 (MQ=255)
                                                                                          gcgATGAAGCTGCTGACCGTGATGATAGAGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGATTCGCCCGGTTGAACCAGCGACGCGGCAAACTCTGCGAGTTCACTCtt                                                        >  2:1549739/1‑149 (MQ=255)
                                                                                              tGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGc                                                    <  2:191410/149‑1 (MQ=255)
                                                                                                tAGCTGCTGCCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGCGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAAATCGGCGAGTTCACGCTTCAGcct                                                  >  2:779884/1‑147 (MQ=255)
                                                                                                tAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGt                                                  >  1:250101/1‑149 (MQ=255)
                                                                                                 aGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGta                                                 >  2:1568066/1‑149 (MQ=255)
                                                                                                    tgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAg                                              <  1:368336/149‑1 (MQ=255)
                                                                                                      ctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGGTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGAGAGTTCACGCTTCAGCGGATAGtt                                            >  2:107925/1‑149 (MQ=255)
                                                                                                             tgatgaTAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtcatc                                     <  2:961202/149‑1 (MQ=255)
                                                                                                              gatgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtcatca                                    <  1:1384763/149‑1 (MQ=255)
                                                                                                               atgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGATTCGCCCGGTTGAACCAGCGACACGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtcataaa                                   >  2:1042045/1‑146 (MQ=255)
                                                                                                                 caaaGTGACATTATTCTTCTGCTCCCCCCCAGTCCGGGCGAGTAGGGCACTGCTGCAGCCATTATCGATAAAGATGGATTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATAc                                 <  1:796526/146‑1 (MQ=255)
                                                                                                                   tAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACgg                               <  1:7330/149‑1 (MQ=255)
                                                                                                                    aGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACggg                              >  1:1116330/1‑149 (MQ=255)
                                                                                                                    aGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACggg                              >  1:1111591/1‑149 (MQ=255)
                                                                                                                     gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACg                                >  2:872129/1‑146 (MQ=255)
                                                                                                                     gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCGGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACg                                <  1:872129/146‑1 (MQ=255)
                                                                                                                      tGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGtt                            >  2:14978/1‑149 (MQ=255)
                                                                                                                          aTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCgc                        >  1:653853/1‑148 (MQ=255)
                                                                                                                          aTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGa                        <  1:590621/149‑1 (MQ=255)
                                                                                                                               tgccctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATGTTCGATAAAGATGGTTTCGCCCGGTTGAACCGGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtca                   >  1:1204157/8‑149 (MQ=255)
                                                                                                                               tcttctGCTCGCCCAGAGTCCGCGCGAGTAGGGCATTGCTGCGGCCATTTTCGATAAAGATGGGTTCGCCCGGTTGAACCAGCGACGCGGCATCCTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGCGTTTCGACGtca                   <  1:74672/149‑1 (MQ=255)
                                                                                                                                cttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGATCCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcat                  <  2:58597/149‑1 (MQ=255)
                                                                                                                                 ccctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACAAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatc                 >  1:1229690/6‑149 (MQ=255)
                                                                                                                                   ctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCAc               >  1:283084/1‑149 (MQ=255)
                                                                                                                                    tacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACt              >  2:892510/3‑149 (MQ=255)
                                                                                                                                     gCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTa             >  1:1544906/1‑149 (MQ=255)
                                                                                                                                     gCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTa             <  1:976301/149‑1 (MQ=255)
                                                                                                                                     gCTCGCCAAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTa             <  2:856541/149‑1 (MQ=255)
                                                                                                                                      cTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTAt            >  1:468077/1‑149 (MQ=255)
                                                                                                                                       tCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATCGTTTCGACCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTTTTGTGGCTCATCAGACGGGTTTCGACGTCATCACTATc           >  1:940476/1‑149 (MQ=255)
                                                                                                                                       tCGCCAAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCTATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCAGACGTGTATCGACGTCATAACGATc           >  1:1510900/1‑149 (MQ=255)
                                                                                                                                         gCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTGTCGATAAAGATGGTTTCGCCCGGTTGAACCATCGACGCGGCAAACTCGGCGAGTTCACGCTTAAGCGTATAGTTGATCATCATACGGGTTTCGACGTAATCACTATCaa         >  1:1462023/1‑149 (MQ=255)
                                                                                                                                          cccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTAACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAg        >  1:913560/1‑149 (MQ=255)
                                                                                                                                            cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTTTAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCg      <  1:1265809/149‑1 (MQ=255)
                                                                                                                                            cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCg      >  1:1282444/1‑149 (MQ=255)
                                                                                                                                                gTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTGCGATAAAGATGGTTTCGCCCGGTGGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAAc  <  2:504790/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAAC  >  NC_000913/1343892‑1344184

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