Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 630,669 T→A M259K (ATG→AAG)  cstA → carbon starvation protein involved in peptide utilization; APC peptide transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913630,6690TA97.2% 407.1 / 1.8 106M259K (ATG→AAG) cstAcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (54/49);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 106 total aligned reads displayed.

CTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGCC  >  NC_000913/630521‑630816
                                                                                                                                                    |                                                                                                                                                   
cTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGc                                                                                                                                                   >  2:533470/1‑151 (MQ=255)
  gcgcTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCtg                                                                                                                                                 >  2:172654/1‑151 (MQ=255)
       aTCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCtggtg                                                                                                                                              <  2:1259815/149‑1 (MQ=255)
       aTCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTggg                                                                                                                                            >  2:632782/1‑151 (MQ=255)
          tGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGtt                                                                                                                                          >  2:129628/1‑150 (MQ=255)
           gCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGtt                                                                                                                                          <  2:1557767/149‑1 (MQ=255)
           gCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTa                                                                                                                                         <  2:408121/150‑1 (MQ=255)
             gTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGGAGCTGGTGGGTTACgg                                                                                                                                      >  1:918066/1‑151 (MQ=255)
             gTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACgg                                                                                                                                      >  2:1512859/1‑151 (MQ=255)
              tCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGt                                                                                                                                     >  1:181495/1‑151 (MQ=255)
                    gCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTtgt                                                                                                                                >  2:348661/1‑150 (MQ=255)
                      gTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTgg                                                                                                                              >  2:1418508/1‑150 (MQ=255)
                        aTTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTggcgg                                                                                                                           >  1:1102781/1‑151 (MQ=255)
                         ttGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGACACCGTACTATGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTggcgg                                                                                                                           >  1:151891/1‑150 (MQ=255)
                          tGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTggcggc                                                                                                                          <  2:582437/150‑1 (MQ=255)
                          tGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTggcgg                                                                                                                           >  2:579125/1‑149 (MQ=255)
                           ggCGAAGTGTCGGTCATCGGTCTGGTATTCCCGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGcggcgg                                                                                                                        <  1:788960/151‑1 (MQ=255)
                            gCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGGTACGGTTTTGTGGcggcgg                                                                                                                        >  2:1132713/1‑150 (MQ=255)
                               aaGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGtgct                                                                                                                    <  2:1364280/151‑1 (MQ=255)
                                 gtgtCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGtgctg                                                                                                                   <  1:1410449/150‑1 (MQ=255)
                                 gtgtCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGtgct                                                                                                                    <  1:166814/149‑1 (MQ=255)
                                     cGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCgg                                                                                                               <  2:71423/150‑1 (MQ=255)
                                       gTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGGTGCCGgtg                                                                                                             >  2:645750/1‑150 (MQ=255)
                                           tCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGTGGtt                                                                                                        >  2:1600580/1‑151 (MQ=255)
                                            cGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGtt                                                                                                        <  2:214594/150‑1 (MQ=255)
                                            cGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGTgggt                                                                                                        >  2:36134/1‑148 (MQ=255)
                                               tCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGGGGGTAc                                                                                                      >  2:118188/1‑149 (MQ=255)
                                                  ggTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTTCCGGGGTGGGTActgc                                                                                                   >  2:323174/1‑149 (MQ=255)
                                                   gTATTCCTGATTTGCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGGCCTGCGCACCGTGCTTTGACTTTACCGGCGTGCAGTTGGCCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTActgct                                                                                                  <  1:48951/149‑1 (MQ=255)
                                                    tATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTgg                                                                                                >  2:407243/1‑150 (MQ=255)
                                                    tATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTTct                                                                                                     >  1:624050/1‑145 (MQ=255)
                                                    tATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGGTACTGCTgg                                                                                                >  1:550591/1‑150 (MQ=255)
                                                      ttCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTTCCGGGGTGGTTACTGCTGGccc                                                                                             >  2:369153/1‑151 (MQ=255)
                                                         cTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTTCCGGGGTGGGTACTGCTGGCCCcgc                                                                                          >  2:1395919/1‑151 (MQ=255)
                                                           gATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCcgc                                                                                          <  1:879264/149‑1 (MQ=255)
                                                           gATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGGTACTGCTGGcccc                                                                                            >  2:1379958/1‑147 (MQ=255)
                                                            aTTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGGGGGTACTGCTGGGCCcgc                                                                                          >  1:825992/1‑148 (MQ=255)
                                                            aTTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGGTACGGTTTTGTGGCGGCGGGGCTGCCGGGGTGGGTACTGCTGGgccccc                                                                                          >  1:1080893/1‑146 (MQ=255)
                                                              tttCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGa                                                                                      <  2:532992/150‑1 (MQ=255)
                                                              tttCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGGGGTTACTGCTGGcccc                                                                                            >  2:1706669/1‑144 (MQ=255)
                                                                tCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGTGGGTACTGCTGGCCCCCCGTGACt                                                                                    >  1:1583982/1‑150 (MQ=255)
                                                                    cATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCCGGGGTGCCGGTGTGGGTACTGCTGGCCCCGCGGGACttac                                                                                 >  1:757518/1‑146 (MQ=255)
                                                                      ttATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGGTGCCGGGGTGGGTACTGCTGGGCCCCCGGGGCTACctc                                                                               >  1:1241592/1‑149 (MQ=255)
                                                                       tATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACctctc                                                                             <  1:1046279/150‑1 (MQ=255)
                                                                            ctGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGGGGGTACTGCTGGGCCCCCGGGGCTACCTCTCTAcc                                                                         >  1:906374/1‑149 (MQ=255)
                                                                            ctGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGGGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCtt                                                                       >  2:430694/1‑151 (MQ=255)
                                                                              ggcggcTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCt                                                                        <  2:683985/148‑1 (MQ=255)
                                                                               gcggcTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                                                    <  1:1198993/151‑1 (MQ=255)
                                                                                     gggTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAAt                                                              <  2:658909/151‑1 (MQ=255)
                                                                                         ggCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCgg                                                           <  2:364366/150‑1 (MQ=255)
                                                                                            aGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACt                                                       <  1:1238619/151‑1 (MQ=255)
                                                                                            aGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCGGCCGGTGTGGGTACTTCGTGCCCCGCGTGACTACCTCTCTACCTTTCTGAAAATCGGGACt                                                       >  1:47960/1‑151 (MQ=255)
                                                                                              aaaGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAt                                                     <  2:1080978/151‑1 (MQ=255)
                                                                                              aaaGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTTGTTTCTGCCTGCCCCGCGGGGCTACCTCTCTTCCTTCCTGAAAAACgggtata                                                      >  1:272188/1‑146 (MQ=255)
                                                                                                 gTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCCGTGGGGGTTCTGCTTGCCCCCGGTGGGTACCTCTCTACCTTCCTGGAAAACCGGACt                                                       >  1:82070/1‑146 (MQ=255)
                                                                                                 gTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGTGGGTTCTGCTGGCCCCCCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                                  >  1:1359512/1‑151 (MQ=255)
                                                                                                   ccGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgtcg                                                <  1:93922/151‑1 (MQ=255)
                                                                                                   ccGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgtcg                                                <  1:1132625/151‑1 (MQ=255)
                                                                                                   ccGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgtc                                                 <  2:376211/150‑1 (MQ=255)
                                                                                                       ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                                            <  1:1134427/151‑1 (MQ=255)
                                                                                                        cTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGGGTGGGTACTGCTGGCCCCGCGTGGCTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                                               >  1:520065/1‑147 (MQ=255)
                                                                                                         tGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGGTACTGCTGGGCCCGCGGGACTACCCCTCTACCTTCCTGAAAATCGGGGCTATCGTCgg                                               >  1:776290/1‑146 (MQ=255)
                                                                                                          gggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCCGGGTGGTTACTGCTGGGCCCGCGGGGCTACCTCTCTACCCTCCTGAAAAACGGGGCTATCGTCGGTCTgg                                          >  1:1465860/1‑150 (MQ=255)
                                                                                                            gCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCgg                                       <  1:1659705/151‑1 (MQ=255)
                                                                                                             cACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGt                                      <  1:836745/151‑1 (MQ=255)
                                                                                                               ccGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGATTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTGCCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTa                                     <  1:371548/150‑1 (MQ=255)
                                                                                                                   aCTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAt                                <  2:1200274/151‑1 (MQ=255)
                                                                                                                   aCTTTGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGGGGGTACTTCTGGCCCCGCGTGACTACCTCTCTACCTTTCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAt                                >  2:1373283/1‑151 (MQ=255)
                                                                                                                      ttGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAttt                              <  2:826056/150‑1 (MQ=255)
                                                                                                                       tGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGGTTCTGCTGGCCCCGCGGGGCTACCTCTCTACCCTCCTGGAAATCGGGACCATCGTCGGGCTGGCGGGAGGCCtttt                             >  1:325147/1‑150 (MQ=255)
                                                                                                                       tGACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTAATGCTGGCCCCGCGTGGCTACCTCTCTACCTTCCTGAAAATCGGGACTATCGGCGGTCTGGCGGGAGGGAtttt                             >  1:226881/1‑150 (MQ=255)
                                                                                                                        gACTTTACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTg                            <  1:129754/150‑1 (MQ=255)
                                                                                                                            ttACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAtt                         <  1:369126/149‑1 (MQ=255)
                                                                                                                            ttACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAtt                         <  2:181508/149‑1 (MQ=255)
                                                                                                                            ttACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAtt                         <  1:645693/149‑1 (MQ=255)
                                                                                                                            ttACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTAGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAtt                         <  1:1361707/149‑1 (MQ=255)
                                                                                                                             tACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt                       <  2:1426393/150‑1 (MQ=255)
                                                                                                                             tACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGGTGCCGGTGTGGTTACTGCTGGCCCCCCGGGACTACCTCTCTACCCTCCCGAAAAACGGGACCATCGTCGGGCTGGCGGGAGGCATTTTGAtt                         >  2:791242/1‑148 (MQ=255)
                                                                                                                             tACCGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTTCCGGTGGGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGGCTGGCGGTAGGCATTTTtattat                       >  2:1182286/1‑150 (MQ=255)
                                                                                                                                cGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgcgc                   <  1:856057/151‑1 (MQ=255)
                                                                                                                                cGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGGCAtttttgtttttcg                    >  1:512426/1‑141 (MQ=255)
                                                                                                                                cGGCGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGGGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGGGGCATTTTTTTTATgc                     >  1:1144967/1‑149 (MQ=255)
                                                                                                                                   cGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCg                <  2:437727/151‑1 (MQ=255)
                                                                                                                                   cGTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCg                <  1:348634/151‑1 (MQ=255)
                                                                                                                                    gTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGa               <  1:96257/151‑1 (MQ=255)
                                                                                                                                    gTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGGTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGGCATTTTGATTATGCGccc                 >  2:590060/1‑149 (MQ=255)
                                                                                                                                    gTGCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCCGCCGGTGTGGTTACTGCTGGCCCCCCGTGACTACCTCTCTACCTTCCTGAAAATCCGGACTATCGTCGGGCTGGCGGGAAGCATTTTGATTATTGGccc                 >  2:1389374/1‑149 (MQ=255)
                                                                                                                                      gCAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGAc              <  1:1685820/150‑1 (MQ=255)
                                                                                                                                      gCAGTTGACCTGGAAGCTGGCGGGTTACGGTTTTGTGGCGGCGGTGGTGCCGGTGTGGTTAATGCCGGCCCCCCGTGACTACCCCTCTCCCTTCCTGTAAATCGGGACTATCGTCGGTCTGGCGGGAGGCCtttttttttttcccccga               >  1:653965/1‑135 (MQ=255)
                                                                                                                                       cAGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGc            >  2:1352963/1‑151 (MQ=255)
                                                                                                                                        aGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCt           <  2:889225/151‑1 (MQ=255)
                                                                                                                                        aGTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGGGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGGCTGGCGGGAGGCATTTTTGTTATGCGccc                 >  2:777618/1‑145 (MQ=255)
                                                                                                                                         gTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTg          <  1:443225/151‑1 (MQ=255)
                                                                                                                                         gTTGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTg          <  2:1667285/151‑1 (MQ=255)
                                                                                                                                           tGACCTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGa         <  2:1396936/150‑1 (MQ=255)
                                                                                                                                              ccTGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCa      <  1:1535390/150‑1 (MQ=255)
                                                                                                                                                tGGAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCCTTTTTATTTTGCGCCCGGCGCTGAcc       >  1:463323/1‑147 (MQ=255)
                                                                                                                                                 ggAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTTATTATGCGCCCGACGCTGACCATGcc  >  1:76809/1‑151 (MQ=255)
                                                                                                                                                 ggAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATg    <  2:550641/149‑1 (MQ=255)
                                                                                                                                                 ggAAGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGcc  <  2:696438/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGCC  >  NC_000913/630521‑630816

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑