Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,416 G→T 8.5% L272M (CTG→ATG) ‡ gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,4160GT 8.5% 426.0 / 20.6 148L272M (CTG→ATG) ‡gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (66/69);  new base T (3/10);  total (69/79)
Fisher's exact test for biased strand distribution p-value = 8.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 153 total aligned reads displayed.

CATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCA  >  NC_000913/3963267‑3963565
                                                                                                                                                      |                                                                                                                                                     
cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                       <  1:910720/150‑1 (MQ=255)
cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                       >  1:1547028/1‑150 (MQ=255)
   aaaGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGAt                                                                                                                                                    <  1:187268/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)
       cggcgTTGAATAATACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGcc                                                                                                                                                <  2:1212873/150‑1 (MQ=255)
         gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCat                                                                                                                                              <  2:1395180/150‑2 (MQ=255)
         gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCGGATCGCCAg                                                                                                                                              <  2:1307640/150‑1 (MQ=255)
          cgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt                                                                                                                                             <  2:290813/150‑3 (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:397105/1‑150 (MQ=255)
              gAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCTGTAAAAATGGCGATCTGGATCGCCAGACCAc                                                                                                                                         >  2:852268/1‑150 (MQ=255)
               ggatattcGCAGCCGTACGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACt                                                                                                                                        >  1:716264/8‑150 (MQ=255)
               aaTATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttttac                                                                                                                                        <  1:749076/150‑8 (MQ=255)
               aaTATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttttac                                                                                                                                        <  1:1028566/150‑8 (MQ=255)
               aaTATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACt                                                                                                                                        >  2:561537/1‑150 (MQ=255)
               aaTAGTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCGAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAACTGGCGATCCTGATCGCCatttttac                                                                                                                                        <  1:1557209/150‑8 (MQ=255)
                atatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCCGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc                                                                                                                                       >  2:1564599/1‑150 (MQ=255)
                  atTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCTGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCCCTCgg                                                                                                                                     >  1:902828/1‑150 (MQ=255)
                  atTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGTTCATACACTGAATATTCAGTTCGGTAAAAATGGCGCTCAGGATCGCCAGACCACTCgg                                                                                                                                     >  1:710143/1‑150 (MQ=255)
                  atTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg                                                                                                                                     >  1:1333643/1‑150 (MQ=255)
                   ttACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCTGTTCGGTTAAAATGGCGACGAGGATCGCCTGACCACTCggg                                                                                                                                    >  2:229303/1‑150 (MQ=255)
                   ttACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCggg                                                                                                                                    <  1:906373/150‑1 (MQ=255)
                    tACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGa                                                                                                                                   <  2:137398/150‑1 (MQ=255)
                      cGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAac                                                                                                                                 <  2:1141128/150‑1 (MQ=255)
                                cGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGACCCCTAACGcc                                                                                                                        >  2:873445/1‑149 (MQ=255)
                                cGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCa                                                                                                                        <  1:1056624/150‑1 (MQ=255)
                                  gCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCCCTCGGGAACCCTAACGCACg                                                                                                                      >  1:1306342/1‑150 (MQ=255)
                                  gCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACg                                                                                                                      >  1:1595240/1‑150 (MQ=255)
                                   cggCGAATATCCTGCTCGACGACCAGATGCAACTTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGCCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGt                                                                                                                     <  1:320656/148‑1 (MQ=255)
                                      cGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTGAGTTCGGTAAAAATGGCGATCAGGAACGCCAGACGACTAGTGAACCCTACCGCACGttac                                                                                                                  >  1:827779/1‑148 (MQ=255)
                                          tatCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCg                                                                                                              >  1:348785/1‑150 (MQ=255)
                                            tCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGCCCATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTc                                                                                                            <  2:331770/150‑1 (MQ=255)
                                                gCTCGACGGCCATATGCACCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTTAATATACAGTTCGGTAAAAATGGCGATCAGGATCGCCTGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGcc                                                                                                        <  2:1483184/150‑1 (MQ=255)
                                                gCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGAAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGcc                                                                                                        <  2:1700401/150‑1 (MQ=255)
                                                 cTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGccc                                                                                                       >  2:732481/1‑150 (MQ=255)
                                                        gCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                                >  2:777196/1‑150 (MQ=255)
                                                           cgATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTcc                                                                                             <  1:166672/149‑1 (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)
                                                                 aaCAGGCCGTAGACCAGGCCTTCCCGCTGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTGCGTTAAAAATGGCGATCAGGACCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcc                                                                                       <  1:877043/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)
                                                                   cATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCTATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcctc                                                                                     <  2:760242/150‑1 (MQ=255)
                                                                    aTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTTAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcctcc                                                                                    <  2:1033677/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)
                                                                              ccAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCTTACACTGAATATTCAGTTCGGGAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTACCGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGccgcccgc                                                                          >  1:316131/1‑150 (MQ=255)
                                                                              ccAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCACCAGTTCCTCCAGCCGACCGc                                                                          >  1:1419733/1‑150 (MQ=255)
                                                                               cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTACCGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCa                                                                         >  2:625094/1‑150 (MQ=255)
                                                                               cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCa                                                                         >  2:194947/1‑150 (MQ=255)
                                                                                aGGCCTTCACGCTGCGCTCCGCCCGCCAGGGTCATACACTGATTATTCAGTTCGGTAATAATGGCTATCAGGATCTCCAGACCACTCGGGAAGACTTACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaa                                                                        <  1:988287/150‑1 (MQ=255)
                                                                                 ggCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaat                                                                       >  2:1478447/1‑150 (MQ=255)
                                                                                 ggCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaat                                                                       <  2:974576/150‑1 (MQ=255)
                                                                                  gCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatg                                                                      <  2:716264/150‑1 (MQ=255)
                                                                                     ttCACGCAGCGCACCGCCCGCCTGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGCTCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCAGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaat                                                                   <  2:484718/150‑1 (MQ=255)
                                                                                        aCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGGAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGc                                                                >  2:82279/1‑150 (MQ=255)
                                                                                             gggcACCGCCCGCCAGGGGCATACACTGAATATTCAGTTCGGTAAAAATGGCTATCAGGAGCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCt                                                           <  2:842797/148‑1 (MQ=255)
                                                                                             gcgcACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGcg                                                           >  1:535182/1‑149 (MQ=255)
                                                                                                cACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGtt                                                        >  2:514019/1‑150 (MQ=255)
                                                                                                 aCCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGttt                                                       <  1:1081765/150‑1 (MQ=255)
                                                                                                  ccgcccgcCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTc                                                      <  1:1198966/150‑1 (MQ=255)
                                                                                                    gcccgccAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                                    >  1:196087/1‑150 (MQ=255)
                                                                                                     cccgccAGGGTCATACACTGAATATTCAGTTCGCTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAAt                                                   <  2:1068343/150‑1 (MQ=255)
                                                                                                         ccAGTGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgct                                               <  2:801253/150‑1 (MQ=255)
                                                                                                          cAGGGTCATACACTGAATATTCAGTTCGGTAACAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgctg                                              <  2:128433/150‑1 (MQ=255)
                                                                                                          cAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgctg                                              >  2:681821/1‑150 (MQ=255)
                                                                                                            gggTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCa                                            >  1:1202407/1‑150 (MQ=255)
                                                                                                             ggTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCaa                                           >  2:917713/1‑150 (MQ=255)
                                                                                                               tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACt                                         <  1:376615/150‑1 (MQ=255)
                                                                                                               tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACt                                         <  1:1478447/150‑1 (MQ=255)
                                                                                                                  tACACTGAATATTCAGTTCGGTATAAATGGCGATCAGGATCGCCAGACCCCTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAAGGGCTCGCTGTTTCAATTGCTGCAACtttt                                      >  1:217633/1‑150 (MQ=255)
                                                                                                                    cacTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGCCCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTcc                                    >  2:25749/1‑150 (MQ=255)
                                                                                                                    cacTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTcc                                    <  2:348785/150‑1 (MQ=255)
                                                                                                                     acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCa                                   >  2:1195904/1‑150 (MQ=255)
                                                                                                                     acTGAATATTCAGTACGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGCGAACACTAACGCACGATCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAAGGAATGGCACGGTGTTTCAATTGCTGCTACTTTTACa                                   >  2:1092508/1‑150 (MQ=255)
                                                                                                                       tGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAgg                                 <  1:732481/150‑1 (MQ=255)
                                                                                                                        gAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAggg                                <  2:1333643/150‑1 (MQ=255)
                                                                                                                          atatTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTa                              >  1:1117925/1‑150 (MQ=255)
                                                                                                                               cAGTTCGGTAAAAATGGCGATCAGGATCGCCATACCACTCGGGAACACTAACGCACGGTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCTCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCg                         <  1:820726/150‑1 (MQ=255)
                                                                                                                                 gTTCGGTAAAAATGGCGATCAGGATCGCCTGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGtt                       <  1:514019/150‑1 (MQ=255)
                                                                                                                                  ttCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGATCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTc                      <  2:72141/150‑1 (MQ=255)
                                                                                                                                  ttCGGGAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTc                      <  1:116298/150‑1 (MQ=255)
                                                                                                                                      gTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGCCCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc                  >  1:359696/1‑150 (MQ=255)
                                                                                                                                          gtctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCatccatcc              >  1:1089850/4‑150 (MQ=255)
                                                                                                                                             tGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCt           >  2:1366651/1‑150 (MQ=255)
                                                                                                                                                 gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGCCCGCAATGAATGGCTCGCTGGTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCAGCCATCCCCTGTGc       >  2:361712/1‑150 (MQ=255)
                                                                                                                                                 gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGc       >  2:1256273/1‑150 (MQ=255)
                                                                                                                                                     aGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCcatc   <  1:21897/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCA  >  NC_000913/3963267‑3963565

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: