Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,181,285 A→T 5.3% K64* (AAG→TAG) ‡ yegV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,2850AT 5.3% 274.1 / 3.2 92K64* (AAG→TAG) ‡yegVputative kinase
Reads supporting (aligned to +/- strand):  ref base A (48/39);  new base T (2/3);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAG  >  NC_000913/2181137‑2181434
                                                                                                                                                    |                                                                                                                                                     
cACGCGTCAGTCGGTGATGGTCGTCGGCGCGGAGGTCATTGACGTGATCGCCGACGCTTATGCCGTCCCCTGGCGTGGGTGAGATATCGAACTGAAACCGAAGAGGCGTAACGTTGGCGGCTGCGACCTGAATATTTCCGTGGCGTGaaa                                                                                                                                                      >  2:1045332/1‑150 (MQ=255)
  cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                                    <  2:613173/150‑1 (MQ=255)
    catCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgcgc                                                                                                                                                  >  2:1224131/3‑150 (MQ=255)
      tCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGCCGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCt                                                                                                                                                >  2:79267/1‑150 (MQ=255)
       cAGTCGGTGATGGGCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTc                                                                                                                                               <  2:171465/150‑1 (MQ=255)
       cAGTCGGTGATGGGCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTc                                                                                                                                               <  2:1681839/150‑1 (MQ=255)
         gTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCgg                                                                                                                                             >  1:539371/1‑150 (MQ=255)
         gTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCgg                                                                                                                                             >  2:140826/1‑150 (MQ=255)
           cGGTGATGGTCGTCGGCGAGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCAGCGCACTGAATATTACCCTGGCGTTAAAGCGCCTCGGCa                                                                                                                                           >  2:164439/1‑150 (MQ=255)
              tGATGGTCGGCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGCaata                                                                                                                                        <  2:182211/150‑4 (MQ=255)
              tGATGGGCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCg                                                                                                                                        <  2:858384/150‑1 (MQ=255)
               gATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGCaatat                                                                                                                                       <  1:1679844/150‑5 (MQ=255)
                 tGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGCaatattc                                                                                                                                     <  1:416534/150‑7 (MQ=255)
                  ggTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGACGTGGCGTTAAAGCGCCTCGGCATCGAAGc                                                                                                                                    >  2:1375546/1‑150 (MQ=255)
                   gtcgtcGGCGCGGCGGCCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGATAAAGCGCCTCGGCATCGAAGCg                                                                                                                                   >  1:103871/1‑150 (MQ=255)
                    gcggcggcgcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCgg                                                                                                                                  <  2:620083/146‑1 (MQ=255)
                      gtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGt                                                                                                                                <  1:182018/150‑1 (MQ=255)
                          gcgcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATg                                                                                                                            <  2:1059773/150‑1 (MQ=255)
                            gcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGcc                                                                                                                          <  2:1562276/150‑1 (MQ=255)
                                cggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGc                                                                                                                      <  1:569838/150‑1 (MQ=255)
                                cggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGc                                                                                                                      >  1:377567/1‑150 (MQ=255)
                                    cATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCt                                                                                                                  <  1:1100243/150‑1 (MQ=255)
                                     aTTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTc                                                                                                                 >  2:1695758/1‑150 (MQ=255)
                                            tGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAgg                                                                                                          >  2:532165/1‑150 (MQ=255)
                                            tGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAgg                                                                                                          <  1:1347428/150‑1 (MQ=255)
                                             gATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAggg                                                                                                         >  1:1546227/1‑150 (MQ=255)
                                                 gCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCgtg                                                                                                     >  1:1361442/1‑150 (MQ=255)
                                                  ccgccgCTGATGCCCTCCCCTGGCGTGGGGGCGAGATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCgtgt                                                                                                    <  1:322828/150‑1 (MQ=255)
                                                    gACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTgg                                                                                                  <  2:1038940/150‑1 (MQ=255)
                                                     aCGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTggg                                                                                                 <  1:1525815/150‑1 (MQ=255)
                                                        cTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCgg                                                                                              >  1:100926/1‑150 (MQ=255)
                                                          tATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGACACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGGTAAAGCGCCGCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgag                                                                                            <  1:32788/150‑1 (MQ=255)
                                                          tATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgag                                                                                            >  1:471012/1‑150 (MQ=255)
                                                             gCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatg                                                                                         >  2:71226/1‑150 (MQ=255)
                                                              cccTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCGCGGTCAGGGCGTGTGGGCGGAgatga                                                                                        >  2:152524/1‑150 (MQ=255)
                                                              cccTCCCCTGGCGTGGGTGCGATATCGAACCGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatga                                                                                        >  1:488483/1‑150 (MQ=255)
                                                                cTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCggg                                                                                                                                 <  1:253354/107‑1 (MQ=255)
                                                                cTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCggg                                                                                                                                 >  2:253354/1‑107 (MQ=255)
                                                                cTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGAtt                                                                                      <  1:1105541/150‑1 (MQ=255)
                                                                cTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGAtt                                                                                      >  2:84630/1‑150 (MQ=255)
                                                                cTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGAtt                                                                                      <  1:668842/150‑1 (MQ=255)
                                                                cTCCCCTGGCGTGGGGGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGAtt                                                                                      <  2:180674/150‑1 (MQ=255)
                                                                 gccccTCGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGGTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTc                                                                                     <  2:323244/149‑1 (MQ=255)
                                                                  ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                                    <  1:199039/150‑1 (MQ=255)
                                                                   cccTGGCGTGGGTCCGAAATCGAAATGAAAAAGCAGAGTGTAAACTTTGGCGGCTGCGCACGGAATATTGCCTTGGCGTTAAAGCGCATCGGCATCGAAGCGGGTAATGAATTGCCTCTCAGTCAGGCCGCGTGGCCGCAGAGGTTTccc                                                                                   >  2:1720938/1‑148 (MQ=255)
                                                                    ccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGAGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCa                                                                                  >  2:23372/1‑150 (MQ=255)
                                                                    ccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGATGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCa                                                                                  >  2:17305/1‑150 (MQ=255)
                                                                        gCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCGGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCATCCg                                                                              <  1:1631617/150‑1 (MQ=255)
                                                                            gggTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCATGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATg                                                                          <  1:593050/150‑1 (MQ=255)
                                                                                gCGATCTCGAACTGAAACAGCAGAGCGGTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCaa                                                                      <  1:164439/150‑1 (MQ=255)
                                                                                     atCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAagag                                                                 >  2:1733260/1‑150 (MQ=255)
                                                                                       cGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAggg                                                               >  2:1626053/1‑150 (MQ=255)
                                                                                        gAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGc                                                              <  2:705763/150‑1 (MQ=255)
                                                                                           cTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTa                                                           <  2:1189661/150‑1 (MQ=255)
                                                                                             gAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAt                                                         >  1:12670/1‑150 (MQ=255)
                                                                                             gAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAt                                                         >  2:1665942/1‑150 (MQ=255)
                                                                                              aaaCAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGGAACCGGAGGGCAAAAGAGGGCTTAATc                                                        >  1:345294/1‑150 (MQ=255)
                                                                                              aaaCAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATc                                                        <  2:889727/150‑1 (MQ=255)
                                                                                              aaaCAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGAGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATc                                                        >  1:805517/1‑150 (MQ=255)
                                                                                                aCAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGTGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAg                                                      <  1:140826/150‑1 (MQ=255)
                                                                                                      gagCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGAt                                                <  1:553627/150‑1 (MQ=255)
                                                                                                      gagCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGGGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGAt                                                >  2:912095/1‑150 (MQ=255)
                                                                                                       agCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATc                                               <  1:1224131/150‑1 (MQ=255)
                                                                                                          gTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGAt                                            >  2:1597222/1‑150 (MQ=255)
                                                                                                           ttAACGTTGGCGGCTGCGCACTGAATATTTCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATa                                           <  1:1549051/150‑1 (MQ=255)
                                                                                                              aCGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGGTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACg                                        <  1:152524/150‑1 (MQ=255)
                                                                                                                gTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGcc                                      >  2:1298590/1‑150 (MQ=255)
                                                                                                                gTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGcc                                      <  2:539371/150‑1 (MQ=255)
                                                                                                                gTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGcc                                      >  2:564297/1‑150 (MQ=255)
                                                                                                                 ttGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCg                                     >  2:1271240/1‑150 (MQ=255)
                                                                                                                   ggcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGaa                                   >  1:1604052/1‑150 (MQ=255)
                                                                                                                   ggcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGaa                                   >  2:688500/1‑150 (MQ=255)
                                                                                                                    gcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGAATCGAAGCGGGTAATGCCTTGCCTCTCGGTCAGGGCGTGTGGGCGGATATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAg                                  >  2:261011/1‑150 (MQ=255)
                                                                                                                     cggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAgg                                 >  2:116267/1‑150 (MQ=255)
                                                                                                                          gcgcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATa                            >  2:64217/1‑150 (MQ=255)
                                                                                                                            gcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGCGGGCTTAATCAGCCTGATCTATAACGCCGAAGGTGATAAc                          >  2:170952/1‑150 (MQ=255)
                                                                                                                                 gAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTg                     <  1:1597222/150‑1 (MQ=255)
                                                                                                                                 gAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTg                     <  1:1557726/150‑1 (MQ=255)
                                                                                                                                  aaTATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTgg                    <  2:1546227/150‑1 (MQ=255)
                                                                                                                                  aaTATAGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTgg                    <  2:297262/150‑1 (MQ=255)
                                                                                                                                     gatGCCGAGGCGCTTTAACGCCACGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGgtgt                 >  1:1408430/3‑150 (MQ=255)
                                                                                                                                     gatGCCGAGGCGCTTTAACGCCACGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGgtgt                 >  1:518008/3‑150 (MQ=255)
                                                                                                                                        gCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTg              >  2:601905/1‑150 (MQ=255)
                                                                                                                                        gCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTg              >  2:669993/1‑150 (MQ=255)
                                                                                                                                        gCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTg              >  1:924244/1‑150 (MQ=255)
                                                                                                                                        gCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGT‑GGCCGAGGTGAATAGCAACCGGATGGCAAAAGAGGGCATAAACAACCAGATCGAAAACCCCGAATGGGAATACGGCTTgtttaagg             >  2:750083/1‑144 (MQ=255)
                                                                                                                                          cGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGgc            <  1:1626053/150‑1 (MQ=255)
                                                                                                                                               cGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATACCGGCTGGTGTCTGGCGCTgg       >  2:1535061/1‑150 (MQ=255)
                                                                                                                                               cGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTgg       <  2:645637/150‑1 (MQ=255)
                                                                                                                                               cGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTCCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTgg       >  2:1145414/1‑150 (MQ=255)
                                                                                                                                               cGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCATGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGCTAACGCCTAAGGTGATAACGGCTGGGGTCTGGCGCggg       >  2:1535066/1‑150 (MQ=255)
                                                                                                                                                gTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGt      >  1:1545623/1‑150 (MQ=255)
                                                                                                                                                  tAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTg    >  2:768082/1‑150 (MQ=255)
                                                                                                                                                    aaGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAg  <  1:1587056/150‑1 (MQ=255)
                                                                                                                                                    aaGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAg  >  1:1630492/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAG  >  NC_000913/2181137‑2181434

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: ‑