Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,077,153 A→C 5.6% I248S (ATC→AGC)  PP_2687 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,077,1530AC 5.6% 270.5 / 4.3 92I248S (ATC→AGC) PP_2687hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/43);  new base C (2/4);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAGTCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCG  >  NC_002947/3077005‑3077299
                                                                                                                                                    |                                                                                                                                                  
gAAGTCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGAt                                                                                                                                                   <  1:2039827/150‑1 (MQ=255)
  aGTCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGGCCCTCGATGc                                                                                                                                                 >  2:1721498/1‑150 (MQ=255)
   gTCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCa                                                                                                                                                >  2:122522/1‑150 (MQ=255)
    tCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCat                                                                                                                                               >  2:866213/1‑149 (MQ=255)
    tCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAg                                                                                                                                               <  1:2140211/150‑1 (MQ=255)
     cGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAgg                                                                                                                                              <  2:2086648/150‑1 (MQ=255)
      gCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGc                                                                                                                                             >  1:477138/1‑150 (MQ=255)
      gCTTAGCCAGGGATGCAACCCTGAATCTAGGAGTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGc                                                                                                                                             >  1:1270800/1‑150 (MQ=255)
       cTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGCCCCCCCCCCGTTCCTCGATGCAGGCa                                                                                                                                            >  1:456923/1‑150 (MQ=255)
       cTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCa                                                                                                                                            >  2:848865/1‑150 (MQ=255)
       cTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCa                                                                                                                                            >  1:1907005/1‑150 (MQ=255)
        ttAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCAt                                                                                                                                           >  1:35996/1‑150 (MQ=255)
         tAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATg                                                                                                                                          <  2:671084/150‑1 (MQ=255)
            ccAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCag                                                                                                                                       >  1:246363/1‑150 (MQ=255)
               gggATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGTATTAATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcga                                                                                                                                    <  1:760767/150‑1 (MQ=255)
               gggATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcga                                                                                                                                    <  2:1648064/150‑1 (MQ=255)
               gggATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcga                                                                                                                                    <  1:866213/150‑1 (MQ=255)
                  aTGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCagcgagcg                                                                                                                                 <  1:200672/150‑1 (MQ=255)
                   tGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGAGCGCGTTCAGTTCGGCGGCCAGGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCAGGCagcgggagc                                                                                                                                >  2:798268/1‑145 (MQ=255)
                    gCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGcgcg                                                                                                                               <  1:1490871/150‑1 (MQ=255)
                    gCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGcgcg                                                                                                                               <  1:922785/150‑1 (MQ=255)
                        cccTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCAc                                                                                                                           <  1:706047/150‑1 (MQ=255)
                           tGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTc                                                                                                                        <  2:304307/150‑1 (MQ=255)
                              aTCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCAc                                                                                                                     <  2:477138/150‑1 (MQ=255)
                               tCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCAcc                                                                                                                    >  1:230777/1‑150 (MQ=255)
                               tCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCAGTACTGCCCGGCAGGCCTGCAAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCAcc                                                                                                                    >  2:1029554/1‑150 (MQ=255)
                                cGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCAccc                                                                                                                   <  1:1479943/150‑1 (MQ=255)
                                 gAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCa                                                                                                                                                >  2:989429/1‑120 (MQ=255)
                                 gAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCa                                                                                                                                                <  1:989429/120‑1 (MQ=255)
                                 gAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCg                                                                                                                  >  2:2037322/1‑150 (MQ=255)
                                  aGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTCCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGc                                                                                                                 <  1:1721498/150‑1 (MQ=255)
                                       ttCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAAcag                                                                                                            <  2:768243/150‑1 (MQ=255)
                                           tcgATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACagcagc                                                                                                        >  2:608303/1‑150 (MQ=255)
                                                       tAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGAAAAGATCTCAACGGCTTGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACAcc                                                                                            <  2:1270800/150‑1 (MQ=255)
                                                          cgcAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACAccggc                                                                                         >  1:2234730/1‑150 (MQ=255)
                                                           gcAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAg                                                                                                                                               >  2:1604185/1‑95 (MQ=255)
                                                           gcAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAg                                                                                                                                               <  1:1604185/95‑1 (MQ=255)
                                                             aGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACAccggccgg                                                                                      <  1:2391662/150‑1 (MQ=255)
                                                               gcgcgcGTTCAGTTCTGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGANCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGgc                                                                                    <  1:798268/150‑1 (MQ=255)
                                                                         aGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTgcg                                                                          <  2:1831825/150‑1 (MQ=255)
                                                                          gTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCg                                                                                                  >  2:241299/1‑125 (MQ=255)
                                                                          gTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCg                                                                                                  <  1:241299/125‑1 (MQ=255)
                                                                           ttCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCACGCCAGCGATGATCGCGTCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTgcggc                                                                        >  1:1748413/1‑150 (MQ=255)
                                                                           ttCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTgcggc                                                                        >  1:1230296/1‑150 (MQ=255)
                                                                               gcggcCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGa                                                                    >  1:535307/1‑150 (MQ=255)
                                                                                      tGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCt                                                             >  1:1087981/1‑150 (MQ=255)
                                                                                         gcTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGt                                                          >  2:558209/1‑150 (MQ=255)
                                                                                           tGGGCCGCCAGCCCGGTACTCCCCGGCAGGCCTGCCAGGTTCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGTGCGCGGCACGTCCACCCGCAACAGGAGCAATTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTg                                                        <  2:173766/150‑1 (MQ=255)
                                                                                            gggCGGCACGCCCGGTACTGCCCGGCAGACCTGCCAGGATCAGGTCCAGTTCCTCGATGAAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGa                                                       >  2:48569/1‑150 (MQ=255)
                                                                                             ggcggcACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGaa                                                      <  1:608303/150‑1 (MQ=255)
                                                                                              gcggcACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAg                                                     >  2:1035161/1‑150 (MQ=255)
                                                                                                  cACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCtt                                                 >  1:240685/1‑150 (MQ=255)
                                                                                                         ggTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGa                                          <  2:456923/150‑1 (MQ=255)
                                                                                                         ggTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGa                                          <  1:1035161/150‑1 (MQ=255)
                                                                                                         ggTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGa                                          <  1:122522/150‑1 (MQ=255)
                                                                                                          gTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAg                                         <  1:1029554/150‑1 (MQ=255)
                                                                                                              tGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCACGCATGCAGCAAGCGCGGCACGTCCACCCGCAACAGTAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCa                                     <  1:1144470/150‑1 (MQ=255)
                                                                                                                     cAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAgcgcgcg                              <  2:35996/150‑1 (MQ=255)
                                                                                                                        gCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCgcag                           >  1:2140334/1‑150 (MQ=255)
                                                                                                                          cTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGa                                                                    >  2:408382/1‑107 (MQ=255)
                                                                                                                          cTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGa                                                                    <  1:408382/107‑1 (MQ=255)
                                                                                                                           tGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGaa                                                      >  1:1861680/1‑120 (MQ=255)
                                                                                                                           tGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGaa                                                      <  2:1861680/120‑1 (MQ=255)
                                                                                                                           tGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGcagcag                        <  1:848865/150‑1 (MQ=255)
                                                                                                                            gCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGgcgc                                                                      >  2:2356519/1‑103 (MQ=255)
                                                                                                                            gCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGgcgc                                                                      <  1:2356519/103‑1 (MQ=255)
                                                                                                                             ccAGGACCAGGTACAGTTACTCGATGCCGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAggg                      <  2:2322954/150‑1 (MQ=255)
                                                                                                                               aGGATCAGGTCCAGTTCCTCGATGCGGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGgcg                    >  2:1642082/1‑150 (MQ=255)
                                                                                                                               aGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGgcg                    >  2:1100618/1‑150 (MQ=255)
                                                                                                                                 gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGa                                                                    >  2:1782729/1‑100 (MQ=255)
                                                                                                                                 gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGa                                                                    <  1:1782729/100‑1 (MQ=255)
                                                                                                                                 gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAgg                       <  1:442841/145‑1 (MQ=255)
                                                                                                                                 gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAgg                       >  2:442841/1‑145 (MQ=255)
                                                                                                                                 gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCt                  >  2:1847666/1‑150 (MQ=255)
                                                                                                                                  cgatcttCCCGTTCTGCGCTGTAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTc                 <  2:1642402/144‑1 (MQ=255)
                                                                                                                                   tCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTcc                >  1:1476332/1‑150 (MQ=255)
                                                                                                                                    cAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCa               <  2:535307/150‑1 (MQ=255)
                                                                                                                                     aGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCAt              >  1:801679/1‑150 (MQ=255)
                                                                                                                                       gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg            <  2:2234730/150‑1 (MQ=255)
                                                                                                                                       gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg            >  1:2398654/1‑150 (MQ=255)
                                                                                                                                       gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg            <  1:1711759/150‑1 (MQ=255)
                                                                                                                                        tCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTgg           <  1:83313/150‑1 (MQ=255)
                                                                                                                                         ccAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCagatcg          >  1:1642402/1‑145 (MQ=255)
                                                                                                                                         ccAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGt          <  2:1782480/150‑1 (MQ=255)
                                                                                                                                            gTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATaga                                                                <  1:481860/93‑1 (MQ=255)
                                                                                                                                            gTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATaga                                                                >  2:481860/1‑93 (MQ=255)
                                                                                                                                             ttCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGtaata      >  1:2359176/1‑150 (MQ=255)
                                                                                                                                                 tCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACcg  >  2:1023957/1‑150 (MQ=255)
                                                                                                                                                   gATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGgc                                                                                    >  1:1690695/1‑66 (MQ=255)
                                                                                                                                                   gATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGgc                                                                                    >  1:1734756/1‑66 (MQ=255)
                                                                                                                                                   gATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGgc                                                                                    <  2:1690695/66‑1 (MQ=255)
                                                                                                                                                   gATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGgc                                                                                    <  2:1734756/66‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAGTCGCTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCG  >  NC_002947/3077005‑3077299

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