Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,470,250 C→A 5.2% D177Y (GAC→TAC)  fkpA ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,470,2500CA 5.2% 870.4 / 15.2 229D177Y (GAC→TAC) fkpAFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (104/113);  new base A (6/6);  total (110/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 233 total aligned reads displayed.

TAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCTTTACC  >  CP009273/3470105‑3470397
                                                                                                                                                 |                                                                                                                                                   
tAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGa                                                                                                                                                   >  1:312426/1‑148 (MQ=255)
 aaGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGAt                                                                                                                                                  <  2:214007/148‑1 (MQ=255)
 aaGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGAt                                                                                                                                                  <  2:2463735/148‑1 (MQ=255)
  aGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATc                                                                                                                                                 <  2:662042/148‑1 (MQ=255)
   gCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCa                                                                                                                                                >  2:3275846/1‑148 (MQ=255)
    ccAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAg                                                                                                                                               <  1:2083943/148‑1 (MQ=255)
      aGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTAGATCAGCg                                                                                                                                             <  2:1082038/148‑1 (MQ=255)
       gTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGAGAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGt                                                                                                                                            <  1:3214493/148‑1 (MQ=255)
        ttCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTa                                                                                                                                           >  2:3063652/1‑148 (MQ=255)
         tCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTAc                                                                                                                                          <  1:2431890/148‑1 (MQ=255)
          cTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTAcc                                                                                                                                         <  1:1103813/148‑1 (MQ=255)
           tggtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCt                                                                                                                                        <  2:2906859/148‑1 (MQ=255)
           tggtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCt                                                                                                                                        <  2:1645623/148‑1 (MQ=255)
           tggtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCt                                                                                                                                        <  1:1976221/148‑1 (MQ=255)
            ggtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCtt                                                                                                                                       >  2:2619708/1‑148 (MQ=255)
            ggtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCtt                                                                                                                                       >  2:792825/1‑148 (MQ=255)
             gtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCttt                                                                                                                                      <  2:2393810/148‑1 (MQ=255)
             gtggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCttt                                                                                                                                      <  2:3270452/148‑1 (MQ=255)
              tggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTg                                                                                                                                     >  1:1773793/1‑148 (MQ=255)
               ggAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGt                                                                                                                                    <  2:3479570/148‑1 (MQ=255)
                 aaTAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAg                                                                                                                                  >  1:992120/1‑148 (MQ=255)
                    aaCCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTAGATCAGCGTACCTTTGTAGtta                                                                                                                               >  2:840141/1‑147 (MQ=255)
                    aaCCAGTTTGATCTTACCGCATGTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCCAAGATCAGCGAACCTTTGTAGtta                                                                                                                               >  2:3190002/1‑147 (MQ=255)
                       cAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACt                                                                                                                            >  1:1482653/1‑148 (MQ=255)
                         gTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTac                                                                                                                          >  2:877381/1‑148 (MQ=255)
                          tttGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTaca                                                                                                                         >  1:2879185/1‑148 (MQ=255)
                            tGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTacaac                                                                                                                       <  2:2221826/148‑1 (MQ=255)
                              aTCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAg                                                                                                                     <  2:702235/148‑1 (MQ=255)
                               tCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTacaacact                                                                                                                    >  1:1635989/1‑146 (MQ=255)
                               tCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAGt                                                                                                                    >  1:3394884/1‑148 (MQ=255)
                               tCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGt                                                                                                                    <  1:1164193/148‑1 (MQ=255)
                                cTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTa                                                                                                                   >  1:157432/1‑148 (MQ=255)
                                 ttACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTAt                                                                                                                  <  1:2552181/148‑1 (MQ=255)
                                 ttACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTAt                                                                                                                  <  1:2516587/148‑1 (MQ=255)
                                   aCCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCg                                                                                                                <  1:1423808/148‑1 (MQ=255)
                                   aCCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCg                                                                                                                <  1:2374821/148‑1 (MQ=255)
                                       ccTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGt                                                                                                            <  2:3179621/148‑1 (MQ=255)
                                       ccTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATAGCTGt                                                                                                            >  2:248458/1‑148 (MQ=255)
                                         tttCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCt                                                                                                          >  2:183234/1‑148 (MQ=255)
                                           tCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCttt                                                                                                        >  1:882817/1‑148 (MQ=255)
                                            cTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTc                                                                                                       >  2:3463287/1‑148 (MQ=255)
                                            cTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTc                                                                                                       >  2:513672/1‑148 (MQ=255)
                                             ttGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCg                                                                                                      <  1:2673009/148‑1 (MQ=255)
                                              tGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCgg                                                                                                     >  1:1057318/1‑148 (MQ=255)
                                               gATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGt                                                                                                    <  1:1323675/148‑1 (MQ=255)
                                                 tGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGc                                                                                                  <  2:1244755/148‑1 (MQ=255)
                                                 tGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGc                                                                                                  >  2:2081014/1‑148 (MQ=255)
                                                  gTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCt                                                                                                 >  1:3515999/1‑148 (MQ=255)
                                                   ttcttcAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCtt                                                                                                <  1:689448/148‑1 (MQ=255)
                                                   ttcttcAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCtt                                                                                                <  1:2831358/148‑1 (MQ=255)
                                                    tcttcAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTc                                                                                               >  1:356058/1‑148 (MQ=255)
                                                      ttcAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGc                                                                                             <  1:1432348/148‑1 (MQ=255)
                                                          gACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCttt                                                                                         >  1:1546110/1‑148 (MQ=255)
                                                           aCCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACCACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTa                                                                                        >  2:2377213/1‑148 (MQ=255)
                                                             cTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTAcc                                                                                      >  1:2878266/1‑148 (MQ=255)
                                                             cTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTAcc                                                                                      >  2:2942317/1‑148 (MQ=255)
                                                                cTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGc                                                                                   <  1:1955737/148‑1 (MQ=255)
                                                                  gTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCtt                                                                                 >  1:1224482/1‑148 (MQ=255)
                                                                  gTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCtt                                                                                 <  2:2538765/148‑1 (MQ=255)
                                                                  gTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCtt                                                                                 >  2:2187579/1‑148 (MQ=255)
                                                                      aaCCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTctac                                                                             <  1:2284255/148‑1 (MQ=255)
                                                                      aaCCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTctac                                                                             <  1:2796593/148‑1 (MQ=255)
                                                                          cGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTAcc                                                                         <  1:3371272/148‑1 (MQ=255)
                                                                           gggATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCt                                                                        <  1:2262877/148‑1 (MQ=255)
                                                                              aTAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGAt                                                                     <  2:1052695/148‑1 (MQ=255)
                                                                                 acacCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATaaa                                                                  <  1:65696/148‑1 (MQ=255)
                                                                                 acacCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATaaa                                                                  >  2:2772128/1‑148 (MQ=255)
                                                                                  cacCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAac                                                                 <  2:1115177/148‑1 (MQ=255)
                                                                                     cGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAaccag                                                              <  1:1310632/148‑1 (MQ=255)
                                                                                     cGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAaccag                                                              >  2:1711482/1‑148 (MQ=255)
                                                                                     cGTCCAGAAGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAaccag                                                              <  1:3298147/148‑1 (MQ=255)
                                                                                      gTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAaccaga                                                             <  1:877380/148‑1 (MQ=255)
                                                                                        ccAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTAGCCTTTACCGGCTTCTACTACCTGATAAaccacacc                                                           >  1:1859625/1‑148 (MQ=255)
                                                                                           gACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGt                                                        <  2:1139452/148‑1 (MQ=255)
                                                                                              ggAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGa                                                     <  1:948812/148‑1 (MQ=255)
                                                                                                  agaaagCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTAGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAagag                                                 >  2:2839265/1‑148 (MQ=255)
                                                                                                   gaaagCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAgg                                                >  1:1703551/1‑148 (MQ=255)
                                                                                                    aaagCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGt                                               <  2:2316574/148‑1 (MQ=255)
                                                                                                        cGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTGTGTAGTTAACTACAACACTAACGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCAACAACCAGATAAACCAGACCAGTATAAGAGGTTTTc                                           >  2:366278/1‑148 (MQ=255)
                                                                                                         ggTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTca                                          <  1:27289/148‑1 (MQ=255)
                                                                                                                  cGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAACGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTtct                                 <  1:1465844/148‑1 (MQ=255)
                                                                                                                      tgtAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCttt                             <  2:2879185/148‑1 (MQ=255)
                                                                                                                       gtAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTg                            >  1:843691/1‑148 (MQ=255)
                                                                                                                            agTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCaaa                       <  2:1637774/148‑1 (MQ=255)
                                                                                                                            agTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCaaa                       >  2:431635/1‑148 (MQ=255)
                                                                                                                              ttGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAAtt                     >  1:2582077/1‑148 (MQ=255)
                                                                                                                                gTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTtc                   <  2:3451338/148‑1 (MQ=255)
                                                                                                                                gTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTtc                   <  2:784073/148‑1 (MQ=255)
                                                                                                                                 tCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTtct                  <  1:3043029/148‑1 (MQ=255)
                                                                                                                                  cGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTtctc                 >  2:2689863/1‑148 (MQ=255)
                                                                                                                                     aCTCTTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTcgcg              >  1:796681/1‑148 (MQ=255)
                                                                                                                                       tctTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTGTGGCAAATTTCTCGCGGt            >  2:3340909/1‑148 (MQ=255)
                                                                                                                                        ctTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGACTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTa           <  1:545464/148‑1 (MQ=255)
                                                                                                                                        ctTTACCGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGACTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTa           <  1:545377/148‑1 (MQ=255)
                                                                                                                                             ccGTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTTATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTTCTCtt      >  1:2151066/1‑148 (MQ=255)
                                                                                                                                             ccGTAGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCtt      <  2:1528741/148‑1 (MQ=255)
                                                                                                                                               gTCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCTTTa    >  1:3079867/1‑148 (MQ=255)
                                                                                                                                                tCGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCTTTAc   >  1:1272250/1‑148 (MQ=255)
                                                                                                                                                 cGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCTTTAcc  >  2:3117115/1‑148 (MQ=255)
                                                                                                                                                 cGATCAGCGTACCTTTGTAGTTAACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCTTTAcc  >  1:470985/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGTCCAGACGGAAAGAAAGCGGTTCACCACGGGTGTAAGAGTTGTCGAACTCTTTACCGTCGATCAGCGTACCTTTGTAGTTCACTACAACAGTATCGCTGTCTTTCGGTGCTTCGCCTTTACCGGCTTCTACTACCTGATAAACCAGACCAGTTGAAGAGGTTTTCACACCTTTCTCTTTGGCAAATTTCTCGCGGTACTCTTTACC  >  CP009273/3470105‑3470397

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