Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,965,670 C→T 5.7% intergenic (‑40/‑91) rhlB ← / → trxA ATP‑dependent RNA helicase/thioredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,6700CT 5.7% 199.1 / 2.5 72intergenic (‑40/‑91)rhlB/trxAATP‑dependent RNA helicase/thioredoxin 1
Reads supporting (aligned to +/- strand):  ref base C (29/38);  new base T (2/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

CCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTC  >  NC_000913/3965522‑3965814
                                                                                                                                                     |                                                                                                                                                
ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGc                                                                                                                                                  <  2:698976/149‑1 (MQ=255)
ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGc                                                                                                                                                  <  1:300774/149‑1 (MQ=255)
ccTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGc                                                                                                                                                  <  2:808663/149‑1 (MQ=255)
    aaTGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGTGTAt                                                                                                                                              >  2:389060/1‑149 (MQ=255)
    aaTGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGTGTAt                                                                                                                                              >  2:772978/1‑149 (MQ=255)
     aTGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGTGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTAt                                                                                                                                              >  2:29787/1‑149 (MQ=255)
      tGGGCTTACAGTTATGAAACCCTTTTTTTTCAAGGGGTTCTACAACCTTCTGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATcc                                                                                                                                            <  1:928710/149‑1 (MQ=255)
        ggCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGGCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCgg                                                                                                                                          <  1:88221/149‑1 (MQ=255)
         gagTACATTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTTAATGTGTTTTCGTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTAGGAAAGCGTATCCGGt                                                                                                                                         <  2:351097/147‑1 (MQ=255)
              cAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGTGCGAAGTCGGAAAACTTCTGTTCTGTTAATTGTGTTTTGCGCATAGTCTGGTAGCATATATGCTTACTAATGCTTTACAATAGGGTATCCGGtcattt                                                                                                                                    >  2:451773/1‑144 (MQ=255)
                 ttATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCTAAGTCGTAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATaaa                                                                                                                                 <  1:699994/149‑1 (MQ=255)
                 ttATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATaaa                                                                                                                                 >  2:208598/1‑149 (MQ=255)
                 ttATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATaaa                                                                                                                                 >  1:843714/1‑149 (MQ=255)
                     gAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCa                                                                                                                             >  2:141321/1‑149 (MQ=255)
                       aaCCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAAc                                                                                                                           >  1:916992/1‑149 (MQ=255)
                          ccTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCtt                                                                                                                        <  1:241133/149‑1 (MQ=255)
                           cTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTTCTCATAGTGTGGTAGAAGAGCAGCGTACTATTGCGTTACGTTAGCGTATCCGGTGATATAAAGTCCTCCttt                                                                                                                       >  2:974213/1‑149 (MQ=255)
                              tttttttCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAg                                                                                                                     <  2:791702/149‑1 (MQ=255)
                               ttttttCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGt                                                                                                                    <  2:785134/149‑1 (MQ=255)
                                  tttCAAGGGCTTCTACAACCTTCGGATGCAGCGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTgg                                                                                                                 <  1:601473/149‑1 (MQ=255)
                                   ttCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTTCTCATATTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACATATAGTTGGt                                                                                                                >  2:178419/1‑149 (MQ=255)
                                      aaGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTaa                                                                                                             >  1:75307/1‑149 (MQ=255)
                                       aGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAAt                                                                                                            >  1:505480/1‑149 (MQ=255)
                                        gggCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATg                                                                                                           <  1:208598/149‑1 (MQ=255)
                                         ggCTTCTACAACCTTCGGATGCAGGGCGTAGTCGGTAAACTTGTGTTCTGTTAAATGTGTTTTGCTCTTAGTGTTGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGt                                                                                                          <  1:741036/149‑1 (MQ=255)
                                            ttCTACAACCTTCGGATGCAGGGCGAAGTCGCAAAACTTCTTGTCTGTTAAATGTGTTTAGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAACGCGTATCCGGTGAAAGAAAGTCAACCTTTAGATAGGTAATGTTac                                                                                                       >  2:583908/1‑149 (MQ=255)
                                             gcTACAACCTGCGGATCCAGGGCGAAGTCGGAAACCTTCTGTTCTGTTAACTGTGTTTTGCTCATAGTGTGGGAGAATATCAGCTAACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGGCAACCTTTAGTTGGTTAATGTTaca                                                                                                      <  2:568575/148‑1 (MQ=255)
                                                 cAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGCAAGCGTATCCGGTGAAATAAATTCAACCTTTAGTTGGTTAAAGTTACACcaa                                                                                                  >  2:159283/1‑149 (MQ=255)
                                                 cAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCACCCTTTAGTTGGTTAATGTTACACcac                                                                                                  >  1:602650/1‑148 (MQ=255)
                                                 cAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGATATGCTCATAGTGTGTTAGAATATCAGCATACTATTTCTTTACGAAATCGTATCCGATGAAATAAATTCAACATTAAGTTGATTAATTTTACACcaa                                                                                                  >  1:605459/1‑149 (MQ=255)
                                                    ccTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACcaacaa                                                                                               >  1:470519/1‑149 (MQ=255)
                                                      ttCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCAGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACAAACAACg                                                                                             >  1:644837/1‑149 (MQ=255)
                                                      atcttaTGCAGCGATAAGTCGTAAAAATTCTTTTCTGTTAAATGAGTATTGATCATATTGTGGTAGAAAATCAGCTAACTATTGCTTTACGAAAGCGAATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACg                                                                                             <  1:207514/144‑1 (MQ=255)
                                                        cGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATCGTGTGGTAGAATCTCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGaa                                                                                           <  2:206814/149‑1 (MQ=255)
                                                          gATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAAc                                                                                         <  1:808612/149‑1 (MQ=255)
                                                          gATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGATAATGTTACACCAACAACGAAAc                                                                                         >  2:281887/1‑149 (MQ=255)
                                                             gCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGACACCaa                                                                                      >  1:778293/1‑149 (MQ=255)
                                                                gggCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAacac                                                                                   <  1:654450/149‑1 (MQ=255)
                                                                       gTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGc                                                                            <  2:935983/149‑1 (MQ=255)
                                                                          ggAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTa                                                                         <  1:622511/149‑1 (MQ=255)
                                                                           gAAAACTTCTGTACTGTTAAATGTGTTTTGATCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATATTACACCAACAACGAAACCAACACTCCACGATTAt                                                                        >  1:583421/1‑149 (MQ=255)
                                                                             aaaCTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTc                                                                      <  2:251336/149‑1 (MQ=255)
                                                                                    tgttctgttAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGa                                                               <  1:706462/149‑1 (MQ=255)
                                                                                       tctgttAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGtt                                                            <  2:475120/149‑1 (MQ=255)
                                                                                           ttAAATGTGTTTTTCTCATAGTGTTGTATAATATAATCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCACCACGCCAGGCTTATTCCTGTGGAGTtatat                                                        <  2:602650/149‑1 (MQ=255)
                                                                                           ttAAATGTGTTTTGCTCTTTGTGTGGTAGAATATCAGCATACTATTGGTTTACGAAAGCGTATCCTGTGAAATAAAGTCAACCTTTAGTTGGTTAATTTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTtatat                                                        <  2:605459/149‑1 (MQ=255)
                                                                                                tgttTTTTGCTCATAGTGTTCTAGATTATATGCTTACTATTGCTTTAGTAAAGCGTATCCGGTGAAATACAGTCAACCTTTAGTTGGTTTATGTTACACCAACAACGAACCCAACACGCCAGGCTTATTCCTGTGGAGGTATATATGAg                                                   <  2:69865/149‑1 (MQ=255)
                                                                                                    ttttGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGAt                                               <  1:983879/149‑1 (MQ=255)
                                                                                                      ttGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAACGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATaa                                             <  1:989605/149‑1 (MQ=255)
                                                                                                        gCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATaaaa                                           <  2:843714/149‑1 (MQ=255)
                                                                                                         cTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAat                                          >  2:711092/1‑149 (MQ=255)
                                                                                                           cATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAatta                                        <  2:931270/149‑1 (MQ=255)
                                                                                                            aTAGTGTGGTAGAATATCAGCTTACTAATGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAattat                                       >  1:35487/1‑149 (MQ=255)
                                                                                                                 gtgGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAAAAACGAAACCAACACGCCAGGCTAATTCCTGTGGAGTTATATATGAGCGAAAAAATAAATCAca                                  >  2:503758/1‑148 (MQ=255)
                                                                                                                  tgGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTTTTCACct                                 >  2:100459/1‑149 (MQ=255)
                                                                                                                    gTAGAATATCAGCTTACTATTGCTTTACGAAATCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACctga                               <  2:651386/149‑1 (MQ=255)
                                                                                                                     tAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCAActgac                              >  1:615980/1‑149 (MQ=255)
                                                                                                                        aaTATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACctgactga                           <  1:386114/149‑1 (MQ=255)
                                                                                                                        aaTATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACctgactga                           >  1:627827/1‑149 (MQ=255)
                                                                                                                         atatCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCtgactgac                          <  2:189936/149‑1 (MQ=255)
                                                                                                                           atCAGCTTACTATTGCTTTACGAAATCGTATCCGGTGAAATAAAGTCAAACTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTgacga                        <  1:431713/149‑1 (MQ=255)
                                                                                                                               gCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACATTTTGTTGGTTTATGTTACACCAACAACGAAACCATCACGCCAGGCTTATTCATGTGGAGTTATATATGAGCGATAATATTATTCTCCTGACTGAGGACAGt                    <  2:97031/149‑1 (MQ=255)
                                                                                                                                  tACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGTAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGtttt                 <  2:916992/149‑1 (MQ=255)
                                                                                                                                  tACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGATGGATAATGTTACACCAACAACGAAACCAACACGCCACGCTTAATCCTGAGGAGTTATATAAGAGCGATAAAATTATTCACCTGACTGACGACAGtttt                 <  2:580248/149‑1 (MQ=255)
                                                                                                                                    atATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGa               <  1:1015775/148‑1 (MQ=255)
                                                                                                                                     tATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGGTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGac              <  1:866947/149‑1 (MQ=255)
                                                                                                                                        tGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACg           <  1:974213/149‑1 (MQ=255)
                                                                                                                                          cTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCCACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGa         >  1:474377/1‑149 (MQ=255)
                                                                                                                                           tttACGAAAGTGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGAt        <  1:389060/149‑1 (MQ=255)
                                                                                                                                           tttACGAAAGTGTATCCGGTGAAATAAAGGCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGAt        <  1:772978/149‑1 (MQ=255)
                                                                                                                                           tttACGAAAGCGTATCCGGTGAAAAAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGAt        >  2:949057/1‑149 (MQ=255)
                                                                                                                                              aCGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTa     >  2:433112/1‑149 (MQ=255)
                                                                                                                                               cGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTAc    <  1:541086/149‑1 (MQ=255)
                                                                                                                                                gAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACAGACGACAGTTTTGACACCGATGTACt   >  2:354269/1‑149 (MQ=255)
                                                                                                                                                 aaaGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTc  >  1:15292/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTC  >  NC_000913/3965522‑3965814

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: