Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,788,441 G→C 7.1% Y79* (TAC→TAG)  PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,788,4410GC 7.1% 318.0 / 13.7 99Y79* (TAC→TAG) PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (49/43);  new base C (4/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCC  >  NC_002947/3788300‑3788586
                                                                                                                                             |                                                                                                                                                 
ggTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGc                                                                                                                                           <  1:1620829/150‑1 (MQ=255)
ggTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGc                                                                                                                                           <  2:594996/150‑1 (MQ=255)
 gTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGcc                                                                                                                                          >  2:1058686/1‑150 (MQ=255)
         gACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGcccc                                                                                                                                  >  2:1577817/1‑150 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGTGCCTCAGCCCCGAg                                                                                                                               <  2:1478983/150‑1 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAg                                                                                                                               >  2:26584/1‑150 (MQ=255)
             tGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCGGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGt                                                                                                                              >  2:162953/1‑150 (MQ=255)
              ggACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGtt                                                                                                                             >  1:1259096/1‑150 (MQ=255)
              ggACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGtt                                                                                                                             <  1:114591/150‑1 (MQ=255)
                  gAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAg                                                                                                                         >  2:1220811/1‑150 (MQ=255)
                          gCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCg                                                                                                                 >  2:21466/1‑150 (MQ=255)
                           cACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGa                                                                                                                <  1:344066/150‑1 (MQ=255)
                           cACGGCAGCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGa                                                                                                                <  2:1417741/150‑1 (MQ=255)
                            aCGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGccc                                                                                                                                   <  2:226630/130‑1 (MQ=255)
                            aCGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGccc                                                                                                                                   >  1:226630/1‑130 (MQ=255)
                              ggCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGACGGGAGACAATGAGTTTGTCCGAAGCATCGGCCAGGTCGCCCGCCGCCGCCCCGAGCCGCAGCCAGGCCGActc                                                                                                             >  1:1492472/1‑150 (MQ=255)
                                 aTCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGa                                                                                                          <  1:1577817/150‑1 (MQ=255)
                                   cgtcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACa                                                                                                        <  1:1389890/150‑1 (MQ=255)
                                   cgtcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACa                                                                                                        <  1:251456/150‑1 (MQ=255)
                                    gtcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAg                                                                                                       >  1:1483802/1‑150 (MQ=255)
                                       gAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCg                                                                                                    >  1:736256/1‑150 (MQ=255)
                                             gTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTa                                                                                              >  1:657374/1‑150 (MQ=255)
                                               gCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCa                                                                                                                          >  2:1706821/1‑120 (MQ=255)
                                               gCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCa                                                                                                                          <  1:1706821/120‑1 (MQ=255)
                                               gCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAgc                                                                                            <  1:1668598/150‑1 (MQ=255)
                                                           gggCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGa                                                                                >  2:1882335/1‑150 (MQ=255)
                                                           gggCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGa                                                                                >  2:1882811/1‑150 (MQ=255)
                                                               cTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGCCTCTGACCGACGACGGTAGGGCAGGCCGTCGaaacg                                                                            >  2:7943/1‑148 (MQ=255)
                                                               cTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAgg                                                                            <  2:1492472/150‑1 (MQ=255)
                                                                tGCAGCCGTTGCTGCAGATAGGAATCAACGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAggg                                                                           >  1:50451/1‑150 (MQ=255)
                                                                     ccGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATc                                                                                  >  2:454160/1‑138 (MQ=255)
                                                                     ccGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATc                                                                                  <  1:454160/138‑1 (MQ=255)
                                                                     ccGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGt                                                                      >  1:1484770/1‑150 (MQ=255)
                                                                      cGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                     <  1:1882811/150‑1 (MQ=255)
                                                                      cGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                     <  1:1882335/150‑1 (MQ=255)
                                                                              cAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAgcgc                                                                                          <  2:1708805/121‑1 (MQ=255)
                                                                              cAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAgcgc                                                                                          >  1:1708805/1‑121 (MQ=255)
                                                                               gatctggAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCCGCCTCGACTTGAAGATAGACCGCCCCTCACATCCGAGGCTAGCGCAGTCCATCGAAAGGGGCAGTGGCCTTGGCg                                                            <  1:529063/145‑1 (MQ=255)
                                                                                 agaGGAATCGTCGTTCAGCACGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCggg                                                          <  2:1756965/148‑1 (MQ=255)
                                                                                 aTAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCATTGGCgcg                                                          <  2:1484770/150‑3 (MQ=255)
                                                                                 aTAGGAATCGTCGTTAAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCggg                                                          <  1:1647899/150‑1 (MQ=255)
                                                                                    ggAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtc                                                       <  2:737734/150‑1 (MQ=255)
                                                                                    ggAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtc                                                       >  2:1860934/1‑150 (MQ=255)
                                                                                     gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCggagatc                                                      >  2:529063/1‑145 (MQ=255)
                                                                                     gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGcc                                                                                                                                          <  2:89040/66‑1 (MQ=255)
                                                                                     gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGcc                                                                                                                                          >  1:89040/1‑66 (MQ=255)
                                                                                      aaTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTa                                                                                              >  1:21492/1‑109 (MQ=255)
                                                                                      aaTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTa                                                                                              <  2:21492/109‑1 (MQ=255)
                                                                                      aaTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgat                                                                                                                                       >  1:1672610/1‑65 (MQ=255)
                                                                                      aaTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgat                                                                                                                                       <  2:1672610/68‑4 (MQ=38)
                                                                                       aTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                    >  1:541155/1‑150 (MQ=255)
                                                                                       aTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                    <  2:541155/150‑1 (MQ=255)
                                                                                         cgtcgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCga                                                  <  1:26584/150‑1 (MQ=255)
                                                                                         cgtcgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCga                                                  <  2:1587882/150‑1 (MQ=255)
                                                                                           tcgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTa                                                                                              <  2:636305/104‑1 (MQ=255)
                                                                                           tcgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTa                                                                                              >  1:636305/1‑104 (MQ=255)
                                                                                           tcgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCgaga                                                >  2:1395164/1‑150 (MQ=255)
                                                                                                  gCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGt                                         <  2:736256/150‑1 (MQ=255)
                                                                                                  gCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGACCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGt                                         >  1:325052/1‑150 (MQ=255)
                                                                                                   cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTa                                        >  1:1765725/1‑150 (MQ=255)
                                                                                                    ccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAgg                                                                            >  2:232661/1‑113 (MQ=255)
                                                                                                    ccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAgg                                                                            <  1:232661/113‑1 (MQ=255)
                                                                                                      gcgcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAACGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTc                                     <  1:1058686/150‑1 (MQ=255)
                                                                                                        gcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCt                                   <  1:21466/150‑1 (MQ=255)
                                                                                                        gcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCt                                   <  2:657374/150‑1 (MQ=255)
                                                                                                        gcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCt                                   <  1:1220811/150‑1 (MQ=255)
                                                                                                        gcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCt                                   <  2:1765725/150‑1 (MQ=255)
                                                                                                          aGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                    <  2:940414/131‑1 (MQ=255)
                                                                                                          aGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                    >  1:940414/1‑131 (MQ=255)
                                                                                                            gCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGAt                               >  2:377135/1‑150 (MQ=255)
                                                                                                             ccAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATc                              >  1:977554/1‑150 (MQ=255)
                                                                                                              cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATcc                             >  2:1568227/1‑150 (MQ=255)
                                                                                                              cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATcc                             >  1:882550/1‑150 (MQ=255)
                                                                                                              cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATcc                             >  2:300471/1‑150 (MQ=255)
                                                                                                                  gaagGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATCGTAGCGCAGGCCCTAGAAAGGGGCAGTGCCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCg                         >  2:490949/1‑150 (MQ=255)
                                                                                                                       gTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGACATCGAAAGGGGCAGTGGCCTTGGCggg                                                          <  1:1689033/112‑1 (MQ=255)
                                                                                                                       gTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGACATCGAAAGGGGCAGTGGCCTTGGCggg                                                          >  2:1689033/1‑112 (MQ=255)
                                                                                                                        tttGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                          >  2:173551/1‑127 (MQ=255)
                                                                                                                        tttGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                          <  1:173551/127‑1 (MQ=255)
                                                                                                                          tGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGccc                 <  2:882550/150‑1 (MQ=255)
                                                                                                                           gTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGcccc                <  1:1568227/150‑1 (MQ=255)
                                                                                                                            tCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCt               >  1:1248608/1‑150 (MQ=255)
                                                                                                                            tCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCt               >  2:83937/1‑150 (MQ=255)
                                                                                                                            tCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCt               >  2:678370/1‑150 (MQ=255)
                                                                                                                             ccGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTg              <  1:83937/150‑1 (MQ=255)
                                                                                                                              cGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTgg             >  1:1111385/1‑150 (MQ=255)
                                                                                                                               gAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGc            >  2:155011/1‑150 (MQ=255)
                                                                                                                               gAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGc            <  2:1250888/150‑1 (MQ=255)
                                                                                                                               gAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGCCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGc            >  1:1250888/1‑150 (MQ=255)
                                                                                                                                 ggctgaggCCCGCTACCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCt          >  2:540349/7‑150 (MQ=255)
                                                                                                                                   cATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                     >  1:1769975/1‑89 (MQ=255)
                                                                                                                                   cATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                     <  2:1769975/89‑1 (MQ=255)
                                                                                                                                   cATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCAATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc        >  1:1549825/1‑150 (MQ=255)
                                                                                                                                     tCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcccc      >  2:452054/1‑150 (MQ=255)
                                                                                                                                     gaggCCCGCTACCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCGcccc      >  2:242913/3‑150 (MQ=255)
                                                                                                                                       ggCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccgg    >  2:423323/1‑150 (MQ=255)
                                                                                                                                        gCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGaaa                                                                              >  2:663944/1‑75 (MQ=255)
                                                                                                                                        gCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGaaa                                                                              <  1:663944/75‑1 (MQ=255)
                                                                                                                                         cccGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccgacc  <  1:155011/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCC  >  NC_002947/3788300‑3788586

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