Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,926,036 A→C 5.4% E230A (GAA→GCA)  queF → 7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,0360AC 5.4% 428.0 / 9.9 141E230A (GAA→GCA) queF7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)
Reads supporting (aligned to +/- strand):  ref base A (73/60);  new base C (4/4);  total (77/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 141 total aligned reads displayed.

TTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTG  >  NC_000913/2925889‑2926184
                                                                                                                                                   |                                                                                                                                                    
ttGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAAc                                                                                                                                                    <  2:1146652/150‑1 (MQ=255)
   tCAGCCACCTGCTTAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCTAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGt                                                                                                                                                 <  1:1824860/150‑1 (MQ=255)
    cAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGAATCGTGGACGCCAACTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTg                                                                                                                                                >  1:1357400/1‑150 (MQ=255)
      gCCACGTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACATGGTCTCAGTCCGTCATCACCACGAGTTCCACGCACAGTgag                                                                                                                                              >  2:350937/1‑148 (MQ=255)
      gCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCg                                                                                                                                              <  2:1875748/150‑1 (MQ=255)
       ccACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGt                                                                                                                                             <  2:182346/150‑1 (MQ=255)
        cACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTg                                                                                                                                            >  2:531001/1‑150 (MQ=255)
        cACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTg                                                                                                                                            <  1:1601875/150‑1 (MQ=255)
          ccTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGa                                                                                                                                          >  1:1873131/1‑150 (MQ=255)
          ccTGCGGAAATCAACCTGCCTGATCACCCATCCACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTAACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGa                                                                                                                                          <  2:107740/150‑1 (MQ=255)
            tgctgAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAAc                                                                                                                                        <  1:1526514/150‑1 (MQ=255)
               tgAAATCAAACTGCCTGATCACCCATCACCCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGCCGCCAAATTGACAGAGAAAGACTGCTCCGTTACCTGGTCTCATTCCGTCACCACAACGAGTTCCACGATCAGTGCGTGGAACGCa                                                                                                                                     >  1:404902/1‑150 (MQ=255)
               tgAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcg                                                                                                                                     <  2:2006695/150‑4 (MQ=255)
                gAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCAt                                                                                                                                    >  1:1585835/1‑150 (MQ=255)
                gAAATCAAACTGCCAGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcgt                                                                                                                                    <  1:1326656/150‑5 (MQ=255)
                  aaTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcgttg                                                                                                                                  <  1:848692/150‑7 (MQ=255)
                  aaTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGCGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCCCGAACAGTGAGTGGAACGCTTCt                                                                                                                                  >  1:300105/1‑150 (MQ=255)
                  aaTCAAACTGCCCGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCCCAACGAGTTCCACGAACAGTGCGTGGAACGCATCt                                                                                                                                  >  1:1194380/1‑150 (MQ=255)
                        aCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATg                                                                                                                            <  1:2099321/150‑1 (MQ=255)
                        aCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATg                                                                                                                            <  1:1782062/150‑1 (MQ=255)
                           gCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGAcc                                                                                                                         <  2:1670296/150‑1 (MQ=255)
                            ccTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCt                                                                                                                        >  1:531734/1‑150 (MQ=255)
                            ccTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCt                                                                                                                        >  1:1097488/1‑150 (MQ=255)
                              tGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                                                                      >  2:277279/1‑150 (MQ=255)
                                   aCCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGc                                                                                                                 >  2:1516144/1‑150 (MQ=255)
                                   aCCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGc                                                                                                                 <  1:894199/150‑1 (MQ=255)
                                    cccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGCAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCAAGTTTAATGCCCTGTTACGCt                                                                                                                >  2:1976266/1‑150 (MQ=255)
                                    cccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCt                                                                                                                >  1:2022097/1‑150 (MQ=255)
                                    cccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCt                                                                                                                >  2:1976261/1‑150 (MQ=255)
                                    cacATCAAACAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGTAACGCATCTTTAATGACCTGTTACGCt                                                                                                                <  1:1274771/148‑1 (MQ=255)
                                      cATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGATGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                              <  1:1935659/150‑1 (MQ=255)
                                       aTCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCt                                                                                                             >  2:230685/1‑150 (MQ=255)
                                        tCAACCTGATTGGGGTTCGCTCCCAATTCAGTATCGTGGACTCCAAATTGACAGAGAAAAACTGGTGCGTTACCTGGTCTCATTGCGTCATCACAACGAGTTCCACGAACAGTGAGTGGAACGCATCTTTAATGACCTGTTACGCTTCTg                                                                                                            <  1:1019189/150‑1 (MQ=255)
                                         cAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgc                                                                                                           >  2:219799/1‑150 (MQ=255)
                                         cAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTAAGCTTCTgc                                                                                                           >  2:215931/1‑150 (MQ=255)
                                          aaCCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcc                                                                                                          <  1:658690/150‑1 (MQ=255)
                                          aaCCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAAATGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcc                                                                                                          >  1:1211822/1‑150 (MQ=255)
                                           aCCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                                                                         >  1:89817/1‑150 (MQ=255)
                                            ccAGATTGGGGTTCCCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTACACGAACAGTGCGTGGAACGCATCTTTAATGACCTTATACGCTTCTgcaac                                                                                                        >  2:1700496/1‑147 (MQ=255)
                                              aGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGTCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccagcc                                                                                                      >  2:1814533/1‑150 (MQ=255)
                                               gATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGGCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccagcca                                                                                                     <  2:1553562/150‑1 (MQ=255)
                                                aTTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGccagccag                                                                                                    <  2:1182560/150‑1 (MQ=255)
                                                   ggggTTGGGTCCAATTTCAGTATCGTTGACTCCAAATTGACAGAGATAATCAGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaa                                                                                                 <  1:113655/150‑1 (MQ=255)
                                                     ggTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaaaa                                                                                               >  2:2033062/1‑150 (MQ=255)
                                                        tCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTg                                                                                            <  1:940269/150‑1 (MQ=255)
                                                        tCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCGGCGTTACCTGGGCTCATTCCGTCATCACACCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTACTGACCTGTTACGCTTCTGCCAGCCAGAAAAATTg                                                                                            <  1:375275/150‑1 (MQ=255)
                                                           cTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCTGAAAAATTGAGc                                                                                         <  1:1286259/150‑1 (MQ=255)
                                                            tCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCg                                                                                        <  2:695775/150‑1 (MQ=255)
                                                              cAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGtt                                                                                      >  1:1791042/1‑150 (MQ=255)
                                                                  ttCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACg                                                                                  <  1:397955/150‑1 (MQ=255)
                                                                        aTCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGGCTCATTCCGTCATCAAAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGtt                                                                            <  2:1293965/150‑1 (MQ=255)
                                                                         tCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTta                                                                           <  1:845446/150‑1 (MQ=255)
                                                                          cGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGGTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTtat                                                                          <  2:638648/150‑1 (MQ=255)
                                                                            tGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATAc                                                                        <  2:1789414/150‑1 (MQ=255)
                                                                             ggACGCCAAATTGACAGAGAAACACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGCCCTGTTACGCTTCTGCCAGCCAGAAAAAGTGAGCGTTTACGCACGTTATAcc                                                                       >  1:1054479/1‑150 (MQ=255)
                                                                                agacAAATTGACAGAGAAAAACTTCTTCGTTTCCTGGTCTCATTCCGTCATAACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCATCCAGAAAAATTGAGCGTTTACGCACGTTATACCcgt                                                                    <  2:1194380/147‑1 (MQ=255)
                                                                                  ccAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGATCAGTGCGTGGAACGCATCTTTAATGAACTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACTCTCGTTATACCcgtcg                                                                  >  1:1350191/1‑150 (MQ=255)
                                                                                       ttGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCg                                                             <  2:1898376/150‑1 (MQ=255)
                                                                                        tGACAGAGAAAAACTGCTTCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCgg                                                            <  2:1336507/150‑1 (MQ=255)
                                                                                         gACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCggg                                                           >  1:1471826/1‑149 (MQ=255)
                                                                                         gACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGt                                                           >  1:1545029/1‑150 (MQ=255)
                                                                                          aCAGAGAAAAACTGCTGCGTTACCTGGTCTTATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTc                                                          <  1:2086315/150‑1 (MQ=255)
                                                                                           cAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCt                                                         >  2:2062640/1‑150 (MQ=255)
                                                                                            agagAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCTTCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTg                                                        <  2:2106991/150‑1 (MQ=255)
                                                                                             gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCTTGGCGGTCTgg                                                       >  2:1295909/1‑150 (MQ=255)
                                                                                             gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg                                                       >  1:1051516/1‑150 (MQ=255)
                                                                                             gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg                                                       <  2:531734/150‑1 (MQ=255)
                                                                                               taaaaaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGAc                                                     <  2:2094598/149‑1 (MQ=255)
                                                                                               gaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGAc                                                     <  2:1495248/150‑1 (MQ=255)
                                                                                               gaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGAc                                                     >  1:1123973/1‑150 (MQ=255)
                                                                                               gaAAAACTGCTGCGTTACCTGGTCTCAATCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCAATAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGAACGTTATACCCGTCGTGGCGTTCTAGAc                                                     >  2:1646404/1‑150 (MQ=255)
                                                                                                aaaaaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACTAACAGTGCGTGGAACGCATCTTTAATGTCCTGTTACGCTTCTGCCAGCAAGAAAAATTTAGCGTTTACGGACTTTATAACCGTAGTTGCGGTCTGGACa                                                    >  2:160578/1‑150 (MQ=255)
                                                                                                  aaaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACAtt                                                  <  1:1722526/150‑1 (MQ=255)
                                                                                                     ctgctgCGGTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAAc                                               <  2:2045708/150‑1 (MQ=255)
                                                                                                         tgcGTTACCTGGTCTTATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTACCCCgg                                           >  2:1731439/1‑149 (MQ=255)
                                                                                                         tgcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACATGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGt                                           >  1:616301/1‑150 (MQ=255)
                                                                                                          gcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTg                                          >  1:276697/1‑150 (MQ=255)
                                                                                                          gcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTg                                          >  1:2127233/1‑150 (MQ=255)
                                                                                                               aCCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCa                                     >  1:1748030/1‑150 (MQ=255)
                                                                                                                 gtGGTCTCATTCCGTCTTCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACGTGTTACGGTTGTGCCAGCGAGAAAAATAGAGCGTTTGCGCACGTGTTACCCGTCGTGGCGGTCTGGACATTAACCGGTGGCGAAGt                                   >  1:143446/2‑150 (MQ=255)
                                                                                                                   ggTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTaa                                 <  1:1850637/150‑1 (MQ=255)
                                                                                                                      ctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAg                              >  2:1097008/1‑150 (MQ=255)
                                                                                                                        cATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCCGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTTGCGCAGTAATAGCg                            >  2:81479/1‑150 (MQ=255)
                                                                                                                             cGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACTCACGTTATACCCGTCTTGGCGGTATGGACATTAACCCGTGGCGCAGTAATAGCGAtttt                       >  1:1964602/1‑150 (MQ=255)
                                                                                                                             cGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGAtttt                       <  1:1897892/150‑1 (MQ=255)
                                                                                                                              gTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTg                      <  2:1873131/150‑1 (MQ=255)
                                                                                                                              gTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTg                      >  2:763945/1‑150 (MQ=255)
                                                                                                                                catcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCTGTAATAGCGATTTTGTc                    >  1:270120/1‑150 (MQ=255)
                                                                                                                                catcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTc                    <  1:1516144/150‑1 (MQ=255)
                                                                                                                                    acaACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCAt                >  2:1104492/1‑150 (MQ=255)
                                                                                                                                      gatgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCg              >  2:1151616/6‑150 (MQ=255)
                                                                                                                                      aaCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCg              <  1:692482/150‑1 (MQ=255)
                                                                                                                                         gcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAAAAGCGATTTTGTCCCATCGAcc           >  2:571973/3‑150 (MQ=255)
                                                                                                                                         gAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGAcc           <  1:453404/150‑1 (MQ=255)
                                                                                                                                         gAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGAcc           >  2:1908327/1‑150 (MQ=255)
                                                                                                                                          cgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCCTCGTGGAGGTCTGGACATTAACACGTGGCGCAGTAATAGAGATATTGTCCCATCGACca          >  2:1372616/2‑150 (MQ=255)
                                                                                                                                          aGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACca          >  2:1941437/1‑150 (MQ=255)
                                                                                                                                           gTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGCCcac         >  1:202823/1‑150 (MQ=255)
                                                                                                                                               cACGAACAGTGCGTGGAACGCATCTTTAATGAACTGTTACGCTTCTGCCAGCCAGAAAAATAGAGCGTTAACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCCTCGACAACAAGa     >  1:321880/1‑150 (MQ=255)
                                                                                                                                                  gAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTACCATCGACCACAATACTg  >  2:2007508/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGTCAGCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTG  >  NC_000913/2925889‑2926184

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