Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,166 A→T 5.8% intergenic (‑67/+166) thiC ← / ← rsd phosphomethylpyrimidine synthase/regulator of sigma D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,1660AT 5.8% 617.3 / 18.3 189intergenic (‑67/+166)thiC/rsdphosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base A (77/101);  new base T (6/5);  total (83/106)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 192 total aligned reads displayed.

GGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGC  >  NC_000913/4196017‑4196315
                                                                                                                                                     |                                                                                                                                                     
ggAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCa                                                                                                                                                       >  1:246650/1‑150 (MQ=255)
  aaaGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAgg                                                                                                                                                     <  2:901335/150‑1 (MQ=255)
  aaaGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAgg                                                                                                                                                     <  2:679713/150‑1 (MQ=255)
   aaGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAggg                                                                                                                                                    <  1:757689/150‑1 (MQ=255)
     ggCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                                                                                                                                                  <  2:188110/150‑1 (MQ=255)
     ggCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                                                                                                                                                  >  1:601077/1‑150 (MQ=255)
      gCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcg                                                                                                                                                 >  2:1070131/1‑150 (MQ=255)
          tGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgc                                                                                                                                             <  1:881333/150‑1 (MQ=255)
           gCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcg                                                                                                                                            <  2:3834/150‑1 (MQ=255)
            ccTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcga                                                                                                                                           >  1:1242278/1‑150 (MQ=255)
            ccTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcga                                                                                                                                           <  2:1046116/150‑1 (MQ=255)
             cTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcgat                                                                                                                                          >  1:101980/1‑150 (MQ=255)
             cTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcgatgcgat                                                                                                                                          >  2:488933/1‑150 (MQ=255)
              ttCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATa                                                                                                                                         >  1:1265755/1‑150 (MQ=255)
               tCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCCCTGCGgcgatc                                                                                                                                              >  2:20961/1‑143 (MQ=255)
                 cAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCa                                                                                                                                      <  1:713601/150‑1 (MQ=255)
                  aGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAg                                                                                                                                     <  2:1157718/150‑1 (MQ=255)
                  aGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAg                                                                                                                                     <  1:722611/150‑1 (MQ=255)
                      tgtCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTa                                                                                                                                 >  1:389136/1‑150 (MQ=255)
                      tgtCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTa                                                                                                                                 <  1:303787/150‑1 (MQ=255)
                        tCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcga                                                                                                                                                <  1:871333/133‑1 (MQ=255)
                         cGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAAct                                                                                                                              >  1:385729/1‑150 (MQ=255)
                         cGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAAct                                                                                                                              >  1:1399921/1‑150 (MQ=255)
                          gATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActc                                                                                                                             <  2:1265755/150‑1 (MQ=255)
                          gATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActc                                                                                                                             >  1:791527/1‑150 (MQ=255)
                          gATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAActc                                                                                                                             >  1:362279/1‑150 (MQ=255)
                              aaaTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGt                                                                                                                         >  1:548438/1‑150 (MQ=255)
                               aaTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGtt                                                                                                                        >  2:425620/1‑150 (MQ=255)
                                aTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTc                                                                                                                       <  2:1242278/150‑1 (MQ=255)
                                aTGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTc                                                                                                                       <  1:927081/150‑1 (MQ=255)
                                   ttGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCt                                                                                                                    <  2:1407666/150‑1 (MQ=255)
                                   ttGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCt                                                                                                                    >  1:145645/1‑150 (MQ=255)
                                    tGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCtt                                                                                                                   >  2:1379453/1‑150 (MQ=255)
                                        cccGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCa                                                                                                               <  2:303362/150‑1 (MQ=255)
                                        cccGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCa                                                                                                               <  2:865280/150‑1 (MQ=255)
                                         ccGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAg                                                                                                              <  2:1273382/150‑1 (MQ=255)
                                            ggCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTa                                                                                                           <  2:1044430/150‑1 (MQ=255)
                                            ggCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTa                                                                                                           <  2:358099/150‑1 (MQ=255)
                                             gcgcGTTGTTCGCGGCGTGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTTGCTTGTCAGACGACGGATGAATCAAGAGACGATCGCCGCAGGTGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAt                                                                                                          <  2:45449/150‑1 (MQ=255)
                                              cgcgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCtt                                                                                                                   <  1:1196668/140‑1 (MQ=255)
                                              cgcgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAtt                                                                                                         <  1:1083350/150‑1 (MQ=255)
                                                cgTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAg                                                                                                       >  1:202382/1‑150 (MQ=255)
                                                  ttgttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGcc                                                                                                     <  2:951923/150‑1 (MQ=255)
                                                     ttCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGa                                                                                                  >  1:698381/1‑150 (MQ=255)
                                                       cgcgGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATc                                                                                                >  1:267388/1‑150 (MQ=255)
                                                       cgcgGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATc                                                                                                >  2:1236831/1‑150 (MQ=255)
                                                           gcggGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGt                                                                                            <  1:688040/150‑1 (MQ=255)
                                                            ggggTCAGTTTTGTTGCAGACATATCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGtt                                                                                           <  2:771103/149‑1 (MQ=255)
                                                            cggGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGtt                                                                                           <  2:1086000/150‑1 (MQ=255)
                                                                 cAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCcgcg                                                                                      <  2:601077/150‑1 (MQ=255)
                                                                 cAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCcgcg                                                                                      >  1:1174749/1‑150 (MQ=255)
                                                                 cAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCcgcg                                                                                      <  1:1070131/150‑1 (MQ=255)
                                                                     tttGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTACGCGGAtg                                                                                  <  2:404495/150‑2 (MQ=255)
                                                                       tgttgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATccc                                                                                <  2:574952/150‑1 (MQ=255)
                                                                         ttgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTGGCCTGATCAGGTTCCGCGGATCCCGa                                                                              <  1:488933/150‑1 (MQ=255)
                                                                         ttgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGa                                                                              >  2:623437/1‑150 (MQ=255)
                                                                          tgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGaa                                                                             >  1:413194/1‑150 (MQ=255)
                                                                              gACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATaaa                                                                         >  1:291862/1‑150 (MQ=255)
                                                                               aCATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAAc                                                                        <  1:802967/150‑1 (MQ=255)
                                                                                cATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACg                                                                       <  2:1198264/150‑1 (MQ=255)
                                                                                cATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACg                                                                       >  1:21044/1‑150 (MQ=255)
                                                                                 aTAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACgg                                                                      <  1:834629/150‑1 (MQ=255)
                                                                                  tAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCCATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTCATCAGGTTCCGCGGATCCCGAATAAACGGt                                                                     >  1:1267894/1‑150 (MQ=255)
                                                                                   aGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGTAGATTATATATTGTTACCTTCGATTGTTTTAGTCTGATAATGTTTTTCGGATCACGAATaaattttc                                                                    <  2:1259436/150‑6 (MQ=255)
                                                                                     cTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                                  >  1:207678/1‑150 (MQ=255)
                                                                                     cTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                                  >  1:954127/1‑150 (MQ=255)
                                                                                      tCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCa                                                                 >  2:573003/1‑150 (MQ=255)
                                                                                        aTTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGc                                                               >  1:341996/1‑150 (MQ=255)
                                                                                         ttCCAAAAAGTTAAGGACGTGGCTTGTCAGAGGACGGATGAAGCAAGAGACGATCGCCGCATGTGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc                                                              <  1:639877/150‑1 (MQ=255)
                                                                                         ttCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc                                                              >  2:1396933/1‑150 (MQ=255)
                                                                                          tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCa                                                             >  1:49678/1‑150 (MQ=255)
                                                                                           ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAg                                                            <  2:791527/150‑1 (MQ=255)
                                                                                            cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGt                                                           <  2:886260/150‑1 (MQ=255)
                                                                                                aaGTTAAGGACTTGGCTTGTCAGAAGACGGATGAAGCAAGAGACGATCGCCGCAGGGGAGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCAGCGGATCCAGAATAAACGGTCTCAGCCAGTTAAg                                                       >  2:512739/1‑150 (MQ=255)
                                                                                                aaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAg                                                       >  1:250748/1‑150 (MQ=255)
                                                                                                 aGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGTCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGTTATTAGCCTGATCAGGTTCCGCGTATCCAGAATAAACGGTCTCAGCCAGTTAAgg                                                      <  2:1399921/150‑1 (MQ=255)
                                                                                                  gTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGc                                                     <  2:1332678/150‑1 (MQ=255)
                                                                                                      aGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTc                                                 <  1:425620/150‑1 (MQ=255)
                                                                                                      aGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTc                                                 <  1:53832/150‑1 (MQ=255)
                                                                                                       ggACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTcc                                                <  2:555385/150‑1 (MQ=255)
                                                                                                             ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAg                                          >  2:281108/1‑150 (MQ=255)
                                                                                                             ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAg                                          <  2:733958/150‑1 (MQ=255)
                                                                                                              gCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGa                                         <  2:337488/150‑1 (MQ=255)
                                                                                                                  gTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtt                                     >  2:161049/1‑150 (MQ=255)
                                                                                                                   tCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTa                                    <  1:899053/150‑1 (MQ=255)
                                                                                                                    cAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                   >  1:1402920/1‑150 (MQ=255)
                                                                                                                    aaGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCTCAGGTATTAGCCTCATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGACAGTTAAGGCACTCCGACAAGAAATTAg                                   >  1:1048706/2‑150 (MQ=255)
                                                                                                                     aGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGAACTCCGACAAGAAATTAga                                  >  2:1301944/1‑149 (MQ=255)
                                                                                                                          acgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc                             <  2:229088/150‑1 (MQ=255)
                                                                                                                              aTGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAAt                         <  2:709605/150‑1 (MQ=255)
                                                                                                                                gAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg                       <  2:134053/150‑1 (MQ=255)
                                                                                                                                 aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg                      >  2:1389084/1‑150 (MQ=255)
                                                                                                                                   gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc                    <  2:131147/150‑1 (MQ=255)
                                                                                                                                    cAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCa                   >  1:311849/1‑150 (MQ=255)
                                                                                                                                      agagACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGc                                                     <  1:136262/114‑1 (MQ=255)
                                                                                                                                      agagACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGc                                                     >  2:136262/1‑114 (MQ=255)
                                                                                                                                      agagACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAtt                 <  1:1389084/150‑1 (MQ=255)
                                                                                                                                       tatcgcaTCGCCCCTGCGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTg                >  1:452602/7‑150 (MQ=255)
                                                                                                                                         gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGaa              >  2:246814/1‑150 (MQ=255)
                                                                                                                                                     aGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgc  <  1:1379453/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGC  >  NC_000913/4196017‑4196315

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