Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,185,540 C→T 71.2% P64L (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,185,5400CT71.2% 123.1 / 73.5 99P64L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base T (32/36);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGC  >  NC_000913/4185395‑4185684
                                                                                                                                                 |                                                                                                                                                
gagTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTa                                                                                                                                              <  1:169573/150‑1 (MQ=255)
gagTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGGTCATGGTATTTCGGAGAAGGTAAAAAGCCGGAAACCATCACCTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTa                                                                                                                                              <  1:241107/150‑1 (MQ=255)
gagTTTGATGCGATCAAAATTGCTCTGGCTGCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTa                                                                                                                                              <  2:930346/150‑1 (MQ=255)
  gTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTaaa                                                                                                                                            >  2:937681/1‑150 (MQ=255)
      gATGCTATCAAAATTGCTCTGGGTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGAt                                                                                                                                        <  1:120114/150‑1 (MQ=255)
       aTGCGATCAAAATTGCTCTGGCTTCGCCATACATGATCCGTTCATGGTATTTCGGGGAAGTTAAAAAGCCGGAAAACATCAACTACCGTACGTTCAAACCAGAACTTGACGTCCTTATCTGCGCCCGTATCTTTGGGCGGGTAAAAGAtt                                                                                                                                       <  1:1134847/150‑1 (MQ=255)
        tGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTa                                                                                                                                      <  2:1157604/150‑1 (MQ=255)
          cGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACg                                                                                                                                    <  1:1037071/150‑1 (MQ=255)
          cGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTGTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACg                                                                                                                                    >  1:804205/1‑150 (MQ=255)
                   ttGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCATTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                                           >  2:1289336/1‑150 (MQ=255)
                     gCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGTAACCATCAACTACCGTACGTTGAAACCAGATCGTGACTGCCTTTTCTGCGTCCGTTTCTTTTGGCTGGTAAAAGATTACGAGTGCCTGTGc                                                                                                                         <  2:1177465/150‑1 (MQ=255)
                     gCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGc                                                                                                                         >  1:434287/1‑150 (MQ=255)
                     gCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCtgttc                                                                                                                         <  2:948225/150‑3 (MQ=255)
                     gCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGc                                                                                                                         <  1:732385/150‑1 (MQ=255)
                      ctctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTGTCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCg                                                                                                                        <  2:475363/150‑1 (MQ=255)
                      ctctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTTCGAGTGCCTGTGCg                                                                                                                        >  2:301564/1‑150 (MQ=255)
                      ctctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCg                                                                                                                        >  1:1215982/1‑150 (MQ=255)
                      ctctGGCGTCGCCAGACATGATCCGTTCATGGTCTTTCGGTTACGTTAAAAAGCCGGAACCCATCAACTACCGTACGTTCAAACCAGACCGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCg                                                                                                                        <  2:329386/150‑1 (MQ=255)
                       tctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGGCGGAAACCTTCAACTACCGTACGTTCTAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCgg                                                                                                                       <  1:515519/150‑1 (MQ=255)
                       tctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCgg                                                                                                                       >  2:62740/1‑150 (MQ=255)
                         tGGCTTCTCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAATAAGCCGGTAACCATCAACTACCGTACGTTCAAACTATAACGTGACGGCCTTTTCTGCGCCCTTATCTTTTGGCTGGTAAAAGATTACGAGTGCCTGTGCGGtt                                                                                                                     <  1:61220/150‑2 (MQ=255)
                          ggCTTCGCCAGACATGATCCGTTGATGGTCTTTCGGTGATGTTAAAAAGCCGGAAACCATCAACTACAGTACGTTCAAACCAGAACGTGACGGCCTTGTCTGCGCGCGTATCTGTGGGCGGGTATAAGATTACGAGTGCCTGTGCGGTaa                                                                                                                    >  2:61694/1‑150 (MQ=255)
                          ggCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTaa                                                                                                                    >  1:515723/1‑150 (MQ=255)
                           gCTGCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGATACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCAGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAg                                                                                                                   <  2:802989/150‑1 (MQ=255)
                            cTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCATGTGCGGTAAGt                                                                                                                  >  1:627949/1‑150 (MQ=255)
                               cGCCAGACATGTTCGGTTCATGGTCTTTCGGTGAAGTTAAAACGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACa                                                                                                               <  1:1257947/150‑1 (MQ=255)
                               cGCCAGACATGATCCGTTCAAGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACa                                                                                                               <  1:498394/150‑1 (MQ=255)
                                gCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGCCGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACaa                                                                                                              >  2:828897/1‑150 (MQ=255)
                                 ccAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCATACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAg                                                                                                             <  1:931533/150‑1 (MQ=255)
                                 ccAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCTTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAg                                                                                                             >  2:538012/1‑150 (MQ=255)
                                 ccAGACATGATCCGTTCATGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAgc                                                                                                            >  2:99726/1‑150 (MQ=255)
                                    gACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAgcgc                                                                                                          >  2:372001/1‑150 (MQ=255)
                                      cATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCt                                                                                                        >  1:1051368/1‑150 (MQ=255)
                                      cATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCt                                                                                                        >  1:49920/1‑150 (MQ=255)
                                       gtGATCCGTTCATGGTATTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCGTTTCTGCGCCCGAATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCTGGAAGAACAAGAGCCTg                                                                                                       >  2:913444/2‑150 (MQ=255)
                                       aTGATCCGTTCAAGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTg                                                                                                       >  2:513684/1‑150 (MQ=255)
                                        tGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGa                                                                                                      >  2:1263393/1‑150 (MQ=255)
                                           tCCGTTCATGGTCTTTCGGGGAAGTTAAAAAGCCGGAAACCATCAACTCCCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGAGCCTGAAac                                                                                                   <  2:978631/150‑1 (MQ=255)
                                              gTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTAAGTTCAAAACAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTTGGCCGGAAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAaccacc                                                                                                >  1:1156990/1‑150 (MQ=255)
                                                 cATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTg                                                                                             >  1:393002/1‑150 (MQ=255)
                                                  aTGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTgg                                                                                            <  2:1031586/150‑1 (MQ=255)
                                                      tCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTc                                                                                        <  2:1090073/150‑1 (MQ=255)
                                                        tttCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATGACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCAt                                                                                      >  1:7731/1‑150 (MQ=255)
                                                         ttCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATc                                                                                     <  2:287874/150‑1 (MQ=255)
                                                                aaGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTgaga                                                                              <  1:847192/150‑1 (MQ=255)
                                                                 aGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTGCAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGaa                                                                             >  1:298244/1‑150 (MQ=255)
                                                                 aGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGAGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGATAAGTACAAGCGCCTGACACACCGTGCCGTCATCTGTGAGaa                                                                             >  1:105575/1‑150 (MQ=255)
                                                                     aaaaaGCCGGAAACCAGCAACTACCGTACGTGCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGAGTACGAGTGCCTGTGCGGTACGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgc                                                                         <  1:369260/150‑1 (MQ=255)
                                                                       aaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcgg                                                                       <  2:952459/150‑1 (MQ=255)
                                                                           ccGGTAACCATCATCTACCGTACGTTCAACCCAGAACGGGACGGCCTTTTCTGCGCCCGGTTCTTGGGGCTGGTTTAAGATGACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCAGCTGTGAGAAGTgcggcggt                                                                   <  1:846551/150‑3 (MQ=255)
                                                                           ccGGAAACCATCAACTTCCGTACGTTCAATCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGtt                                                                   <  2:105308/150‑1 (MQ=255)
                                                                            cGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTg                                                                  <  2:271828/150‑1 (MQ=255)
                                                                            cGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTg                                                                  <  1:685181/150‑1 (MQ=255)
                                                                                 aCCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACACGCGCCTGAACCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTg                                                             >  2:349831/1‑150 (MQ=255)
                                                                                  ccATCAACTACCGTACGTTCACACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTTAGTACAAGCGACTGAAACACCGTGGCGTCATCTGTGAGAAGTGCTGCGTTGACGTgc                                                            >  2:245886/1‑149 (MQ=255)
                                                                                  ccATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGa                                                            >  1:561678/1‑150 (MQ=255)
                                                                                  ccATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGAGTACGAGTGACTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGa                                                            >  2:451166/1‑150 (MQ=255)
                                                                                      cAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCa                                                        <  2:991540/150‑1 (MQ=255)
                                                                                      cAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGAGTTGAAGTGACCCa                                                        >  1:856180/1‑150 (MQ=255)
                                                                                       aaCTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGTGAAGTTCGGCGTTGAAGTGCCCCAg                                                       >  2:205844/1‑150 (MQ=255)
                                                                                        aCTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGa                                                      >  1:744053/1‑150 (MQ=255)
                                                                                         cTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGAc                                                     <  1:172580/150‑1 (MQ=255)
                                                                                            ccGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTaa                                                  >  1:144006/1‑150 (MQ=255)
                                                                                             cGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTaaa                                                 >  2:1231448/1‑150 (MQ=255)
                                                                                              gTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAg                                                >  1:976318/1‑150 (MQ=255)
                                                                                              gTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGAGACCCAGACTAAAg                                                >  2:381156/1‑150 (MQ=255)
                                                                                              gTACGTGCAACCCAGAACGTGACGGCCTTTTATGCGCCCGTATCTTTGGGCTTGTTAAAGATTACGATTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGTCGTCATCTGTGAGTATTGCGGCGTTGTAGTTACCCAGACTAAAg                                                <  2:166508/150‑1 (MQ=255)
                                                                                               tACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGt                                               <  1:758419/150‑1 (MQ=255)
                                                                                                 cGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTAc                                             <  1:800787/150‑1 (MQ=255)
                                                                                                 cGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTAc                                             <  2:405221/150‑1 (MQ=255)
                                                                                                  gTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGCCTAAAGTACg                                            >  1:1127488/1‑150 (MQ=255)
                                                                                                   ttCAAACGAGAAAGTGACGGCCTCTTCTGAGCCCGTATCTTTTGGCTGTTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGc                                           <  2:834875/150‑1 (MQ=255)
                                                                                                   ttCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAACCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGCCCCAGACTAAAGTACGc                                           >  2:588447/1‑150 (MQ=255)
                                                                                                         ccAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGAGTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAg                                     >  1:985433/1‑150 (MQ=255)
                                                                                                          cAGAACGCGACGGCCTTTTCTGCGCCCGTATCTTTGCGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGCCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGc                                    <  1:704503/150‑1 (MQ=255)
                                                                                                           aGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCTGTAAGTACAATCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCg                                   <  2:59747/150‑1 (MQ=255)
                                                                                                           aGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCg                                   <  2:915916/150‑1 (MQ=255)
                                                                                                                    cGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACCCCGTGGCGTCATCTGTGAGAAGTGCGGCGTAGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCca                          >  2:333929/1‑150 (MQ=255)
                                                                                                                    cGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCca                          >  2:440277/1‑150 (MQ=255)
                                                                                                                      gCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCcaca                        <  1:861723/150‑1 (MQ=255)
                                                                                                                      gCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCcaca                        <  1:710931/150‑1 (MQ=255)
                                                                                                                       cctagtCTGCGCCCTTATCTTTGGGCGGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGAGCCTGAAACACCGGGGAGGCATCTGTGAGAAGTGCGGCGTTGAAGTTACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACAt                       <  2:1248212/145‑1 (MQ=255)
                                                                                                                        cTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATc                      <  2:144006/150‑1 (MQ=255)
                                                                                                                        cTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATc                      >  2:1096342/1‑150 (MQ=255)
                                                                                                                         ttttCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCg                     <  2:1161598/150‑1 (MQ=255)
                                                                                                                             cTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACt                 >  2:1186924/1‑150 (MQ=255)
                                                                                                                              tGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTg                >  2:749184/1‑150 (MQ=255)
                                                                                                                              tGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAATCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTg                >  2:1035506/1‑150 (MQ=255)
                                                                                                                               gcgcCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTgg               <  2:597209/150‑1 (MQ=255)
                                                                                                                                cgcCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGc              <  1:409883/150‑1 (MQ=255)
                                                                                                                                 gcCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCt             <  2:434287/150‑1 (MQ=255)
                                                                                                                                  cccGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCtt            <  1:372001/150‑1 (MQ=255)
                                                                                                                                  cccGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGAGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCTTGAGCGTATGGGCCACATCGAACTGGCtt            >  1:1208224/1‑150 (MQ=255)
                                                                                                                                    cGTATCTTTGGGACGGTAAAAGATTACGAGTGCCTGTGCGGTTAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGGTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCttac          <  2:925728/150‑3 (MQ=255)
                                                                                                                                      tATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGACGTGAGCGTATGGGCCACATCGAACTGGCTTcccc        >  2:408162/1‑150 (MQ=255)
                                                                                                                                      tATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTcccc        <  2:1051368/150‑1 (MQ=255)
                                                                                                                                        tCTTTGGGCCGGTAATAGATTACGAGTGCCTGTGCGGTTAGTACAAGCGCCTGTAACAGCGTGGCGGCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAATGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGa      <  2:737603/150‑1 (MQ=255)
                                                                                                                                           ttGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTg   <  1:937681/150‑1 (MQ=255)
                                                                                                                                            tGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCACCGACTgc  >  2:34096/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGC  >  NC_000913/4185395‑4185684

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: