Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 266,175 A→G Y127H (TAC→CAC)  yafY ← lipoprotein, inner membrane; degP regulator; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,1750AG99.0% 379.2 / ‑3.9 100Y127H (TAC→CAC) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (45/54);  total (45/55)
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.66e-01

TAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCC  >  NC_000913/266030‑266322
                                                                                                                                                 |                                                                                                                                                   
tAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCa                                                                                                                                                 >  1:211799/1‑150 (MQ=255)
 aaaGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCaa                                                                                                                                                >  1:1295496/1‑150 (MQ=255)
     cTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTa                                                                                                                                            <  1:809098/150‑1 (MQ=255)
       cgcgGGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTAttt                                                                                                                                         <  1:1225368/151‑1 (MQ=255)
            gCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCtt                                                                                                                                    >  1:1317804/1‑151 (MQ=255)
              aGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCttt                                                                                                                                   <  1:266539/150‑1 (MQ=255)
              aGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTg                                                                                                                                  >  2:753401/1‑151 (MQ=255)
                tACCCCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTgg                                                                                                                                 <  1:401004/150‑1 (MQ=255)
                tACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGt                                                                                                                                >  1:209793/1‑151 (MQ=255)
                  cacCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGtat                                                                                                                              >  2:780648/1‑151 (MQ=255)
                  cacCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGtat                                                                                                                              <  1:753479/151‑1 (MQ=255)
                    cccGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGtata                                                                                                                             <  2:172546/150‑1 (MQ=255)
                    cccGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGtata                                                                                                                             <  1:234496/150‑1 (MQ=255)
                    cccGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGtata                                                                                                                             <  2:646518/150‑1 (MQ=255)
                    cccGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATAc                                                                                                                            <  1:150262/151‑1 (MQ=255)
                        ccAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGt                                                                                                                        <  2:819453/151‑1 (MQ=255)
                        ccAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGt                                                                                                                        >  2:233616/1‑151 (MQ=255)
                          aGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTcc                                                                                                                      <  1:1192449/151‑1 (MQ=255)
                           gCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCg                                                                                                                     <  2:624467/151‑1 (MQ=255)
                            cATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGt                                                                                                                    >  1:1379857/1‑151 (MQ=255)
                            cATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGt                                                                                                                    >  2:694497/1‑151 (MQ=255)
                               aaaCCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCa                                                                                                                  <  1:678387/150‑1 (MQ=255)
                                aaCCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCAcc                                                                                                                >  1:212488/1‑151 (MQ=255)
                                   cGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCAcctcc                                                                                                             <  2:1057655/151‑1 (MQ=255)
                                     aCTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTg                                                                                                           >  1:419237/1‑151 (MQ=255)
                                     aCTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTg                                                                                                           <  2:272119/151‑1 (MQ=255)
                                         aGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTAt                                                                                                       <  1:711506/151‑1 (MQ=255)
                                          gTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTAtt                                                                                                      <  2:667769/151‑1 (MQ=255)
                                          gTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTAtt                                                                                                      >  2:770469/1‑151 (MQ=255)
                                            gAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTcc                                                                                                    >  1:915055/1‑151 (MQ=255)
                                            gAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTcc                                                                                                    >  2:1452427/1‑151 (MQ=255)
                                             aaaGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTcc                                                                                                    >  1:970741/1‑150 (MQ=255)
                                             aaaGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCa                                                                                                   <  2:159950/151‑1 (MQ=255)
                                             aaaGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCa                                                                                                   <  1:1323477/151‑1 (MQ=255)
                                              aaGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTcc                                                                                                    >  2:1403857/1‑149 (MQ=255)
                                              aaGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTcc                                                                                                    >  1:1410173/1‑149 (MQ=255)
                                              aaGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAg                                                                                                  <  2:1022963/151‑1 (MQ=255)
                                                 tATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGttt                                                                                               >  2:826403/1‑151 (MQ=255)
                                                  aTTGTCATTCTTGTCTCCTTTTTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTg                                                                                              <  1:314813/151‑1 (MQ=255)
                                                  aTTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTg                                                                                              <  2:788489/151‑1 (MQ=255)
                                                    tGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGAt                                                                                            >  2:1317605/1‑151 (MQ=255)
                                                       cATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGtt                                                                                         <  1:71166/151‑1 (MQ=255)
                                                       cATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGtt                                                                                         <  2:985470/151‑1 (MQ=255)
                                                       cATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGtt                                                                                         >  2:458625/1‑151 (MQ=255)
                                                        aTTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGtt                                                                                         >  2:1369739/1‑150 (MQ=255)
                                                        aTTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTa                                                                                        <  1:234560/151‑1 (MQ=255)
                                                         ttCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTAt                                                                                       >  2:1413187/1‑151 (MQ=255)
                                                         ttCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTAt                                                                                       >  1:408264/1‑151 (MQ=255)
                                                          tCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATc                                                                                      >  2:1433467/1‑151 (MQ=255)
                                                           cTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCg                                                                                     >  2:528797/1‑151 (MQ=255)
                                                                  ccTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATac                                                                              >  2:989222/1‑151 (MQ=255)
                                                                  ccTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATa                                                                               <  2:212464/150‑1 (MQ=255)
                                                                   cTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATac                                                                              >  1:1174063/1‑150 (MQ=255)
                                                                   cTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATa                                                                               <  2:1057026/149‑1 (MQ=255)
                                                                      tGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATacac                                                                            >  2:636090/1‑149 (MQ=255)
                                                                      tGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATacac                                                                            >  1:387266/1‑149 (MQ=255)
                                                                       gTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAgc                                                                         >  2:367488/1‑151 (MQ=255)
                                                                         cGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAgcg                                                                        <  1:760449/150‑1 (MQ=255)
                                                                             gttTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATtc                                                                   >  1:660648/1‑151 (MQ=255)
                                                                             gttTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATtc                                                                   <  2:1465503/151‑1 (MQ=255)
                                                                              ttttATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATtc                                                                   >  2:805615/1‑150 (MQ=255)
                                                                                 tATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATtcatc                                                                >  1:472571/1‑150 (MQ=255)
                                                                                 tATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATtcatc                                                                <  1:872880/150‑1 (MQ=255)
                                                                                  aTTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATtcatca                                                               <  2:1270105/150‑1 (MQ=255)
                                                                                  aTTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAg                                                              <  2:1175393/151‑1 (MQ=255)
                                                                                   ttCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGt                                                             >  1:1372421/1‑151 (MQ=255)
                                                                                     ccTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTgg                                                           >  1:81528/1‑151 (MQ=255)
                                                                                     ccTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTgg                                                           >  2:444158/1‑151 (MQ=255)
                                                                                        gATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCt                                                         <  2:93441/150‑1 (MQ=255)
                                                                                        gATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTg                                                        >  2:563623/1‑151 (MQ=255)
                                                                                          ttGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTgg                                                       >  1:1305853/1‑150 (MQ=255)
                                                                                              cAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCGTTGGTATACGGGTCCGTCACCTCCTGGTATTCCGGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACgg                                                   <  1:1445804/150‑1 (MQ=255)
                                                                                               aGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGtt                                                 >  2:1044360/1‑151 (MQ=255)
                                                                                                gggTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGttt                                                <  2:1182567/151‑1 (MQ=255)
                                                                                                    tCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTc                                            >  2:1208865/1‑151 (MQ=255)
                                                                                                    tCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTc                                            >  1:709509/1‑151 (MQ=255)
                                                                                                     cGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTcc                                           <  1:1413294/151‑1 (MQ=255)
                                                                                                       aGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTg                                         <  2:211775/151‑1 (MQ=255)
                                                                                                       aGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTg                                         <  2:1416038/151‑1 (MQ=255)
                                                                                                         ggTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCa                                       <  1:887865/151‑1 (MQ=255)
                                                                                                          gTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCat                                      <  1:1433572/151‑1 (MQ=255)
                                                                                                             tCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATc                                   >  1:363345/1‑151 (MQ=255)
                                                                                                              cAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATc                                   <  2:800452/150‑1 (MQ=255)
                                                                                                               aGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATcc                                  >  2:1152954/1‑150 (MQ=255)
                                                                                                               aGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCa                                 <  2:1295391/151‑1 (MQ=255)
                                                                                                                      gCCCCATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGgt                           <  1:60692/150‑1 (MQ=255)
                                                                                                                       ccccATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGgtag                         <  2:213894/151‑1 (MQ=255)
                                                                                                                       ccccATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGgta                          <  2:842195/150‑1 (MQ=255)
                                                                                                                        cccATCAGTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGgtagt                        <  1:241124/151‑1 (MQ=255)
                                                                                                                               gTTTGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACt                 <  1:978023/151‑1 (MQ=255)
                                                                                                                                 ttGCCGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTg                <  2:381237/150‑1 (MQ=255)
                                                                                                                                    ccGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCAt            <  1:226431/151‑1 (MQ=255)
                                                                                                                                    ccGTCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCAt            <  2:209770/151‑1 (MQ=255)
                                                                                                                                       tCATGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCt         >  1:481907/1‑151 (MQ=255)
                                                                                                                                         aTGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCg       <  1:780727/151‑1 (MQ=255)
                                                                                                                                         aTGGGCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCg       <  2:1267084/151‑1 (MQ=255)
                                                                                                                                           gggCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCa     <  1:1448692/151‑1 (MQ=255)
                                                                                                                                           gggCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCa     <  2:583290/151‑1 (MQ=255)
                                                                                                                                            ggCGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAg    <  1:1458604/151‑1 (MQ=255)
                                                                                                                                              cGTGCTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGcc  >  1:284119/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAAAGCTCGCGAGCAGTACACCCGCCAGCATAAACCGACTTAGTGAAAGTATTGTCATTCTTGTCTCCTTTGTCGTTGTTTTATTCCTGATTGTCAGGGTTCGAGGGTATCAGTAGCTGCCCCATCAGTTTGCCGTCATGGGCGTACTCAAAGTATTTTTCTTTGGTATACGGGTCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCC  >  NC_000913/266030‑266322

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