Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,011,032 T→C 6.2% L17P (CTG→CCG)  bcsG → endoglucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,011,0320TC 6.2% 225.8 / 2.7 75L17P (CTG→CCG) bcsGendoglucanase
Reads supporting (aligned to +/- strand):  ref base T (32/38);  new base C (2/3);  total (34/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 = 9.98e-01

CTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGCAACTG  >  NC_002947/3010883‑3011177
                                                                                                                                                     |                                                                                                                                                 
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                                                                                                                     >  2:1459902/1‑150 (MQ=255)
                              cGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCaaa                                                                                                                     <  1:1459902/150‑1 (MQ=255)
                                  gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGaa                                                                                                                                       <  1:1519312/128‑1 (MQ=255)
                                  gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGaa                                                                                                                                       >  2:1519312/1‑128 (MQ=255)
                                  gCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGACCAGCTCCAAAGGAGGAACCCGAGACTGCTTCATGCCGCGAACAGCCCCGCTGTACTTGCCCGCGTGCCGTGGCTGGGGGGGTGGAACCTCTATTTTTTTGCAAAGttt                                                                                                                 >  2:1189426/1‑150 (MQ=255)
                                   ccTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                                 >  2:1235375/1‑117 (MQ=255)
                                   ccTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTgg                                                                                                                                                 <  1:1235375/117‑1 (MQ=255)
                                     ttCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAgg                                                                                                                                             <  1:469936/119‑1 (MQ=255)
                                     ttCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAgg                                                                                                                                             >  2:469936/1‑119 (MQ=255)
                                       cccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCACCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTgatgg                                                                                                                                      >  1:1188372/1‑120 (MQ=18)
                                       cccccccGTTACCTGAAACCGATGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCACCACCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTgatgg                                                                                                                                      <  2:1188372/124‑5 (MQ=18)
                                        ccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgctg                                                                                                           >  2:702399/1‑150 (MQ=255)
                                         cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctct                                                                                                                                  <  2:787619/126‑1 (MQ=255)
                                         cccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACctct                                                                                                                                  >  1:787619/1‑126 (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                                                                                                                              <  2:951699/105‑1 (MQ=255)
                                                                  gTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttt                                                                                                                              >  1:951699/1‑105 (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                                                                                                                          >  1:1468876/1‑90 (MQ=255)
                                                                                     gCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttttt                                                                                                                          <  2:1468876/90‑1 (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                                    <  2:422668/150‑1 (MQ=255)
                                                                                                               caccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCttt                                    >  1:422668/1‑150 (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:1111540/150‑1 (MQ=255)
                                                                                                                              ggTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCt                     <  1:598857/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)
                                                                                                                                                     |                                                                                                                                                 
CTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGGCTGCGTATTGCCCGGCAACTG  >  NC_002947/3010883‑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: