Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,011,029 G→T 5.3% G16V (GGA→GTA)  bcsG → endoglucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,011,0290GT 5.3% 248.4 / 2.7 76G16V (GGA→GTA) bcsGendoglucanase
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base T (2/2);  total (35/41)
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.81e-01

GCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGCAACTG  >  NC_002947/3010881‑3011177
                                                                                                                                                    |                                                                                                                                                    
gcCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGa                                                                                                                                                     >  2:894761/1‑150 (MQ=255)
  cTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACt                                                                                                                                                   >  1:1801983/1‑150 (MQ=255)
   tgtgTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTg                                                                                                                                                  >  1:1561929/1‑150 (MQ=255)
    gtgtTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                                 <  2:30288/150‑1 (MQ=255)
     tgtTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTggg                                                                                                                                                >  2:1937351/1‑150 (MQ=255)
         cGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTCGCCTGGACTGGGAggg                                                                                                                                            <  1:1042009/150‑1 (MQ=255)
          ggCATTACGTGGCGTGTTGCGCCGTCGCGTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCCCCCCCCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCCCCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGt                                                                                                                                           <  1:120249/150‑1 (MQ=255)
             aTTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCAGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGa                                                                                                                                        >  1:873225/1‑150 (MQ=255)
               tACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAAc                                                                                                                                      >  2:1461827/1‑150 (MQ=255)
                aCGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAAcc                                                                                                                                     <  1:1911882/150‑1 (MQ=255)
                   tGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctct                                                                                                                                  >  2:1378024/1‑150 (MQ=255)
                           tGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttttt                                                                                                                          <  2:1801983/150‑1 (MQ=255)
                            gcgcCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTg                                                                                                                         <  2:861871/150‑1 (MQ=255)
                                cGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCaaa                                                                                                                     <  1:1459902/150‑1 (MQ=255)
                                cGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCaaa                                                                                                                     >  2:1459902/1‑150 (MQ=255)
                                    gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGaa                                                                                                                                       <  1:1519312/128‑1 (MQ=255)
                                    gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGaa                                                                                                                                       >  2:1519312/1‑128 (MQ=255)
                                    gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGACCAGCTCCAAAGGAGGAACCCGAGACTGCTTCATGCCGCGAACAGCCCCGCTGTACTTGCCCGCGTGCCGTGGCTGGGGGGGTGGAACCTCTATTTTTTTGCAAAGttt                                                                                                                 >  2:1189426/1‑150 (MQ=255)
                                     ccTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                                 <  1:1235375/117‑1 (MQ=255)
                                     ccTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                                 >  2:1235375/1‑117 (MQ=255)
                                       ttCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAgg                                                                                                                                             >  2:469936/1‑119 (MQ=255)
                                       ttCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAgg                                                                                                                                             <  1:469936/119‑1 (MQ=255)
                                         cccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCACCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTgatgg                                                                                                                                      >  1:1188372/1‑120 (MQ=18)
                                         cccccccGTTACCTGAAACCGATGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCACCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTgatgg                                                                                                                                      <  2:1188372/124‑5 (MQ=18)
                                          ccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgctg                                                                                                           >  2:702399/1‑150 (MQ=255)
                                           cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctct                                                                                                                                  >  1:787619/1‑126 (MQ=255)
                                           cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctct                                                                                                                                  <  2:787619/126‑1 (MQ=255)
                                           cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTgg                                                                                                          <  2:926077/150‑1 (MQ=255)
                                               cGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTcgc                                                                                                      <  1:702399/150‑1 (MQ=255)
                                                    ccTGAAACCGCCGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACAGCGGCGCACTCTGTACTGGCCCGCTGGCCTGGACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                                                 <  1:1304832/150‑1 (MQ=255)
                                                     cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                                >  1:577301/1‑150 (MQ=255)
                                                     cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                                                >  1:573694/1‑150 (MQ=255)
                                                      tGAAACCGCTGGGCTTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGc                                                                                               >  1:251976/1‑150 (MQ=255)
                                                            cGCTGGGCGTCCGCTGCCGCGCTGCACCCCAGGAACCCCACAATGACTAATGCCCCCATCAGCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTc                                                                                         <  2:577301/150‑1 (MQ=255)
                                                            cGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTc                                                                                         <  1:910304/150‑1 (MQ=255)
                                                             gCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCg                                                                                        >  2:1563081/1‑150 (MQ=255)
                                                              cTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAANCTCTATNTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGa                                                                                       <  1:1356143/150‑1 (MQ=255)
                                                                   cGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTcc                                                                                  <  1:412646/150‑1 (MQ=255)
                                                                    gTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttt                                                                                                                              >  1:951699/1‑105 (MQ=255)
                                                                    gTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttt                                                                                                                              <  2:951699/105‑1 (MQ=255)
                                                                     tACGCCGCCGCGCTGCACCCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAg                                                                                <  1:1042391/150‑1 (MQ=255)
                                                                           gccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTg                                                                          <  2:91903/150‑1 (MQ=255)
                                                                            ccgcGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCCCCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGc                                                                         <  1:750151/150‑1 (MQ=255)
                                                                              gcgcTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTgat                                                                                                                                        >  1:1924025/1‑83 (MQ=9)
                                                                              gcgcTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTgat                                                                                                                                        <  2:1924025/85‑3 (MQ=9)
                                                                                    cACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCtg                                                                 <  2:1561929/150‑1 (MQ=255)
                                                                                    cACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGACATGCTc                                                                                         >  1:1655024/1‑126 (MQ=255)
                                                                                    cACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGACATGCTc                                                                                         <  2:1655024/126‑1 (MQ=255)
                                                                                      cGCCAGGAACCCCAGAATGACTAATACCACCATCCCCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCCCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCtgtt                                                               <  1:1236981/150‑1 (MQ=255)
                                                                                       gCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttttt                                                                                                                          <  2:1468876/90‑1 (MQ=255)
                                                                                       gCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttttt                                                                                                                          >  1:1468876/1‑90 (MQ=255)
                                                                                            gAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGc                                                         <  2:251976/150‑1 (MQ=255)
                                                                                              aCCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGTTg                                                       <  1:1937351/150‑1 (MQ=255)
                                                                                                 ccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGt                                                    >  2:549260/1‑150 (MQ=255)
                                                                                                  cAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGtt                                                   >  2:795296/1‑150 (MQ=255)
                                                                                                         aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCtt                                            <  1:1378024/150‑1 (MQ=255)
                                                                                                           tAATACCACCATCTCCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGt                                          <  1:977071/150‑1 (MQ=255)
                                                                                                             aTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                              <  2:985551/144‑1 (MQ=255)
                                                                                                             aTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                                              >  1:985551/1‑144 (MQ=255)
                                                                                                                ccaccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCtt                                     >  2:1111540/1‑150 (MQ=255)
                                                                                                                 caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                                    >  1:422668/1‑150 (MQ=255)
                                                                                                                 caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                                    <  2:422668/150‑1 (MQ=255)
                                                                                                                  accaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTg                                   >  2:806942/1‑150 (MQ=255)
                                                                                                                     aTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCATAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCg                                >  2:1129408/1‑150 (MQ=255)
                                                                                                                      tCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCgg                               <  1:549260/150‑1 (MQ=255)
                                                                                                                               cGGTACTGGCCCGCTGGCCTGGACTGGGAGGTTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTTCGTTTCTGCTTGTCTCTTTGCCGGGCAAGTGGc                      >  2:233499/1‑150 (MQ=255)
                                                                                                                               cGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGc                      <  1:894761/150‑1 (MQ=255)
                                                                                                                                ggTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCt                     <  1:598857/150‑1 (MQ=255)
                                                                                                                                ggTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCt                     <  1:1111540/150‑1 (MQ=255)
                                                                                                                                         ccGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTTTTTTGTTTCTTCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTCCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGc            >  2:560142/1‑150 (MQ=255)
                                                                                                                                         ccGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGc            >  2:33034/1‑150 (MQ=255)
                                                                                                                                           gCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCACAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGccc          <  1:1189426/150‑1 (MQ=255)
                                                                                                                                                 cTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGc       >  2:1536118/1‑147 (MQ=255)
                                                                                                                                                 cTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGc       <  1:1536118/147‑1 (MQ=255)
                                                                                                                                                  tGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGc                                                                         <  1:1623128/80‑1 (MQ=255)
                                                                                                                                                  tGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGc                                                                         >  2:1623128/1‑80 (MQ=255)
                                                                                                                                                   ggACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGCAActg  >  1:817171/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGCAACTG  >  NC_002947/3010881‑3011177

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