Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,196,103 C→T 8.1% M217I (ATG→ATA hldE ← heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,196,1030CT 8.1% 301.5 / 8.4 98M217I (ATG→ATAhldEheptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base C (47/43);  new base T (4/4);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTA  >  NC_000913/3195957‑3196250
                                                                                                                                                  |                                                                                                                                                   
aCCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGTATATGCAGCGGCGCTTTACCCGGTTGCAGCCGCGACATACCCTGTTCGGACCGGGTCACTAACAGCGCCGAGAGTTCGTAATCGGCAATCACTTTCATg                                                                                                                                                  >  2:1177059/1‑150 (MQ=255)
 ccGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCGTTACCCGGGTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGc                                                                                                                                                 <  2:1070010/150‑1 (MQ=255)
    tgtCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCcgc                                                                                                                                              >  2:990955/1‑150 (MQ=255)
       cGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCt                                                                                                                                           <  1:1029170/150‑1 (MQ=255)
        gCCCGCACCGGTAACGGCATCCACTTCCTGCGCGTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTc                                                                                                                                          <  2:1051917/150‑1 (MQ=255)
          ccGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTcaa                                                                                                                                        <  1:461730/150‑1 (MQ=255)
          ccGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGGTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTcaa                                                                                                                                        <  1:455373/150‑1 (MQ=255)
           cGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTcaac                                                                                                                                       <  1:1678537/150‑1 (MQ=255)
           cGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTcaac                                                                                                                                       <  1:732307/150‑1 (MQ=255)
           cGCACCGGTAACGTCATACACATCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTcaac                                                                                                                                       <  2:346598/150‑1 (MQ=255)
               ccGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAAtc                                                                                                                                   <  2:582603/150‑1 (MQ=255)
                   tAACGTCATCCACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAAACGGCAATCAGTTTCATGCCGCGCTCAACAATCtctt                                                                                                                               <  2:1282029/150‑1 (MQ=255)
                    aaCGTCATACACTTCCTGCGCTTGGGTGGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCtcttc                                                                                                                              <  2:1480632/150‑1 (MQ=255)
                     aCGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCtcttct                                                                                                                             <  1:1143093/150‑1 (MQ=255)
                      cGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGAGATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAAGCAGTTTCATGCCGCGCTCAACGAtctcagctt                                                                                                                            >  1:276054/1‑145 (MQ=255)
                      cGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTcttctt                                                                                                                            >  2:1324178/1‑150 (MQ=255)
                      cGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTcttctt                                                                                                                            >  1:211286/1‑150 (MQ=255)
                       ggcATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCttcttc                                                                                                                           <  1:595634/148‑1 (MQ=255)
                           tACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCTCTATACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGAAATCAGTTTCATGACGCGCTCAACAATCTCTACATCGGta                                                                                                                       >  1:351886/1‑149 (MQ=255)
                           tACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTc                                                                                                                       >  2:426589/1‑150 (MQ=255)
                           tACACTTCCTGCGCTTGGGGTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTc                                                                                                                       <  1:931325/150‑1 (MQ=255)
                                 tCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTTCCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGGCTTCCAt                                                                                                                 >  2:1707751/1‑150 (MQ=255)
                                   cTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAAGCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTGTTTCAttt                                                                                                               >  2:457336/1‑150 (MQ=255)
                                      cgcTTGGGTTGGAATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTAcc                                                                                                            <  1:1052046/150‑1 (MQ=255)
                                          tGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGaca                                                                                                        <  2:674192/150‑1 (MQ=255)
                                            ggTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGacaac                                                                                                      <  1:491029/150‑1 (MQ=255)
                                             gTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGACTTACATTTACCGacaaca                                                                                                     <  2:9270/150‑1 (MQ=255)
                                               gggCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCTTGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGc                                                                                                   <  2:1162184/149‑1 (MQ=255)
                                                 gCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCtt                                                                                                 <  1:118435/150‑1 (MQ=255)
                                                  cATATGCAGCGGCGCTTTACCTGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTc                                                                                                <  2:404540/150‑1 (MQ=255)
                                                   atatGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCa                                                                                               <  1:835849/150‑1 (MQ=255)
                                                    tatGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCaa                                                                                              >  2:1557836/1‑150 (MQ=255)
                                                     atGCAGCGTCGCTTTACCCGGGTGCAGCAGCGACATACACTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTTATCGGCAATCAGTTTCATGCAGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCaaa                                                                                             <  2:709709/150‑1 (MQ=255)
                                                     atGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCaaa                                                                                             <  1:1547381/150‑1 (MQ=255)
                                                      tGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAAt                                                                                            <  1:544576/150‑1 (MQ=255)
                                                        cAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTc                                                                                          >  1:1306818/1‑150 (MQ=255)
                                                         aGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTcc                                                                                         >  2:191231/1‑150 (MQ=255)
                                                            ggcgCTTTACCCGGTTGCAGCAGCGACATTCCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCgag                                                                                      <  1:1502181/150‑1 (MQ=255)
                                                                     cccGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCGGAGAGCTTCGGc                                                                             >  1:629554/1‑150 (MQ=255)
                                                                     cccGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGc                                                                             >  1:1656681/1‑150 (MQ=255)
                                                                         gTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTa                                                                         <  2:600695/150‑1 (MQ=255)
                                                                                   cGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTATCAGCGTAgc                                                               <  1:1483542/150‑1 (MQ=255)
                                                                                      cATCCCCTGTTCGGAACGGGTCACGACCAGAGCCGAGAGTTCGTACTCGGCAATCAGTTTCATGCCGCGCTCAACAATCGCTTCTTCGGTCTTACATTTACCGACAACAGCGTCAAATTCCGAGCCAGTCGTCGGTAACAGCGTAgcgac                                                            >  2:261375/1‑148 (MQ=255)
                                                                                      cATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACATCAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGcgcc                                                            >  2:1299938/1‑150 (MQ=255)
                                                                                        tACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTTCATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGcgccgc                                                          >  2:27297/1‑150 (MQ=255)
                                                                                            cTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTTCATTTACCGACAACAGCTGCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTa                                                      >  2:105899/1‑150 (MQ=255)
                                                                                             tGTTCGGAACGGGTCACTAACAGAGCCGTGAGTTCGTAATCGGCAATCAGTTTCATTCCGCGCTCAACAATCTCTTCTTCGGGCTTACATTTACCGACAACAGCTTCAAATTCCGTGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAg                                                     <  2:1299868/150‑1 (MQ=255)
                                                                                                 cGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCt                                                 >  1:547916/1‑150 (MQ=255)
                                                                                                  ggAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTc                                                >  1:1603784/1‑150 (MQ=255)
                                                                                                    accGGGTCACAAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCaa                                              <  2:342170/148‑1 (MQ=255)
                                                                                                        gggcACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTATTCTTCTGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATc                                          <  1:1140790/147‑1 (MQ=255)
                                                                                                        ggTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATc                                          >  2:1580681/1‑150 (MQ=255)
                                                                                                         gTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCg                                         >  2:127426/1‑150 (MQ=255)
                                                                                                         gTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCg                                         >  2:123019/1‑150 (MQ=255)
                                                                                                          tCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCgg                                        >  1:660960/1‑150 (MQ=255)
                                                                                                           cACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGGCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGt                                       <  1:303902/150‑1 (MQ=255)
                                                                                                            aCTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATCCCGAGCTCAACAATCTCTTCTTCGGTCTTACATTTTCCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGAGCTCAAAATCGGTa                                      >  1:292479/1‑150 (MQ=255)
                                                                                                             cTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGtga                                     >  1:321174/1‑148 (MQ=255)
                                                                                                               aaCAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGATTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCGCAAAATCGGTACCt                                   >  1:847663/1‑150 (MQ=255)
                                                                                                                 cAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCttt                                 <  2:1017543/150‑1 (MQ=255)
                                                                                                                 cAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCttt                                 >  2:776287/1‑150 (MQ=255)
                                                                                                                    agCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTgg                              <  2:1084182/150‑1 (MQ=255)
                                                                                                                     gCCGAGAGTTCGTACTCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAGCAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGGAGCGCTCAAAATCGGTACCTTGTGGa                             >  1:274172/1‑150 (MQ=255)
                                                                                                                     gCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTGCCTTTTGGa                             <  1:1467435/150‑1 (MQ=255)
                                                                                                                     gCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCATCATCGGTCTTACATTTACCGACATCAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGa                             <  1:204740/150‑1 (MQ=255)
                                                                                                                      ccGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGat                            >  2:207366/1‑150 (MQ=255)
                                                                                                                      ccGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGat                            >  2:385334/1‑150 (MQ=255)
                                                                                                                      ccGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGat                            >  1:814368/1‑150 (MQ=255)
                                                                                                                      ccGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGat                            >  1:892906/1‑150 (MQ=255)
                                                                                                                       gatcttttCGTAATAGGAAATCAGTTTTATGCCGCGATCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           <  2:288107/144‑1 (MQ=255)
                                                                                                                       gatctgTTCGTAATCTGCAATCAGTTTTATTCCGCGATCAAAAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTAGTCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           <  2:291869/145‑1 (MQ=255)
                                                                                                                       cGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTTTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           >  2:307360/1‑150 (MQ=255)
                                                                                                                         agagTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatcaa                         <  2:1494119/150‑1 (MQ=255)
                                                                                                                            gatCGTAATCGGCAATCAGTTTTATGCCGCGCTCAACAATCTCTTCTTCCGTCTTACATTTACCGACAACACCTTCAAATTCCTAGAGATTCGGCCTTAACAGCGTAGCGCCGCGCTAGCGCTCAAAATCGGTACCTTATGGatcagatc                      >  1:291869/3‑146 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTTATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatcagatc                      >  1:288107/1‑146 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTTATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           >  1:1432337/1‑145 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTTATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           <  2:1432337/145‑1 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTTATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           <  2:1425770/145‑1 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTTATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatc                           >  1:1425770/1‑145 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTTGatcaatca                      <  1:1324178/150‑1 (MQ=255)
                                                                                                                            gTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGatcaatca                      >  1:361631/1‑150 (MQ=255)
                                                                                                                             ttCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAg                     >  2:333686/1‑150 (MQ=255)
                                                                                                                                gTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCAc                  >  1:1659198/1‑150 (MQ=255)
                                                                                                                                  aaTCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTACGAGAGATTCGGAGTTAACAGCGTAGCGCCGCGGTAGCGCGCAAAATCGGTACCTTTTTGATCAATCAGCACCg                >  1:1069869/1‑150 (MQ=255)
                                                                                                                                   aTCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCgg               <  1:426589/150‑1 (MQ=255)
                                                                                                                                    tCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGa              <  1:191231/150‑1 (MQ=255)
                                                                                                                                    tCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGa              >  1:1185786/1‑150 (MQ=255)
                                                                                                                                    tCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGa              >  2:1683073/1‑150 (MQ=255)
                                                                                                                                     cGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGaa             >  2:568893/1‑150 (MQ=255)
                                                                                                                                     cGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGaa             >  1:1018486/1‑150 (MQ=255)
                                                                                                                                        cAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCACTCAGCACCGGAacac          >  1:696110/1‑150 (MQ=255)
                                                                                                                                           tCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCg       >  2:166015/1‑150 (MQ=255)
                                                                                                                                              gTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCTTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCCGCACCGGAACCCCCGCtt    >  1:850673/1‑150 (MQ=255)
                                                                                                                                              gTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCtt    >  1:1468536/1‑150 (MQ=255)
                                                                                                                                               tttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCttt   <  1:1408380/150‑1 (MQ=255)
                                                                                                                                               tttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCttt   <  2:660960/150‑1 (MQ=255)
                                                                                                                                                ttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTa  >  1:833347/1‑150 (MQ=255)
                                                                                                                                                ttCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTa  <  2:808285/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTCGGAACGGGTCACTAACAGAGCCGAGAGTTCGTAATCGGCAATCAGTTTCATGCCGCGCTCAACAATCTCTTCTTCGGTCTTACATTTACCGACAACAGCTTCAAATTCCGAGAGATTCGGCGTTAACAGCGTAGCGCCGCGGTAGCGCTCAAAATCGGTACCTTTTGGATCAATCAGCACCGGAACACCCGCTTTA  >  NC_000913/3195957‑3196250

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