Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 553,889 T→C 5.1% I18V (ATC→GTC) ‡ lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913553,8890TC 5.1% 561.8 / 14.2 175I18V (ATC→GTC) ‡lpxHUDP‑2,3‑diacylglucosamine pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (79/87);  new base C (4/5);  total (83/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 178 total aligned reads displayed.

GACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTT  >  NC_000913/553741‑554035
                                                                                                                                                    |                                                                                                                                                  
gACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                    >  4:602292/1‑149 (MQ=255)
     cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCg                                                                                                                                               <  3:724357/149‑1 (MQ=255)
     cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCg                                                                                                                                               >  3:949834/1‑149 (MQ=255)
     cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCg                                                                                                                                               >  2:591378/1‑149 (MQ=255)
        tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGtt                                                                                                                                            <  3:1116567/149‑1 (MQ=255)
        tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGtt                                                                                                                                            <  2:885659/149‑1 (MQ=255)
        tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGtt                                                                                                                                            <  1:553850/149‑1 (MQ=255)
         ttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                           <  4:1151747/149‑2 (MQ=255)
           gATCGCCGCCGTATCTTGCGATGGAGTGGGTTAGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAACATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGttcttc                                                                                                                                        <  4:1406932/149‑1 (MQ=255)
            aTCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGttcttc                                                                                                                                        >  4:1137874/1‑149 (MQ=255)
            aTCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGtt                                                                                                                                            <  4:110145/145‑1 (MQ=255)
            aTCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGtt                                                                                                                                            >  3:110145/1‑145 (MQ=255)
             tCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTcc                                                                                                                                       >  1:1422225/1‑149 (MQ=255)
              cgccgccgcCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCa                                                                                                                                      <  1:737190/149‑1 (MQ=255)
                    cgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCaga                                                                                                                                <  2:245999/149‑1 (MQ=255)
                    cgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCaga                                                                                                                                <  3:141062/149‑1 (MQ=255)
                      ccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAgaga                                                                                                                              <  2:857381/149‑1 (MQ=255)
                      ccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCAATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAgaga                                                                                                                              >  3:308915/1‑149 (MQ=255)
                       cATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGAt                                                                                                                             >  2:749149/1‑149 (MQ=255)
                           ttGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAg                                                                                                                         <  1:604041/149‑1 (MQ=255)
                            tGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGa                                                                                                                        <  2:1029299/149‑1 (MQ=255)
                             gCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                       <  2:505805/149‑1 (MQ=255)
                             gCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                       <  4:308530/149‑1 (MQ=255)
                              cGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATc                                                                                                                      <  4:538666/149‑1 (MQ=255)
                               gATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCt                                                                                                                     >  3:894256/1‑149 (MQ=255)
                                aTGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTg                                                                                                                    <  4:897852/149‑1 (MQ=255)
                                aTGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTg                                                                                                                    <  2:1349132/149‑1 (MQ=255)
                                       gggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaaa                                                                                                             <  4:332756/149‑1 (MQ=255)
                                       gggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaaa                                                                                                             <  3:1092616/149‑1 (MQ=255)
                                         gTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAaga                                                                                                           <  2:1406183/149‑1 (MQ=255)
                                             ggATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtgt                                                                                                       <  1:17117/149‑1 (MQ=255)
                                             ggATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtgt                                                                                                       >  4:768895/1‑149 (MQ=255)
                                              gATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTc                                                                                                      <  4:1303633/149‑1 (MQ=255)
                                                   tcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCAc                                                                                                 >  2:1231494/1‑149 (MQ=255)
                                                     gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGa                                                                                               <  1:1184879/149‑1 (MQ=255)
                                                     gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGa                                                                                               <  3:698139/149‑1 (MQ=255)
                                                      tcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGAt                                                                                              >  2:1172763/1‑149 (MQ=255)
                                                        gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                            >  2:1510343/1‑149 (MQ=255)
                                                        gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                            <  1:1679295/149‑1 (MQ=255)
                                                          cAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACt                                                                                          >  3:1486974/1‑149 (MQ=255)
                                                           aaTCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTc                                                                                         <  1:312432/149‑1 (MQ=255)
                                                            aTCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCg                                                                                        >  4:402173/1‑149 (MQ=255)
                                                            aTCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCg                                                                                        <  3:534779/149‑1 (MQ=255)
                                                                aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                    >  4:222512/1‑149 (MQ=255)
                                                                    ttCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGt                                                                                <  3:1158105/149‑1 (MQ=255)
                                                                    ctCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGt                                                                                <  1:1231494/148‑1 (MQ=255)
                                                                     tCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTc                                                                               <  2:1679080/149‑1 (MQ=255)
                                                                       aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                             >  3:1549827/1‑149 (MQ=255)
                                                                       aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                             >  3:1451539/1‑149 (MQ=255)
                                                                        aaCAGATGGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACg                                                                            >  2:943661/1‑149 (MQ=255)
                                                                                gCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                    <  4:1486974/149‑1 (MQ=255)
                                                                                gCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                    <  4:1717270/149‑1 (MQ=255)
                                                                                gCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                    <  3:1329061/149‑1 (MQ=255)
                                                                                  cAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGa                                                                  >  2:1131572/1‑149 (MQ=255)
                                                                                    aGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATa                                                                <  2:670642/149‑1 (MQ=255)
                                                                                      aaTATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAATATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAAc                                                              <  3:894107/149‑1 (MQ=255)
                                                                                        tataCAGCGCGTCGGCCTTGCGGGTTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGCCACGCTTTAAATGATAACGt                                                            <  3:1137874/149‑1 (MQ=255)
                                                                                              gcgcGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTcttctt                                                      >  2:1095949/1‑149 (MQ=255)
                                                                                              gcgcGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTcttctt                                                      >  2:519605/1‑149 (MQ=255)
                                                                                                gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                    <  1:636771/149‑1 (MQ=255)
                                                                                                gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                    <  1:1172763/149‑1 (MQ=255)
                                                                                                gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAATGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                    <  1:1131572/149‑1 (MQ=255)
                                                                                                 cgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTg                                                   <  1:1317510/149‑1 (MQ=255)
                                                                                                 cgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCCTTTg                                                   <  1:1097405/149‑1 (MQ=255)
                                                                                                  gTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTgg                                                  <  2:1592489/149‑1 (MQ=255)
                                                                                                           gCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTg                                         <  1:1578662/149‑1 (MQ=255)
                                                                                                             gggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGt                                       >  3:447732/1‑149 (MQ=255)
                                                                                                             gggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGt                                       <  4:423975/149‑1 (MQ=255)
                                                                                                             gggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGt                                       >  2:1626046/1‑149 (MQ=255)
                                                                                                             gggCTTCCACCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCGCGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCAAGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGt                                       >  2:220688/1‑149 (MQ=255)
                                                                                                              ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCaa                                                                     >  3:1472337/1‑118 (MQ=255)
                                                                                                                    ccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAcc                                >  1:967720/1‑149 (MQ=255)
                                                                                                                    ccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAcc                                >  3:1057282/1‑149 (MQ=255)
                                                                                                                     cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                               >  4:1436335/1‑149 (MQ=255)
                                                                                                                     cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                               >  2:104321/1‑149 (MQ=255)
                                                                                                                     cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                               <  2:557615/149‑1 (MQ=255)
                                                                                                                      ccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGc                              >  3:1095360/1‑149 (MQ=255)
                                                                                                                        gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCaa                                                                     >  3:869590/1‑108 (MQ=255)
                                                                                                                        gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCaa                                                                     <  4:869590/108‑1 (MQ=255)
                                                                                                                        gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGt                                       <  2:1688061/138‑1 (MQ=255)
                                                                                                                          tAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACg                          >  1:1755267/1‑149 (MQ=255)
                                                                                                                            aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAc                                 <  1:337456/140‑1 (MQ=255)
                                                                                                                            aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAc                                 >  2:337456/1‑140 (MQ=255)
                                                                                                                            aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                        >  4:714546/1‑149 (MQ=255)
                                                                                                                               aaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCgg                     <  4:29940/149‑1 (MQ=255)
                                                                                                                                aCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                    >  4:1038603/1‑149 (MQ=255)
                                                                                                                                  gCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTg                  >  4:956454/1‑149 (MQ=255)
                                                                                                                                     gAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCaa               >  4:1291051/1‑149 (MQ=255)
                                                                                                                                      agaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCa                >  2:194090/3‑147 (MQ=255)
                                                                                                                                      agaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCa                <  1:194090/145‑1 (MQ=255)
                                                                                                                                       aaaCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAaca             >  3:192131/1‑149 (MQ=255)
                                                                                                                                        aaCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAacac            <  3:1291051/149‑1 (MQ=255)
                                                                                                                                         aCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACAcc           >  4:1555209/1‑149 (MQ=255)
                                                                                                                                             gcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTg       >  4:391840/1‑149 (MQ=255)
                                                                                                                                              cggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGa      <  3:505759/149‑1 (MQ=255)
                                                                                                                                              cggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGa      <  1:1226960/149‑1 (MQ=255)
                                                                                                                                               ggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAt     <  3:566850/149‑1 (MQ=255)
                                                                                                                                                  gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtttt  <  1:1485163/149‑1 (MQ=255)
                                                                                                                                                  gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtttt  >  1:140648/1‑149 (MQ=255)
                                                                                                                                                  gATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAtttt  <  3:768895/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTT  >  NC_000913/553741‑554035

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