Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,943,054 G→A 5.1% L295F (CTC→TTC)  PP_3478 ← secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,943,0540GA 5.1% 243.2 / 2.6 79L295F (CTC→TTC) PP_3478secretion protein
Reads supporting (aligned to +/- strand):  ref base G (38/37);  new base A (2/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

CTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAG  >  NC_002947/3942908‑3943201
                                                                                                                                                  |                                                                                                                                                   
cTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCAt                                                                                                                                                  >  1:555953/1‑150 (MQ=255)
 tCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATg                                                                                                                                                 <  2:336147/150‑1 (MQ=255)
  ccATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATgg                                                                                                                                                >  1:619195/1‑150 (MQ=255)
   cATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGt                                                                                                                                               >  2:620611/1‑150 (MQ=255)
     tCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCt                                                                                                                                             >  1:470353/1‑150 (MQ=255)
           ccACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCa                                                                                                                                       <  2:470353/150‑1 (MQ=255)
           ccACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGGGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCa                                                                                                                                       >  2:46768/1‑150 (MQ=255)
            cACTTCCAGCACCACCCCGGGGTTGGACTGGTCCTGCGACTGCTGCAGATTCTCGGTCATGCCGACCTCTCCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCAt                                                                                                                                      <  1:603711/150‑1 (MQ=255)
              cTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAACATGGTcaa                                                                                                                                            <  1:1168892/142‑3 (MQ=255)
              cTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAACATGGTcaa                                                                                                                                            >  2:1168892/1‑140 (MQ=255)
              cTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTACAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGt                                                                                                                                    <  1:1250348/150‑1 (MQ=255)
               ttCCAGCACCACCTCGGGGTTGGACTGGTGATGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCGGACATGTCCCTCATGGTCGGCGGGTTGAGGCGCGCGGCAACGAACACGTCGCGGGTCTTCAGCATGGGGTTGACGATGtt                                                                                                                                   >  2:54840/1‑150 (MQ=255)
               ttCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGGGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                                                                                   <  2:1523317/150‑1 (MQ=255)
                             cGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                                                                                     <  2:1438870/134‑1 (MQ=255)
                             cGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                                                                                     >  1:1438870/1‑134 (MQ=255)
                                 gTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                                                                                 <  2:1330954/150‑1 (MQ=255)
                                    ggACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGAt                                                                                                              <  1:991760/150‑1 (MQ=255)
                                      ccatctCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGAGGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGc                                                                                                            <  1:1413477/145‑1 (MQ=255)
                                        tGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                          <  1:1316056/150‑1 (MQ=255)
                                        tGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                          <  1:987809/150‑1 (MQ=255)
                                          gTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGt                                                                                                        >  1:568274/1‑150 (MQ=255)
                                           tCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTg                                                                                                                                           >  2:597931/1‑114 (MQ=255)
                                           tCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTg                                                                                                                                           <  1:597931/114‑1 (MQ=255)
                                              tGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAgg                                                                                                    <  1:230615/150‑1 (MQ=255)
                                               gCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGt                                                                                                   >  2:600728/1‑150 (MQ=255)
                                                       agcagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTg                                                                                           >  2:1283936/1‑150 (MQ=255)
                                                       agcagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTg                                                                                           >  1:1327581/1‑150 (MQ=255)
                                                         cagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgcg                                                                                         >  1:140847/1‑150 (MQ=255)
                                                                tcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGCGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAcc                                                                                  >  1:1117332/1‑150 (MQ=255)
                                                                    ccATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGtt                                                                              <  2:619195/150‑1 (MQ=255)
                                                                        gcgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                                                                          <  2:1117332/150‑1 (MQ=255)
                                                                         cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTgg                                                                         <  1:1283936/150‑1 (MQ=255)
                                                                          gcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGt                                                                        >  2:1116421/1‑150 (MQ=255)
                                                                           cACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTa                                                                       <  2:59443/150‑1 (MQ=255)
                                                                           cACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTa                                                                       >  2:598781/1‑150 (MQ=255)
                                                                            aCCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAg                                                                      <  2:555953/150‑1 (MQ=255)
                                                                               gcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTcc                                                                   <  1:598781/150‑1 (MQ=255)
                                                                                           tCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTgc                                                       >  2:1201814/1‑150 (MQ=255)
                                                                                               gcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTgcggcg                                                   >  2:1582186/1‑150 (MQ=255)
                                                                                                 aTGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGaa                                                 <  1:1471467/150‑1 (MQ=255)
                                                                                                    gTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCg                                              <  1:399223/150‑1 (MQ=255)
                                                                                                    gTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCg                                              <  1:600728/150‑1 (MQ=255)
                                                                                                     tCAGGGTGTTGAGGCGCTCGTCAACGAGCACGTCGTGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCgg                                             >  1:1603758/1‑150 (MQ=255)
                                                                                                     tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCgg                                             <  1:620611/150‑1 (MQ=255)
                                                                                                     tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCgg                                             >  1:388783/1‑150 (MQ=255)
                                                                                                        gggTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGt                                          >  2:218237/1‑150 (MQ=255)
                                                                                                         ggTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTa                                         >  2:1199825/1‑150 (MQ=255)
                                                                                                             ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAAc                                                                                              >  2:982979/1‑93 (MQ=255)
                                                                                                             ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAAc                                                                                              >  2:983174/1‑93 (MQ=255)
                                                                                                             ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAAc                                                                                              <  1:983174/93‑1 (MQ=255)
                                                                                                             ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAAc                                                                                              <  1:982979/93‑1 (MQ=255)
                                                                                                             ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGAcc                                     <  2:140847/150‑1 (MQ=255)
                                                                                                               gAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCAt                                   >  1:748319/1‑150 (MQ=255)
                                                                                                               gAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCAt                                   >  1:234873/1‑150 (MQ=255)
                                                                                                                aGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTACTTGGTCTTCTGCGGCGAATCGGGGTAGACCATg                                  >  1:756788/1‑150 (MQ=255)
                                                                                                                        gTCAGCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGAGGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGtgtt                          <  1:437710/150‑1 (MQ=255)
                                                                                                                        gTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGtgtt                          <  1:46768/150‑1 (MQ=255)
                                                                                                                           aaCGAGCACGTCGTGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTc                       <  2:1603758/150‑1 (MQ=255)
                                                                                                                           aaCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGt                                                                                                        <  1:210297/69‑1 (MQ=255)
                                                                                                                           aaCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGt                                                                                                        >  2:210297/1‑69 (MQ=255)
                                                                                                                           aaCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                                                                                      <  1:1473221/71‑1 (MQ=255)
                                                                                                                           aaCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                                                                                      >  2:1473221/1‑71 (MQ=255)
                                                                                                                            aCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCg                      <  2:234873/150‑1 (MQ=255)
                                                                                                                              gAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGtt                    >  1:818590/1‑150 (MQ=255)
                                                                                                                                    gTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAt              >  1:639972/1‑150 (MQ=255)
                                                                                                                                         gggTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtcttct         >  2:211411/1‑150 (MQ=255)
                                                                                                                                         gggTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtcttct         >  2:1077229/1‑150 (MQ=255)
                                                                                                                                          ggTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTg        >  2:411728/1‑150 (MQ=255)
                                                                                                                                           gTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCTGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTgg       <  1:1116421/150‑1 (MQ=255)
                                                                                                                                            tCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                            >  2:1623341/1‑112 (MQ=255)
                                                                                                                                            tCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                            <  1:1623341/112‑1 (MQ=255)
                                                                                                                                            tCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGgt                             <  1:1515066/127‑1 (MQ=255)
                                                                                                                                            tCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGgt                             >  2:1515066/1‑127 (MQ=255)
                                                                                                                                             gttGAACATGGTCAAGACCATGCTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTa                                                                       <  1:1035108/83‑1 (MQ=37)
                                                                                                                                             gttGAACATGGTCAAGACCATGCTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTa                                                                       >  2:1035108/2‑84 (MQ=37)
                                                                                                                                              ttGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGAc                                      >  2:336992/1‑116 (MQ=255)
                                                                                                                                              ttGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGAc                                      <  1:336992/116‑1 (MQ=255)
                                                                                                                                               tGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCa   <  2:818590/150‑1 (MQ=255)
                                                                                                                                                gAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAg  >  1:375277/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAG  >  NC_002947/3942908‑3943201

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: ‑