Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,304,847 T→G N584T (AAC→ACC)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,304,8470TG97.5% 316.0 / ‑1.2 81N584T (AAC→ACC) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (40/39);  total (41/40)
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.00e-01

TGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGG  >  NC_000913/4304701‑4304993
                                                                                                                                                  |                                                                                                                                                  
tGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTg                                                                                                                                                  >  1:850533/1‑149 (MQ=255)
  gCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGt                                                                                                                                              <  1:1152824/151‑1 (MQ=255)
    tttCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGCAAGGGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCATTCAGCGTCAGGGTGAGGt                                                                                                                                              >  2:833064/1‑149 (MQ=255)
     ttcttcGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGtt                                                                                                                                             <  2:189688/149‑1 (MQ=255)
      tcttcGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGtt                                                                                                                                             <  1:980540/148‑1 (MQ=255)
      tcttcGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTAt                                                                                                                                           <  2:561708/150‑1 (MQ=255)
      tcttcGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTAt                                                                                                                                           >  1:818984/1‑150 (MQ=255)
       cttcGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCAAGCTATCGTTCAGCGTCAGGGTGAGGTTATc                                                                                                                                          >  1:1154222/1‑150 (MQ=255)
         tcGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATc                                                                                                                                          <  1:846985/148‑1 (MQ=255)
           gCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGATAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGt                                                                                                                                     >  1:268594/1‑151 (MQ=255)
             ttttACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTg                                                                                                                                   >  2:1072570/1‑151 (MQ=255)
                    aGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCatat                                                                                                                            <  2:422584/151‑1 (MQ=255)
                    aGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCata                                                                                                                             >  1:1084656/1‑150 (MQ=255)
                    aGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCata                                                                                                                             <  1:957927/150‑1 (MQ=255)
                        ccGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGa                                                                                                                         <  2:597825/150‑1 (MQ=255)
                        ccGACATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGaa                                                                                                                        <  2:924154/151‑1 (MQ=255)
                            cATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTc                                                                                                                    <  2:70203/151‑1 (MQ=255)
                             aTTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTc                                                                                                                    >  1:913436/1‑150 (MQ=255)
                             aTTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCa                                                                                                                   <  2:1247349/151‑1 (MQ=255)
                               tttGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAg                                                                                                                  <  1:302829/150‑1 (MQ=255)
                                 tGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGc                                                                                                                <  2:799922/150‑1 (MQ=255)
                                 tGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCt                                                                                                               >  2:456292/1‑151 (MQ=255)
                                  ggCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCt                                                                                                               <  2:1000581/150‑1 (MQ=255)
                                    cTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGat                                                                                                            >  2:1000519/1‑151 (MQ=255)
                                          cGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATAttttt                                                                                                      >  2:1066283/1‑151 (MQ=255)
                                          cGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATAttttt                                                                                                      <  2:557023/151‑1 (MQ=255)
                                           ggTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATAttttt                                                                                                      >  1:615009/1‑150 (MQ=255)
                                             tCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTAcc                                                                                                   >  2:930114/1‑151 (MQ=255)
                                                 cACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCa                                                                                               <  2:524912/151‑1 (MQ=255)
                                                    ggCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCt                                                                                            <  1:595466/151‑1 (MQ=255)
                                                    ggCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCt                                                                                            <  2:704252/151‑1 (MQ=255)
                                                     gCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCtt                                                                                           <  1:1078297/151‑1 (MQ=255)
                                                      cGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCttt                                                                                          >  2:1191447/1‑151 (MQ=255)
                                                         gCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTcgc                                                                                       <  2:1189690/151‑1 (MQ=255)
                                                         gCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTcgc                                                                                       <  1:1122387/151‑1 (MQ=255)
                                                              tCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCCGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTa                                                                                   >  1:926733/1‑150 (MQ=255)
                                                                ttCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCg                                                                                >  2:768473/1‑151 (MQ=255)
                                                                 tCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATc                                                                                 <  2:110230/149‑1 (MQ=255)
                                                                 tCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATc                                                                                 <  2:992738/149‑1 (MQ=255)
                                                                 tCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATc                                                                                 <  2:943947/149‑1 (MQ=255)
                                                                  cACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAg                                                                              <  1:685020/151‑1 (MQ=255)
                                                                   aCGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAgg                                                                             <  1:1152961/151‑1 (MQ=255)
                                                                     ggCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCg                                                                           <  2:691751/151‑1 (MQ=255)
                                                                     ggCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCg                                                                           <  1:521697/151‑1 (MQ=255)
                                                                     ggCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCg                                                                           <  1:114627/151‑1 (MQ=255)
                                                                               gAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCAt                                                                  >  2:847768/1‑150 (MQ=255)
                                                                               gAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCCGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCAt                                                                  >  2:138707/1‑150 (MQ=255)
                                                                                aaaGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATaa                                                                <  1:104083/151‑1 (MQ=255)
                                                                                   gagGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGt                                                             <  1:550836/151‑1 (MQ=255)
                                                                                    agGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTc                                                            >  2:1259135/1‑151 (MQ=255)
                                                                                     ggCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCg                                                           >  1:403742/1‑151 (MQ=255)
                                                                                        gTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAAc                                                        >  1:616041/1‑151 (MQ=255)
                                                                                            gCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCa                                                    <  1:758452/151‑1 (MQ=255)
                                                                                               tGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCAttt                                                 >  1:865221/1‑151 (MQ=255)
                                                                                               tGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCAtt                                                  <  2:91194/150‑1 (MQ=255)
                                                                                                  ggTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTcg                                               >  1:277522/1‑150 (MQ=255)
                                                                                                    ttGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTcgacc                                            >  1:397940/1‑151 (MQ=255)
                                                                                                     tGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTcgaccg                                           >  1:1085962/1‑151 (MQ=255)
                                                                                                       cAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATCTTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTcgaccgac                                         >  2:828634/1‑151 (MQ=255)
                                                                                                        aGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTcgaccgac                                         >  1:143625/1‑150 (MQ=255)
                                                                                                         gCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACAt                                       >  2:818464/1‑151 (MQ=255)
                                                                                                         gCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACAt                                       <  1:299031/151‑1 (MQ=255)
                                                                                                            ttttCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATccc                                    >  2:737145/1‑151 (MQ=255)
                                                                                                            ttttCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATccc                                    >  1:437139/1‑151 (MQ=255)
                                                                                                                  ggTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAAt                              >  2:261389/1‑151 (MQ=255)
                                                                                                                    tAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATgg                            >  2:279961/1‑151 (MQ=255)
                                                                                                                       ttAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGgtgt                         >  1:607064/1‑151 (MQ=255)
                                                                                                                         aaGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTcc                       >  1:213110/1‑151 (MQ=255)
                                                                                                                         aaGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTcc                       <  2:37771/151‑1 (MQ=255)
                                                                                                                         aaGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTcc                       >  2:298212/1‑151 (MQ=255)
                                                                                                                           gCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCgg                     <  2:1002435/151‑1 (MQ=255)
                                                                                                                            cACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGACGTTGCTCATATTGAAGTTCTGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCggg                    >  1:1003104/1‑151 (MQ=255)
                                                                                                                             aCGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCggg                    <  1:1214496/150‑1 (MQ=255)
                                                                                                                                cTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAAc                >  2:581892/1‑151 (MQ=255)
                                                                                                                                  aTCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCg              >  1:487717/1‑151 (MQ=255)
                                                                                                                                    cGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCgg             <  2:297874/150‑1 (MQ=255)
                                                                                                                                    cGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCgg             >  1:1159222/1‑150 (MQ=255)
                                                                                                                                       tCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCggtgg          <  2:79211/150‑1 (MQ=255)
                                                                                                                                         aGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGc       >  1:1222545/1‑151 (MQ=255)
                                                                                                                                           cGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCg     >  1:444820/1‑151 (MQ=255)
                                                                                                                                              cAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTgg  <  1:1072506/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGG  >  NC_000913/4304701‑4304993

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