Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,415 A→C 8.7% L272R (CTG→CGG) ‡ gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,4150AC 8.7% 382.0 / 26.2 150L272R (CTG→CGG) ‡gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (66/70);  new base C (3/10);  total (70/80)
Fisher's exact test for biased strand distribution p-value = 8.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-02

Only 100 of 152 total aligned reads displayed.

CATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATC  >  NC_000913/3963267‑3963564
                                                                                                                                                     |                                                                                                                                                     
cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                      >  2:376615/1‑150 (MQ=255)
cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                      >  1:1547028/1‑150 (MQ=255)
cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                      <  1:892312/150‑1 (MQ=255)
   aaaGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGAt                                                                                                                                                   <  1:187268/150‑1 (MQ=255)
   aaaGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGAt                                                                                                                                                   <  2:809036/150‑1 (MQ=255)
    aaGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGAtt                                                                                                                                                  <  2:444414/150‑2 (MQ=255)
     aGCGGCGTTTAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCg                                                                                                                                                 <  1:888485/150‑1 (MQ=255)
     aGCGGCGTTGAATAATACGCAGCGTACGGCATCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCATGATCg                                                                                                                                                 <  2:1072834/150‑1 (MQ=255)
       cggcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGcc                                                                                                                                               <  2:520247/150‑1 (MQ=255)
        ggcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                                                                                                                                              <  2:1449921/150‑1 (MQ=255)
         gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCTACATGCCGTATACCAGGCCTTCACGCATCGCACCGCCCGCCAGTGTCATACACTGAATCTTCAGTTTGGTAAATTTGGCGATCCTGATCGCCat                                                                                                                                             <  2:267981/150‑2 (MQ=255)
          cgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGa                                                                                                                                            <  2:910933/150‑1 (MQ=255)
            ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttt                                                                                                                                          <  1:1406653/150‑5 (MQ=255)
             tGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCattttt                                                                                                                                         <  1:911184/150‑6 (MQ=255)
             tGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCa                                                                                                                                         >  2:725469/1‑150 (MQ=255)
              gAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCTGTAAAAATGGCGATCTGGATCGCCAGACCAc                                                                                                                                        >  2:852268/1‑150 (MQ=255)
               ggatattcGCAGCCGTACGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACt                                                                                                                                       >  1:716264/8‑150 (MQ=255)
               aaTATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttttac                                                                                                                                       <  1:1028566/150‑8 (MQ=255)
               aaTATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttttac                                                                                                                                       <  1:749076/150‑8 (MQ=255)
               aaTAGTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCGAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAACTGGCGATCCTGATCGCCatttttac                                                                                                                                       <  1:1557209/150‑8 (MQ=255)
                  atTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg                                                                                                                                    >  2:770278/1‑150 (MQ=255)
                  atTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg                                                                                                                                    >  1:1333643/1‑150 (MQ=255)
                  atTACGCAGCGTACGGCTGCGAATATCCTGCGCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg                                                                                                                                    <  2:230477/150‑1 (MQ=255)
                   ttACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCggg                                                                                                                                   <  1:906373/150‑1 (MQ=255)
                    tACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGa                                                                                                                                  <  2:137398/150‑1 (MQ=255)
                        cAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAacac                                                                                                                              <  2:1394058/150‑1 (MQ=255)
                                cGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCa                                                                                                                       <  1:1056624/150‑1 (MQ=255)
                                  gCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCCCTCGGGAACCCTAACGCACg                                                                                                                     >  1:1306342/1‑150 (MQ=255)
                                  gCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACg                                                                                                                     >  2:1198966/1‑150 (MQ=255)
                                  gCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACg                                                                                                                     >  2:513671/1‑150 (MQ=255)
                                  gCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACg                                                                                                                     >  1:1595240/1‑150 (MQ=255)
                                   cggCGAATATCCTGCTCGACGACCAGATGCAACTTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGCCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGt                                                                                                                    <  1:320656/148‑1 (MQ=255)
                                   cTGCGAATATCCTGCTCGACGGCCTGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACCCACGt                                                                                                                    <  2:875394/150‑1 (MQ=255)
                                      cGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTGAGTTCGGTAAAAATGGCGATCAGGAACGCCAGACGACTAGTGAACCCTACCGCACGttac                                                                                                                 >  1:827779/1‑148 (MQ=255)
                                        aaTATCCTGCTCGACGGCCAGATGCAACCTCCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAg                                                                                                               <  2:1508700/150‑1 (MQ=255)
                                         atatCCTGCTCGCCGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGc                                                                                                              <  2:165112/150‑1 (MQ=255)
                                          tatCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCg                                                                                                             >  1:348785/1‑150 (MQ=255)
                                            tCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGCCCATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTc                                                                                                           <  2:331770/150‑1 (MQ=255)
                                             ccTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACAGTGAATATTCAGATCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCTGGAACACTAACGCACGTTCCAGCGTCa                                                                                                          >  1:208455/1‑150 (MQ=255)
                                              cTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAg                                                                                                         >  2:521835/1‑150 (MQ=255)
                                                gCTCGACGGCCATATGCACCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTTAATATACAGTTCGGTAAAAATGGCGATCAGGATCGCCTGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGcc                                                                                                       <  2:1483184/150‑1 (MQ=255)
                                                gCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGcc                                                                                                       <  1:908550/150‑1 (MQ=255)
                                                 cTCGACGGCCATATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCTATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGccc                                                                                                      <  2:845060/150‑1 (MQ=255)
                                                 cTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGccc                                                                                                      >  2:732481/1‑150 (MQ=255)
                                                     aCGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACACACGTTCCAGCGTCAGCCCGTCa                                                                                                  >  1:246615/1‑150 (MQ=255)
                                                       ggCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAt                                                                                                >  1:965409/1‑150 (MQ=255)
                                                        gCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                               >  2:777196/1‑150 (MQ=255)
                                                           aGATGCAACATGCCGTAGACCCGGCCTTCACGCAGCGCACCGCCCGCCCGGGTAATACACTGAATATTCAGTTCGGTAAAAATGCCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTcc                                                                                            <  2:1432918/150‑1 (MQ=255)
                                                           aGATGCAACATGCCGTAGACCCGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTTAATATTCAGTTCGGTAAAATTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATcccc                                                                                            <  1:666827/150‑1 (MQ=255)
                                                           aGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTcc                                                                                            >  2:489983/1‑150 (MQ=255)
                                                            gATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCa                                                                                           <  1:201066/150‑1 (MQ=255)
                                                                  aCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCCATCTCCAGGtcct                                                                                     >  1:576667/1‑150 (MQ=255)
                                                                  aCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcct                                                                                     >  1:440742/1‑150 (MQ=255)
                                                                  aCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcct                                                                                     <  1:406348/150‑1 (MQ=255)
                                                                    aTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTTAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcctcc                                                                                   <  2:1033677/150‑1 (MQ=255)
                                                                    aTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcctcc                                                                                   <  1:725469/150‑1 (MQ=255)
                                                                    aTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcctcc                                                                                   <  1:1474877/150‑1 (MQ=255)
                                                                       ccGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAGCTCCAGTTCCTCCAGc                                                                                >  2:248285/1‑150 (MQ=255)
                                                                        cGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGTGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATATCCAGTTCCTCCAGcc                                                                               >  2:1317858/1‑150 (MQ=255)
                                                                        cGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGGAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGGCAATCTCCAGTTCCTCCAGcc                                                                               >  2:1223811/1‑150 (MQ=255)
                                                                         gTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCg                                                                              >  2:1656191/1‑150 (MQ=255)
                                                                              ccAGGCCTTCACGCCGCGCACCGCCCGCCAGGTTCATACACTGAATATTCAGTTCGGTAAACATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGc                                                                         <  2:134608/150‑1 (MQ=255)
                                                                              ccAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCACCAGTTCCTCCAGCCGACCGc                                                                         >  1:1419733/1‑150 (MQ=255)
                                                                               cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTACCGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCa                                                                        >  2:625094/1‑150 (MQ=255)
                                                                               cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCa                                                                        >  2:194947/1‑150 (MQ=255)
                                                                               cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCa                                                                        >  2:618829/1‑150 (MQ=255)
                                                                                 ggCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaat                                                                      <  2:974576/150‑1 (MQ=255)
                                                                                 ggCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaat                                                                      >  2:1362281/1‑150 (MQ=255)
                                                                                  gCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatg                                                                     <  2:716264/150‑1 (MQ=255)
                                                                                  gCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatg                                                                     <  1:833093/150‑1 (MQ=255)
                                                                                   ccTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatga                                                                    <  2:1113749/150‑1 (MQ=255)
                                                                                        aCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGc                                                               >  2:457362/1‑150 (MQ=255)
                                                                                        aCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGGAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGc                                                               >  2:82279/1‑150 (MQ=255)
                                                                                          gCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTc                                                             <  2:1246514/150‑1 (MQ=255)
                                                                                             gcgcACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGcg                                                          >  1:535182/1‑149 (MQ=255)
                                                                                                 aCCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGttt                                                      <  1:1081765/150‑1 (MQ=255)
                                                                                                  ccgcccgcCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTc                                                     <  1:1198966/150‑1 (MQ=255)
                                                                                                       cgccAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtg                                                >  2:342424/1‑150 (MQ=255)
                                                                                                       cgccAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtg                                                <  1:188915/150‑1 (MQ=255)
                                                                                                          cAGGGTCATACACTGAATATTCAGTTCGGTAACAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgctg                                             <  2:128433/150‑1 (MQ=255)
                                                                                                            gggTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCa                                           >  1:1202407/1‑150 (MQ=255)
                                                                                                             ggTCATACACTTAATATTCAGTTCGTTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCaa                                          <  1:1648249/150‑1 (MQ=255)
                                                                                                               tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACt                                        <  1:376615/150‑1 (MQ=255)
                                                                                                               tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACt                                        <  1:1478447/150‑1 (MQ=255)
                                                                                                                     acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCa                                  >  2:1195904/1‑150 (MQ=255)
                                                                                                                     acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCa                                  >  2:1705767/1‑150 (MQ=255)
                                                                                                                     acTGAATATTCAGTACGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGCGAACACTAACGCACGATCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAAGGAATGGCACGGTGTTTCAATTGCTGCTACTTTTACa                                  >  2:1092508/1‑150 (MQ=255)
                                                                                                                        gAATATTCATTTCGGTAAAAATGGCGAGCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAggg                               <  2:710143/150‑1 (MQ=255)
                                                                                                                          atatTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTa                             >  1:1117925/1‑150 (MQ=255)
                                                                                                                                aGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGt                       <  1:770278/150‑1 (MQ=255)
                                                                                                                                 gTTCGGTAAAAATGGCGATCAGGATCGCCTGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGtt                      <  1:514019/150‑1 (MQ=255)
                                                                                                                                 gTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGtt                      >  2:1320831/1‑150 (MQ=255)
                                                                                                                                  ttCGGGAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTc                     <  1:116298/150‑1 (MQ=255)
                                                                                                                                      gTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGCCCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc                 >  1:359696/1‑150 (MQ=255)
                                                                                                                                          gtctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCatccatcc             >  1:1089850/4‑150 (MQ=255)
                                                                                                                                           tctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATccc            >  1:697436/3‑150 (MQ=255)
                                                                                                                                                cGATCAGGATCGCCAGACCACTCGGGAACACAAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgtg       <  2:1066460/150‑1 (MQ=255)
                                                                                                                                                 gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGCCCGCAATGAATGGCTCGCTGGTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCAGCCATCCCCTGTGc      >  2:361712/1‑150 (MQ=255)
                                                                                                                                                     aGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCcatc  <  1:21897/150‑1 (MQ=255)
                                                                                                                                                     aGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCcatc  >  2:1288715/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATC  >  NC_000913/3963267‑3963564

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: