Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 553,889 T→C 6.4% 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 6.4% 276.6 / 7.4 94I18V (ATC→GTC) ‡lpxHUDP‑2,3‑diacylglucosamine pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (34/54);  new base C (3/3);  total (37/57)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

GGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTC  >  NC_000913/553740‑554038
                                                                                                                                                     |                                                                                                                                                     
ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                       >  2:722529/1‑150 (MQ=255)
 cacacCGCTTTGAGCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGGCGGCCTTGCGGGCTGCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATc                                                                                                                                                      <  1:1602928/149‑1 (MQ=255)
     ccGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCg                                                                                                                                                  <  2:1672306/150‑1 (MQ=255)
      cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCgg                                                                                                                                                 >  1:1649953/1‑150 (MQ=255)
       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGt                                                                                                                                                <  2:466598/150‑1 (MQ=255)
       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGt                                                                                                                                                <  1:1425118/150‑1 (MQ=255)
       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGt                                                                                                                                                <  1:1418537/150‑1 (MQ=255)
       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGt                                                                                                                                                <  2:682728/150‑1 (MQ=255)
         tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAGACGCAGAAAACCGGCGGTGATCGCCGGttc                                                                                                                                              <  1:1033356/150‑1 (MQ=255)
         tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGttc                                                                                                                                              <  2:750315/150‑1 (MQ=255)
         tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                              <  1:1731311/150‑2 (MQ=255)
         tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGttt                                                                                                                                              <  2:1201835/150‑2 (MQ=255)
           tGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGAGGGCTTCACCCGAGATAATACGCAGAAAACCGGCGGTGATCGCCGGtgctt                                                                                                                                            >  1:1503164/1‑150 (MQ=255)
                 ccgccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGc                                                                                                                                      >  1:1663844/1‑150 (MQ=255)
                   gccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCag                                                                                                                                    <  2:1035636/150‑1 (MQ=255)
                   gccgccATCTTGCGATGGAGTGGGGTGGGATCGTCGTCGCCAATCCATGCTTCCAACAGATCGCCAAGAATAAACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCag                                                                                                                                    <  2:1562121/150‑1 (MQ=255)
                       ccATCTTGCGATGGAGTGGGTTGGGTTCGTCGTCGCCAATCCATGCTTCAATCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGAt                                                                                                                                <  2:415251/150‑1 (MQ=255)
                         aTCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGa                                                                                                                              <  1:1178607/150‑1 (MQ=255)
                          tCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGaa                                                                                                                             >  1:1189869/1‑150 (MQ=255)
                           cTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAg                                                                                                                            <  2:1585872/150‑1 (MQ=255)
                              gCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATc                                                                                                                         <  2:1003176/150‑1 (MQ=255)
                               cGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCt                                                                                                                        <  1:1646099/150‑1 (MQ=255)
                                  tGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCCCGCAGAGATGAAGATCTGCa                                                                                                                     >  1:252122/1‑150 (MQ=255)
                                     aGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATa                                                                                                                  >  1:715647/1‑150 (MQ=255)
                                       tGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCAGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaaa                                                                                                                <  1:692615/150‑1 (MQ=255)
                                        gggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGATGATCTGCAATAAag                                                                                                               >  2:23564/1‑150 (MQ=255)
                                        gggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAag                                                                                                               >  1:1255797/1‑150 (MQ=255)
                                         ggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGCTCGCAGGTTCTACCACGCAGAGATGAAGATCTGCAATAAaga                                                                                                              >  1:1226455/1‑150 (MQ=255)
                                           ttGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgt                                                                                                            <  1:1530639/150‑1 (MQ=255)
                                            tGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtg                                                                                                           <  2:1192998/150‑1 (MQ=255)
                                             gggATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtgt                                                                                                          >  2:508781/1‑150 (MQ=255)
                                                 tcgtcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATCTACAGCGCGGCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGcc                                                                                                      <  2:1167148/150‑1 (MQ=255)
                                                   gtcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCAc                                                                                                    >  1:1230018/1‑150 (MQ=255)
                                                      gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGAt                                                                                                 >  1:673013/1‑150 (MQ=255)
                                                      gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACACAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGAt                                                                                                 >  1:735429/1‑150 (MQ=255)
                                                        cgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                               >  2:1555173/1‑150 (MQ=255)
                                                        cgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                               >  2:911878/1‑150 (MQ=255)
                                                         gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTAc                                                                                              >  2:659271/1‑150 (MQ=255)
                                                          ccAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACt                                                                                             >  1:1415528/1‑150 (MQ=255)
                                                             aTCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGc                                                                                          >  1:640050/1‑150 (MQ=255)
                                                              tCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCt                                                                                         <  2:631152/150‑1 (MQ=255)
                                                                cATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGGTaa                                                                                       >  2:1326320/1‑150 (MQ=255)
                                                                      tCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCa                                                                                 >  1:1397862/1‑150 (MQ=255)
                                                                      tCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCa                                                                                 >  1:1410761/1‑150 (MQ=255)
                                                                       cAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCCCGCTTACTCGCTACCGGTCAc                                                                                >  1:692885/1‑150 (MQ=255)
                                                                         aaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGc                                                                              >  1:1639340/1‑150 (MQ=255)
                                                                          aCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCt                                                                             <  2:1107426/150‑1 (MQ=255)
                                                                             taTCGCCAATAATATACAGCGCGTAGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAATACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTc                                                                          <  2:270009/149‑1 (MQ=255)
                                                                             gATCGCCAAGAATATACAGCGCGTCGGCCGTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTc                                                                          <  2:565376/150‑1 (MQ=255)
                                                                             caTCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTc                                                                          <  1:1065706/149‑1 (MQ=255)
                                                                              aTCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCa                                                                         <  2:647900/150‑1 (MQ=255)
                                                                                ctccAAGAATATACAGCGCGGCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                       <  1:138175/148‑1 (MQ=255)
                                                                                 gCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCTGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACTATTACTCGCTAACGGTCACGCTTTCTATg                                                                      <  1:893544/150‑1 (MQ=255)
                                                                                  ccAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGa                                                                     >  1:1645583/1‑150 (MQ=255)
                                                                                  ccAAGAATATACAGCGCGTCGGCATTGCGGGCTTCCCCCGCTAAAAATCGCAGAAAACCGGCGGTGATCGCCTGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGa                                                                     <  2:265707/150‑1 (MQ=255)
                                                                                   cAAGTATATACAGCGCGTCGGCCTTGCGGGCTGACCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCATTAAAGAGTGTCGCCACGATTACTCGCTAACGGGCACGCTTTCAATGAt                                                                    <  2:534073/150‑1 (MQ=255)
                                                                                        atatACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCATTAAAGAGTGTCGCCACGAGTACTCGCTAACGGTCACGCTTTCTATGATAACGt                                                               >  1:1655171/1‑150 (MQ=255)
                                                                                        atatACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGt                                                               >  1:130715/1‑150 (MQ=255)
                                                                                          ataCAGCGCGTCGGCCTTGCGGGCTTGCCCAGCTAAAAAACGCAGATAACCGGCGGTGATGGCCGGTTCTTCCACGCAGAGTTGAAGATCTGCAATAAAGAGTGTCTCCACGATTTCTCGCTAACGGTCACGCTTTCAAAGATAACGtct                                                             <  2:410262/150‑1 (MQ=255)
                                                                                               gcgcGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCttt                                                        >  2:617309/1‑150 (MQ=255)
                                                                                                 gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTg                                                      <  1:887430/150‑1 (MQ=255)
                                                                                                   gTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAATGGTCACGCTTTCAATGATAACGTCTTCTTTTGGc                                                    <  2:75716/150‑1 (MQ=255)
                                                                                                   gTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGc                                                    <  2:1025672/150‑1 (MQ=255)
                                                                                                         cTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTcc                                              >  2:1044055/1‑150 (MQ=255)
                                                                                                         cTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTcc                                              <  1:1496557/150‑1 (MQ=255)
                                                                                                            gCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTgg                                           <  2:731945/150‑1 (MQ=255)
                                                                                                              gggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCCCGCTTTCAATGATAAGGTCTTCTTTTGGCACGTCCTGGTg                                         >  1:1182751/1‑150 (MQ=255)
                                                                                                               ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGc                                        <  2:1304533/150‑1 (MQ=255)
                                                                                                                   tCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATAc                                    <  1:661382/150‑1 (MQ=255)
                                                                                                                    cccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACCATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGATGAATAcc                                   >  2:1415414/1‑150 (MQ=255)
                                                                                                                     ccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                                  <  2:130715/150‑1 (MQ=255)
                                                                                                                      cccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGc                                 >  2:640255/1‑150 (MQ=255)
                                                                                                                       ccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCt                                <  2:438952/150‑1 (MQ=255)
                                                                                                                        cGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTa                               <  2:266250/150‑1 (MQ=255)
                                                                                                                        cGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTa                               <  2:1515038/150‑1 (MQ=255)
                                                                                                                            aaaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATTATAACGGCTGCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                           >  1:151754/1‑150 (MQ=255)
                                                                                                                            aaaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGAc                           <  1:896999/150‑1 (MQ=255)
                                                                                                                              aaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCg                         <  1:1110570/150‑1 (MQ=255)
                                                                                                                               aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCgg                        <  2:1649953/150‑1 (MQ=255)
                                                                                                                               aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCgg                        <  2:1645583/150‑1 (MQ=255)
                                                                                                                               aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCgg                        >  2:717753/1‑150 (MQ=255)
                                                                                                                                aaCGCAGATAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                       <  1:27611/150‑1 (MQ=255)
                                                                                                                                aaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                       >  1:39191/1‑150 (MQ=255)
                                                                                                                                 aCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGtt                      <  1:722529/150‑1 (MQ=255)
                                                                                                                                  gtggaagaaCCGGCGATCACCGCCGGTTATTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGACACGAATACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTg                     >  1:1151462/8‑150 (MQ=255)
                                                                                                                                   gCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGc                    <  2:715647/150‑1 (MQ=255)
                                                                                                                                   ccAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCATTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGc                    <  2:639752/149‑1 (MQ=255)
                                                                                                                                    cAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGGCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCa                   <  2:880386/150‑1 (MQ=255)
                                                                                                                                     gaagaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCaa                  >  1:400147/5‑150 (MQ=255)
                                                                                                                                      gAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAac                 <  1:1252340/150‑1 (MQ=255)
                                                                                                                                      aagaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAAAAAAGAGTGTCGACACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCAACGACCGGTTGCAac                 >  2:1417376/4‑150 (MQ=255)
                                                                                                                                            cggctgTCATCGCCGGTTCTTCCACTCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCAAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCttt           <  2:1503164/150‑1 (MQ=255)
                                                                                                                                            cggcggTGATCGCCGGTTCTTCCACGCAGAGATGAAGCTCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCttt           <  1:1445276/150‑1 (MQ=255)
                                                                                                                                                 gTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGAttt      <  2:1386175/150‑1 (MQ=255)
                                                                                                                                                    aTCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGt   <  2:27010/150‑1 (MQ=255)
                                                                                                                                                     tCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTCATTTTGTc  >  2:184435/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGATTTTGTC  >  NC_000913/553740‑554038

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