Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,721,925 T→C 5.2% V658A (GTG→GCG)  patZ → peptidyl‑lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,721,9250TC 5.2% 231.9 / 2.6 76V658A (GTG→GCG) patZpeptidyl‑lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base C (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

TCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTG  >  NC_000913/2721779‑2722073
                                                                                                                                                  |                                                                                                                                                    
tCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCa                                                                                                                                                   <  1:844908/150‑1 (MQ=255)
   tGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCTATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGt                                                                                                                                                >  1:1219649/1‑150 (MQ=255)
    gAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGtt                                                                                                                                               >  2:1361342/1‑150 (MQ=255)
    gAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGAGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGtt                                                                                                                                               >  1:142536/1‑150 (MQ=255)
      aGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCTGGTTTc                                                                                                                                             <  2:1479766/150‑1 (MQ=255)
         tCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCaa                                                                                                                                          >  1:1735618/1‑150 (MQ=255)
         tCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCaa                                                                                                                                          >  1:1622709/1‑150 (MQ=255)
         tCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCaa                                                                                                                                          <  2:447171/150‑1 (MQ=255)
           aaGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACt                                                                                                                                        >  2:729435/1‑150 (MQ=255)
           aaGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACt                                                                                                                                        <  1:1326202/150‑1 (MQ=255)
            aGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACtt                                                                                                                                       <  2:76067/150‑1 (MQ=255)
            aGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACtt                                                                                                                                       >  2:165882/1‑150 (MQ=255)
                    tcgCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCg                                                                                                                               >  1:1747473/1‑150 (MQ=255)
                     cgCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGa                                                                                                                              >  1:1030968/1‑150 (MQ=255)
                     cgCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGa                                                                                                                              <  2:1222313/150‑1 (MQ=255)
                       cACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGAtt                                                                                                                            >  2:84099/1‑150 (MQ=255)
                        aCTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCTATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTg                                                                                                                           <  2:1219649/150‑1 (MQ=255)
                        aCTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTg                                                                                                                           <  2:1061415/150‑1 (MQ=255)
                         cTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGc                                                                                                                          <  2:431211/150‑1 (MQ=255)
                            ccgccgCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCg                                                                                                                       <  2:1803171/150‑1 (MQ=255)
                             cgccgcTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCgg                                                                                                                      >  2:1171485/1‑150 (MQ=255)
                                  cTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTc                                                                                                                                >  1:882908/1‑135 (MQ=255)
                                  cTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTc                                                                                                                                <  2:882908/135‑1 (MQ=255)
                                  cTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAAtt                                                                                                                 >  2:1020798/1‑150 (MQ=255)
                                     aaCATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAg                                                                                                              <  1:188500/150‑1 (MQ=255)
                                      aCATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGc                                                                                                             >  1:321331/1‑150 (MQ=255)
                                      aCATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGATACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGc                                                                                                             >  2:1560618/1‑150 (MQ=255)
                                             ccTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCAGAAATTCAGCGTCTGGa                                                                                                      <  1:659392/150‑1 (MQ=255)
                                                  cccGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTc                                                                                                 <  1:174473/150‑1 (MQ=255)
                                                      cTATCTTGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTTGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATcc                                                                                             <  1:1950876/150‑1 (MQ=255)
                                                          cTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTCGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTg                                                                                         >  1:804815/1‑150 (MQ=255)
                                                           tGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGc                                                                                        <  1:1657865/150‑1 (MQ=255)
                                                           tGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCAATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGc                                                                                        >  2:1207323/1‑150 (MQ=255)
                                                             gTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTg                                                                                      <  2:924077/150‑1 (MQ=255)
                                                                 ttCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCtt                                                                                  >  1:719656/1‑150 (MQ=255)
                                                                  tCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTc                                                                                 >  2:604029/1‑150 (MQ=255)
                                                                   cAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCt                                                                                >  2:1570413/1‑150 (MQ=255)
                                                                          tCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTg                                                                         >  2:1365554/1‑150 (MQ=255)
                                                                           cAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGa                                                                        >  2:367817/1‑150 (MQ=255)
                                                                           cAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGa                                                                        >  2:374984/1‑150 (MQ=255)
                                                                           cAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGa                                                                        >  1:1521591/1‑150 (MQ=255)
                                                                            aaaaGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGCGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCCTCCTTTGCTGGCTTCTGGCAGTGaa                                                                       >  2:1930001/1‑150 (MQ=255)
                                                                            aaaaGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGCGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGaa                                                                       >  2:1923080/1‑150 (MQ=255)
                                                                                  aaaaaGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTAcc                                                                 >  1:1706288/1‑150 (MQ=255)
                                                                                      aGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGCGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCgcgc                                                             <  1:1930001/150‑1 (MQ=255)
                                                                                      aGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGCGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCgcgc                                                             <  1:1923080/150‑1 (MQ=255)
                                                                                          tCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGa                                                         <  2:553241/150‑1 (MQ=255)
                                                                                              gcgcgcAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTc                                                     <  2:1138932/150‑1 (MQ=255)
                                                                                                    aGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTg                                               <  1:84099/150‑1 (MQ=255)
                                                                                                            aCGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATAtctc                                       >  1:1584583/1‑150 (MQ=255)
                                                                                                             cGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCg                                      <  1:218436/150‑1 (MQ=255)
                                                                                                             cGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGACGATTGCCCGGAAATTCACCGTCTCGATATTCAACCTTTGCTGGCTTCTGGCACTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCg                                      >  2:125730/1‑150 (MQ=255)
                                                                                                                  aTTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCCGGATATCTCGCCGtt                                 >  2:185789/1‑150 (MQ=255)
                                                                                                                    tGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTg                               <  1:1510505/150‑1 (MQ=255)
                                                                                                                     ggATGTCGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGa                              <  2:804815/150‑1 (MQ=255)
                                                                                                                         gTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGc                          >  2:1153643/1‑150 (MQ=255)
                                                                                                                           tGCAGGCTTGAGCCAGCTTGTGGTGCAGGTTTACTAGTTGATTGTCGATTGACCGGAAATTCAGCGTCTGGATATTCATCATTTGCTGGCTTCTGTCAGTGAATTTACCGCGCTGGATGTCACTCTTTATATCTCGACGTTTTAAGGCGa                        <  2:1706288/150‑1 (MQ=255)
                                                                                                                             cAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTg                                                           <  2:1439250/113‑1 (MQ=255)
                                                                                                                             cAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTg                                                           >  1:1439250/1‑113 (MQ=255)
                                                                                                                              aGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATaa                     <  1:1361342/150‑1 (MQ=255)
                                                                                                                               ggCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGTATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAAc                    >  1:847780/1‑150 (MQ=255)
                                                                                                                               ggCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAAc                    <  1:83347/150‑1 (MQ=255)
                                                                                                                                 cTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGagatatc                  <  1:1020798/150‑7 (MQ=255)
                                                                                                                                  ttGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACgag                 <  1:1570413/150‑1 (MQ=255)
                                                                                                                                   tGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACgaga                <  1:729435/150‑1 (MQ=255)
                                                                                                                                    gAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGagag               >  2:212054/1‑150 (MQ=255)
                                                                                                                                    gAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGagag               <  1:643250/150‑1 (MQ=255)
                                                                                                                                      gCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTGATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTc             >  2:1538196/1‑150 (MQ=255)
                                                                                                                                          gCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCt         <  2:1030968/150‑1 (MQ=255)
                                                                                                                                            ttCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTgg       <  1:1171485/150‑1 (MQ=255)
                                                                                                                                              cTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCa     >  2:1718791/1‑150 (MQ=255)
                                                                                                                                              cTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCa     >  1:169744/1‑150 (MQ=255)
                                                                                                                                                ggTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGt   >  2:1587018/1‑150 (MQ=255)
                                                                                                                                                ggTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGt   >  2:1270329/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTg  >  1:1402737/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCCGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTg  <  1:185789/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTG  >  NC_000913/2721779‑2722073

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