Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,597 C→T P750L (CCG→CTG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5970CT97.8% 292.6 / ‑4.1 92P750L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (48/41);  minor base G (0/2);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTC  >  NC_000913/4187451‑4187740
                                                                                                                                                  |                                                                                                                                               
gCAAACTGAAACCGTGATTAACCGTCACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTgg                                                                                                                                               >  1:811385/1‑149 (MQ=255)
   aaCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATg                                                                                                                                            <  1:188994/149‑1 (MQ=255)
    aCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATgg                                                                                                                                           <  2:800517/149‑1 (MQ=255)
         aaCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCGGGATGGCTCca                                                                                                                                      <  1:557292/149‑1 (MQ=255)
          aCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcat                                                                                                                                     <  1:1028706/149‑1 (MQ=255)
          aCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcat                                                                                                                                     >  2:513251/1‑149 (MQ=255)
           ccGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcatc                                                                                                                                    >  1:857692/1‑149 (MQ=255)
             gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcatcat                                                                                                                                  <  2:134618/149‑1 (MQ=255)
              tGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCatcatc                                                                                                                                 <  2:1280993/149‑1 (MQ=255)
                aTTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGa                                                                                                                               >  1:1300635/1‑149 (MQ=255)
                     ccGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCACCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCCCAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGGTGGCTCCATCATCGAAACGc                                                                                                                          >  2:869050/1‑149 (MQ=255)
                       gTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCa                                                                                                                        >  1:845745/1‑149 (MQ=255)
                       gTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCa                                                                                                                        <  1:166659/149‑1 (MQ=255)
                              tCAGGAAGAGAAGCAGGTTTCCTTCAACAGCACCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACcg                                                                                                                 <  2:1048965/149‑1 (MQ=255)
                                aGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGTGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACcgcg                                                                                                               <  2:731992/149‑1 (MQ=255)
                                aGGAAGAGAAGCAGGTTTCCTTCAACAACATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCATGGTCTGATGGCGAAGCTGGATGGCTCCATCACCGAAACACCACTCCCcccg                                                                                                               >  2:180157/1‑149 (MQ=255)
                                  gAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGTGCTGATGGCGAAGCTGGATGGCTCCATCATCGACACGCCAATCACCGCGaa                                                                                                             >  1:451328/1‑149 (MQ=255)
                                  gAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGaa                                                                                                             >  1:1094952/1‑149 (MQ=255)
                                    agagAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGGACt                                                                                                           >  2:180359/1‑149 (MQ=255)
                                    aaccacgCAAGTTTCCTTCAACAGCATCTACATGATCGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAACCGGCAATCGCCGCCAAcc                                                                                                           >  2:904162/7‑148 (MQ=255)
                                      agaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTc                                                                                                         <  2:837608/149‑1 (MQ=255)
                                             ggTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCTAACTTCCGTGAAg                                                                                                  <  1:926335/149‑1 (MQ=255)
                                              gTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                  <  1:1327682/148‑1 (MQ=255)
                                              gTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                  >  2:1327682/1‑148 (MQ=255)
                                                 tCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAGGGTCt                                                                                              >  1:1096153/1‑149 (MQ=255)
                                                 tCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGACGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCt                                                                                              <  1:388606/149‑1 (MQ=255)
                                                  ccTTCAACAGCATCTACAGGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATCCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTg                                                                                             <  1:470003/149‑1 (MQ=255)
                                                     tCAACCCCATCTACGTGCTGGCCGCCCCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCtttgg                                                                                          <  2:508660/149‑5 (MQ=255)
                                                     tCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                                                                          <  1:492801/149‑1 (MQ=255)
                                                        aCAGCCTCTACAGGATGGCCGTCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTa                                                                                       <  2:376677/149‑1 (MQ=255)
                                                         cAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGACACGCCAATCACCGGGAACTTCGGTGAAGGTCTGAACGTAc                                                                                      >  1:707727/1‑149 (MQ=255)
                                                         cAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTAc                                                                                      >  1:1044257/1‑149 (MQ=255)
                                                           gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                                                                    <  2:731159/149‑1 (MQ=255)
                                                           gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGCCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                                                                    <  1:823118/149‑1 (MQ=255)
                                                                 aCATGATGGCCGACTCCGGTTCGCGTGGTTCTGCGGCACAGATTCGTCAACTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCTAACTTCCGTGAAGGCGGGTACGTACGCCAGTag                                                                              >  2:26508/1‑148 (MQ=255)
                                                                 aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGTACTTCCGGGAAGGGCTGAACGTACTCCAGTAc                                                                              >  2:762876/1‑149 (MQ=255)
                                                                  cATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACt                                                                             >  2:1332293/1‑149 (MQ=255)
                                                                    tgatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCAGAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTc                                                                           >  2:819312/1‑149 (MQ=255)
                                                                       tgGCCGACTCCGGTGCGCGTGGTTCTGCGGACCATTTGCTGTAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtc                                                                        >  1:665018/1‑149 (MQ=255)
                                                                          ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                                     >  2:529803/1‑149 (MQ=255)
                                                                          ccGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCAGCATCGGAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGTACGTACTCCAGTACTTCATcccc                                                                     >  1:524224/1‑149 (MQ=255)
                                                                            gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccac                                                                   <  2:386429/149‑1 (MQ=255)
                                                                            gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccac                                                                   <  2:1094952/149‑1 (MQ=255)
                                                                            gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccac                                                                   <  1:867544/149‑1 (MQ=255)
                                                                            gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTATGAACGTACTCCAGTACTTGATCTccac                                                                   >  1:367443/1‑149 (MQ=255)
                                                                              cgccGGTGCGCGGGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccc                                                                 <  1:1238613/147‑1 (MQ=255)
                                                                               tCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccca                                                                >  1:115137/1‑149 (MQ=255)
                                                                               tCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccca                                                                <  2:665018/149‑1 (MQ=255)
                                                                                 cGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCCCCGCGAACTTCCGTGGAGGTCTGAACGTAATCCAGTACTTCATCTccacgcccg                                                              >  2:341459/1‑144 (MQ=255)
                                                                                 cGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACg                                                              >  1:1032426/1‑149 (MQ=255)
                                                                                       gcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGAACTTCATCTCCACCCACGGTGCTc                                                        <  1:1008826/149‑1 (MQ=255)
                                                                                           ggTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTaa                                                    >  1:539581/1‑149 (MQ=255)
                                                                                            gTTCTGCGGCACAGATTCGGCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTaaa                                                   <  1:177666/149‑1 (MQ=255)
                                                                                             ttCTGCGGCACAGATTCGTCACATTGCTGGTATGCGCTGTCTGATGGCGAAGCTCGATGGCTCCATCATCCAAACGCCAATCACCGCGAACTTCAGTGAAGGTCTAACCGGCCTCCAGTACTGCAACTCCACCCACCGTGCtggaagag                                                  >  2:127206/1‑142 (MQ=255)
                                                                                                   ggCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTgg                                            >  1:614592/1‑149 (MQ=255)
                                                                                                   ggCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTgg                                            >  1:614718/1‑149 (MQ=255)
                                                                                                    gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCGTCTCCCCCCACGGTGCTCGTAAAGGGCTGGc                                           >  2:768506/1‑149 (MQ=255)
                                                                                                    gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAgg                                                 >  1:145442/1‑143 (MQ=255)
                                                                                                    gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAgg                                                 <  2:145442/143‑1 (MQ=255)
                                                                                                    gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                           <  1:364180/149‑1 (MQ=255)
                                                                                                    gCACAGATTCGTCAGCTTGCTGATATGCGTGGTCTGCTGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATGACCGCGAACTTCGGTGAAGGTCTGAACGTACTCCAGTTCTTCATCTCCACCCACGGGGCTCGTAAAGGTCTggt                                           >  2:81757/1‑148 (MQ=255)
                                                                                                         gATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCACGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCGTCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATa                                      >  1:656000/1‑149 (MQ=255)
                                                                                                         gATTCGTCAGCTTGCGGTAATGCGGGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATa                                      <  2:968784/149‑1 (MQ=255)
                                                                                                          aTTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATAc                                     >  1:450364/1‑149 (MQ=255)
                                                                                                             cgccAGCTGGCTGGCATGCGTCGTCCGATCGCGAAGCTGGATCGCCACTACCCGCAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGc                                  <  1:38948/146‑1 (MQ=255)
                                                                                                              ggcAGCTGGCGGGTAGGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCGGCGGGTTTCGTAAAGGTCTGGCGGATACCGCa                                 <  2:967614/147‑1 (MQ=255)
                                                                                                              gTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCa                                 <  1:513251/149‑1 (MQ=255)
                                                                                                              gTCAGCTTGCTGGTAGGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCa                                 <  2:1241592/149‑1 (MQ=255)
                                                                                                               tCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCAc                                >  1:577972/1‑149 (MQ=255)
                                                                                                               tCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGCAAAGGGCTGGCGGATACCGCAc                                >  2:1176274/1‑149 (MQ=255)
                                                                                                               ccAGCTCGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCAc                                <  1:681470/148‑1 (MQ=255)
                                                                                                                 aGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGTCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTg                              >  1:292339/1‑149 (MQ=255)
                                                                                                                 aGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCCGTACTTCGTCTCCCCCCACGGTCCTCGTAAAGGTCTGGCGGATACCGCCCTg                              >  2:128996/1‑149 (MQ=255)
                                                                                                                      gCTGGTATGCTTGGTCATATGGCAGAGCTGGATGGCTCACTCAATGCATAACTAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAAc                         <  1:1323309/149‑1 (MQ=255)
                                                                                                                       cTGGTATGCGTCGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACt                        <  1:1289242/149‑1 (MQ=255)
                                                                                                                        tGGTATGCGTGGTCTGGTGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTg                       <  2:383741/149‑1 (MQ=255)
                                                                                                                        tGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCCCGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTg                       >  1:788098/1‑149 (MQ=255)
                                                                                                                         ggTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGACGGTTTGATTGTTTACGGCGACGTCAGCTCCACCCACGGTGCTCGTAAAGGTCTGCCGGACACCGCACGGAAAACTGc                      <  2:1300635/149‑1 (MQ=255)
                                                                                                                          gTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACTGCg                     >  1:41301/1‑149 (MQ=255)
                                                                                                                            aTGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGaa                   <  2:614718/149‑1 (MQ=255)
                                                                                                                            aTGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGaa                   <  2:614592/149‑1 (MQ=255)
                                                                                                                             tGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGTTCTGGCGGATACCGCACTGAAAACTGCGAAc                  >  1:922351/1‑149 (MQ=255)
                                                                                                                             tGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCCCTGAAAACTGCGacc                  >  1:803466/1‑147 (MQ=255)
                                                                                                                             tGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAAc                  >  1:263280/1‑149 (MQ=255)
                                                                                                                                 tGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCTGATACCGCACTGAAAACTGCGAACTCCg              <  1:819312/149‑1 (MQ=255)
                                                                                                                                 tGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCg              >  2:517835/1‑149 (MQ=255)
                                                                                                                                  ggTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCgg             <  1:378819/149‑1 (MQ=255)
                                                                                                                                      tGATGGCGAAGCTGGATGGCTCCATCATCGAACCGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCCCTGAAAGCTGCGAACTCCGGTGAc         >  2:1096209/1‑149 (MQ=255)
                                                                                                                                      tGATGGCGAAGCTGGATGGCTCAATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAc         <  1:1249642/149‑1 (MQ=255)
                                                                                                                                        aTGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGCACTTCCGGGAAGGTCTGAACGTACTCCAGTACTTCATCTCCAACCACGGTGCTCGTAAAGGTCTGCCGGATACCGCACTGAAAACTGCGAACGCCGGTTAcgt       >  1:837725/1‑147 (MQ=255)
                                                                                                                                          ggCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCTAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGa     >  2:1284138/1‑149 (MQ=255)
                                                                                                                                             gAAGCGGGACGGCTCCAGCATCGAAACGCCAACCACCGCGAACTTCCGCGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAACACTGCGAACTCCGGTTACCTGACtc  <  2:665898/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTC  >  NC_000913/4187451‑4187740

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