Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,440,923 C→A 18.7% G36V (GGC→GTC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,9230CA18.7% 127.9 / 24.7 73G36V (GGC→GTC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (36/23);  new base A (8/6);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

CAGGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCT  >  NC_000913/3440775‑3441071
                                                                                                                                                    |                                                                                                                                                    
cAGGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGc                                                                                                                                                      <  1:630527/149‑1 (MQ=255)
  ggTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAAACTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGTCGCAGTGCGTTACCCAGAGTATGGCCa                                                                                                                                                    >  2:1009992/1‑149 (MQ=255)
      gAGCAGGATTTCCAGGATATCTTCCTGCACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAACGc                                                                                                                                                >  1:811907/1‑149 (MQ=255)
          aGGATTTCCAGGATATCATCCAGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCACTCTCAACCTCGGTCACCGCGCAACCCGGCACCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACg                                                                                                                                            <  1:53588/149‑1 (MQ=255)
            gATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGct                                                                                                                                          >  2:336337/1‑149 (MQ=255)
            gATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGACACCGCGCAACCCGGCATCGAAGAGAGCAGAATACGGCGCAGTCCGTTACCCAGAGTATGGACAAAGCCACGca                                                                                                                                          >  1:886823/1‑148 (MQ=255)
                 ccAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTaaa                                                                                                                                     >  1:936968/1‑149 (MQ=255)
                    ggATATCTTCCTGAACGACTTCTTTTGTGCTGTACTCATGTATTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTAACCCATAGTATGGCCAAAGCCACGCTATAACGGc                                                                                                                                  >  1:671403/1‑149 (MQ=255)
                     gATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCt                                                                                                                                 >  2:890354/1‑149 (MQ=255)
                          cTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAgg                                                                                                                            <  1:912002/149‑1 (MQ=255)
                             ccTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCACCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCTGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                         >  2:585659/1‑149 (MQ=255)
                             ccTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCTATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                         <  1:917546/149‑1 (MQ=255)
                                 aaCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCt                                                                                                                     >  2:839949/1‑149 (MQ=255)
                                   cGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGAATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCTGAGTATGGACAAAGCCCCGCTCTTAAGGCTCAAGGGTCACCGTg                                                                                                                   >  1:601675/1‑149 (MQ=255)
                                    gCCTTCTTTGGTGCTGTACTCTTGTAGTACACAATCAATCTCAACCTCGGTCACCGCGCAACCCGTCATCGATGAGAGCAGAATACTGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTgg                                                                                                                  <  2:753129/149‑1 (MQ=255)
                                     ccTTCTTTGGTGCTGTACTGATGTAGTACACGATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGGTACCCAGAGTATGGCCAAAGCCACGCTGTAAAGGCTCAAGGGTCACCTTGgc                                                                                                                 <  2:439041/149‑1 (MQ=255)
                                       ttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGCTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgt                                                                                                               <  2:120438/149‑1 (MQ=255)
                                         cttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGACTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgagc                                                                                                             >  2:103921/1‑149 (MQ=255)
                                                      cTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCGCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                                <  2:295703/149‑1 (MQ=255)
                                                       tCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTg                                                                                               >  2:407237/1‑149 (MQ=255)
                                                         aTGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                             >  1:552431/1‑149 (MQ=255)
                                                         aTGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGACAGCAGAATACGGCGCAGTGCGTTACCCTGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                             <  1:446978/149‑1 (MQ=255)
                                                          tGTAGTAGACCATCTATCTCAACCTCGGGCACCGCGCAACCCGGGAGCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGGATGGCCACAGCCACGCGCTAAAGGCTCAAGGGTCACCTTGGCGGGCGTCGAACTCACTTGCTc                                                                                            <  2:671403/149‑1 (MQ=255)
                                                           gtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCCCCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                           >  1:133276/1‑149 (MQ=255)
                                                               taCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACATTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                                       >  1:929541/1‑149 (MQ=255)
                                                                 cacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGAGCGTTACCCACAGAAAGCCCACAGCCACGCTCTAAAGGCTCAAGGGTCACCACGGCGTGCGTCGAACGACATTGCTCGAAATCa                                                                                     >  1:184837/1‑149 (MQ=255)
                                                                  acCATCAATCTCAACCTCTGTCACCGCGCAACCCTGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCaa                                                                                    <  1:354498/149‑1 (MQ=255)
                                                                   ccttCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                                   <  1:169742/146‑1 (MQ=255)
                                                                    cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                                  >  2:520168/1‑149 (MQ=255)
                                                                    cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                                  <  1:640195/149‑1 (MQ=255)
                                                                     aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa                                                                                 >  2:160242/1‑149 (MQ=255)
                                                                           ctcAACCTCGGTCACCGCGCACCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAACGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                           >  2:890440/1‑149 (MQ=255)
                                                                               aCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                                       <  1:861965/149‑1 (MQ=255)
                                                                               aCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                                       <  1:336337/149‑1 (MQ=255)
                                                                                     gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTGTTAGaaa                                                                 >  1:54249/1‑149 (MQ=255)
                                                                                      tCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGGCGAAATCACTTGCTCGATATCCACACGGCGCGGTTATAGCAAc                                                                >  1:562907/1‑149 (MQ=255)
                                                                                         ccGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActct                                                             >  2:70957/1‑149 (MQ=255)
                                                                                         ccGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActct                                                             <  2:115381/149‑1 (MQ=255)
                                                                                          cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTg                                                            >  1:138093/1‑149 (MQ=255)
                                                                                               aaCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                       >  2:537436/1‑149 (MQ=255)
                                                                                                aCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGACACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCCCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACGCTGTCACAg                                                      >  2:707584/1‑149 (MQ=255)
                                                                                                 cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTAACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAgc                                                     >  1:871668/1‑148 (MQ=255)
                                                                                                      cATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGTCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCt                                                <  2:645609/149‑1 (MQ=255)
                                                                                                          gATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGCACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAt                                            >  2:76849/1‑149 (MQ=255)
                                                                                                           aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGCACCCTGCAtt                                           >  1:849353/1‑149 (MQ=255)
                                                                                                           aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGCACCCTGCAtt                                           >  2:311477/1‑149 (MQ=255)
                                                                                                              agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgtg                                        <  2:936968/149‑1 (MQ=255)
                                                                                                               gagCAGAATACGGCGCAGTGCGTTACCCAGAGTTTGGCCAAAGACACGCTCTAAAGGCTCATGGGTCACCTTGGCTTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTGACAGAACCCTGCATTgtgt                                       <  2:845690/149‑1 (MQ=255)
                                                                                                               gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCCAGGGTCACCTTGGCGTGGGTCGACCTCACTTGCTCGATATCAACCAGGCGCGGTGTTAGCAACTCTGTCCCAGACCCCTGCATTgggt                                       >  1:920857/1‑149 (MQ=255)
                                                                                                                agCAGAATACGGCGCAGGGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                      <  1:307085/149‑1 (MQ=255)
                                                                                                                  cAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGGACCCTGCATTGTGTcca                                    >  2:675303/1‑148 (MQ=255)
                                                                                                                   aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctc                                   >  2:932275/1‑149 (MQ=255)
                                                                                                                    gAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGCACTCACTTCCTCGATATCAACCAGGCGCGGTTGTAGAAACGCTGTCACAGAACCCTGGATTGGGTGctct                                  >  1:185388/1‑149 (MQ=255)
                                                                                                                     aaTACGGCGCAGTGCGGTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                                 <  2:407068/149‑1 (MQ=255)
                                                                                                                       tACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGGCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                               <  2:418397/149‑1 (MQ=255)
                                                                                                                        aCGGCGAAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAACCTGCATTGTGTCCTCTCttt                              >  1:12664/1‑149 (MQ=255)
                                                                                                                          ggCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGATCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTgg                            <  1:1009992/149‑1 (MQ=255)
                                                                                                                           gcgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGt                           >  2:960488/1‑149 (MQ=255)
                                                                                                                            cgcAGTGCGTTATCCAGAGTATGGACAAATCCACGCTCTAATGTCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCATGCATTGTGTCCTCTCTTTGGTa                          <  2:367827/149‑1 (MQ=255)
                                                                                                                             gcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTAc                         >  2:427569/1‑149 (MQ=255)
                                                                                                                               aGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAACGCTCAAGGGTCCCCTTGGCTTGCGTCGAACTCACTTGCTCGATATCAACCATGAGCCGTTTTAGAAACTCTTTCACACAACCCTGCATTGTGTCCTCTCTTTGctacta                       >  1:535096/1‑149 (MQ=255)
                                                                                                                                gTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTaa                      <  2:152958/149‑1 (MQ=255)
                                                                                                                                 tGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAg                     >  2:356770/1‑149 (MQ=255)
                                                                                                                                  gCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGc                    <  1:890354/149‑1 (MQ=255)
                                                                                                                                    gTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCGCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCtt                  <  1:585659/149‑1 (MQ=255)
                                                                                                                                     ttACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCCCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt                 >  2:553117/1‑149 (MQ=255)
                                                                                                                                       aCCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGCTATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATGGTGTCCTCTCTTTGGTACTAAGCTTTAc               >  1:440678/1‑149 (MQ=255)
                                                                                                                                         ccAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt             >  2:369171/1‑149 (MQ=255)
                                                                                                                                         ccAGAGTATGGACCAAGCCACGCGCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAATGTGTCCTCTCTTTGGTACTAAGCTTTACtt             <  2:572122/149‑1 (MQ=255)
                                                                                                                                          cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg            <  1:840739/149‑1 (MQ=255)
                                                                                                                                               tATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTGCTCTCTTTTGTACTAAGCTTTACTTGGAGTa       >  2:680913/1‑149 (MQ=255)
                                                                                                                                                  ggCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAg    <  1:839949/149‑1 (MQ=255)
                                                                                                                                                   gCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc   >  2:536356/1‑149 (MQ=255)
                                                                                                                                                   gCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc   >  1:319280/1‑149 (MQ=255)
                                                                                                                                                    ccAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt  >  1:224757/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCT  >  NC_000913/3440775‑3441071

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