Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,943,054 G→A 5.4% 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.4% 222.3 / 2.6 74L295F (CTC→TTC) PP_3478secretion protein
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base A (2/2);  total (36/38)
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.70e-01

CACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCT  >  NC_002947/3942905‑3943203
                                                                                                                                                     |                                                                                                                                                     
cACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAg                                                                                                                                                       >  1:638894/1‑150 (MQ=255)
   cTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCAt                                                                                                                                                    >  1:320079/1‑150 (MQ=255)
    tCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATg                                                                                                                                                   >  1:5855/1‑150 (MQ=255)
       aTCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTc                                                                                                                                                <  2:320079/150‑1 (MQ=255)
        tCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCt                                                                                                                                               >  1:1279421/1‑150 (MQ=255)
            cTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAc                                                                                                                                           >  1:1301281/1‑150 (MQ=255)
                  ttCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                                                                                     >  2:1167709/1‑150 (MQ=255)
                  ttCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                                                                                     <  2:1304865/150‑1 (MQ=255)
                  ttCCAGCACCACATCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                                                                                     >  2:205824/1‑150 (MQ=255)
                              cTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                                                                                       <  1:411487/136‑1 (MQ=255)
                              cTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                                                                                       >  2:411487/1‑136 (MQ=255)
                              cTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGtt                                                                                                                         <  2:343578/150‑1 (MQ=255)
                                cGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTgg                                                                                                                       <  2:1301281/150‑1 (MQ=255)
                                      tGGACTGGTCTTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGCCAGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCGCGTCAACGAACACGTCGCGGGTCTTGAGCATGTTCTTGACCATGTTCAACGCGGTGTTGGCGTCGa                                                                                                                 <  2:144099/150‑1 (MQ=255)
                                         aCTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGCTGc                                                                                                              <  1:205824/150‑1 (MQ=255)
                                         aCTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGc                                                                                                              <  1:1167709/150‑1 (MQ=255)
                                             gTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGt                                                                                                          <  2:108193/150‑1 (MQ=255)
                                                cTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                                                                                   <  1:915566/138‑1 (MQ=255)
                                                cTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                                                                                   >  2:915566/1‑138 (MQ=255)
                                                cTGCGACTGCAGCAGCTTCTAGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                                                                                       >  1:1238319/1‑150 (MQ=255)
                                                 tGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAgg                                                                                                      >  2:736281/1‑150 (MQ=255)
                                                 tGCGACTGCAGCAGCTTCTCGGCCAGGCGCACTGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAgg                                                                                                      <  2:685374/150‑1 (MQ=255)
                                                      cTGCAGCAGCTTCTAGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGaa                                                                                                 <  2:1238319/150‑1 (MQ=255)
                                                         cagcagCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCgg                                                                                                                             <  2:932799/119‑1 (MQ=255)
                                                         cagcagCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCgg                                                                                                                             >  1:932799/1‑119 (MQ=255)
                                                           gcagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgc                                                                                            >  1:1416535/1‑150 (MQ=255)
                                                           gcagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgc                                                                                            <  2:1416535/150‑1 (MQ=255)
                                                            cagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgcg                                                                                           >  1:338680/1‑150 (MQ=255)
                                                                  ctcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGTGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                                                                                     <  2:338680/150‑1 (MQ=255)
                                                                    cGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCa                                                                                   <  1:1352279/150‑1 (MQ=255)
                                                                       ccATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTAAGGTAGAACGTGCGCATGACCAGtt                                                                                >  2:787439/1‑150 (MQ=255)
                                                                         aTGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGaa                                                                                                 >  1:165334/1‑131 (MQ=255)
                                                                         aTGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGaa                                                                                                 <  2:165334/131‑1 (MQ=255)
                                                                            cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTgg                                                                           <  2:280734/150‑1 (MQ=255)
                                                                             gcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                            >  1:740790/1‑116 (MQ=255)
                                                                             gcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                            <  2:740790/116‑1 (MQ=255)
                                                                             gcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGt                                                                          >  1:1031843/1‑150 (MQ=255)
                                                                                 cgcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTc                                                                      >  1:293957/1‑150 (MQ=255)
                                                                                  gcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAc                                                                                                                                           <  2:145142/80‑1 (MQ=255)
                                                                                  gcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAc                                                                                                                                           >  1:145142/1‑80 (MQ=255)
                                                                                       cGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAGCGTGCGCATGACCAGTTCCTGGTAGTCCTTGGt                                                                >  2:850176/1‑150 (MQ=255)
                                                                                        gACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTcaa                                                                                                                                              >  1:1341722/1‑69 (MQ=255)
                                                                                        gACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTcaa                                                                                                                                              <  2:1341722/71‑3 (MQ=255)
                                                                                        gACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTAAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGtc                                                               <  1:787439/150‑1 (MQ=255)
                                                                                         aCGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGc                                                                                                              <  2:783084/102‑1 (MQ=255)
                                                                                         aCGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGc                                                                                                              >  1:783084/1‑102 (MQ=255)
                                                                                                gcgcATGGTCAGGGTGTTGAGGCGCTCGTCGACGAACACGTCGCGGGTCTTGAGCATGTTTTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTgcgg                                                       <  2:901904/150‑1 (MQ=255)
                                                                                                gcgcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGtctt                                                             >  2:87044/1‑144 (MQ=255)
                                                                                                gcgcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGtctt                                                             <  1:87044/144‑1 (MQ=255)
                                                                                                gcgcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTgcgg                                                       <  2:1031843/150‑1 (MQ=255)
                                                                                                        tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAACATGGTCAAGACCATGCTCa                                                                                                                                   >  2:796613/1‑66 (MQ=25)
                                                                                                        tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAACATGGTCAAGACCATGCTCa                                                                                                                                   <  1:796613/66‑1 (MQ=25)
                                                                                                           gggTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACATGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggggg                                            >  1:772870/1‑149 (MQ=255)
                                                                                                                ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTg                                                                  <  2:420346/123‑1 (MQ=255)
                                                                                                                ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTg                                                                  >  1:420346/1‑123 (MQ=255)
                                                                                                                ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGAcc                                       >  1:582324/1‑150 (MQ=255)
                                                                                                                  gAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCAt                                     >  2:645896/1‑150 (MQ=255)
                                                                                                                    ggCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCGGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGa                                   >  1:385537/1‑150 (MQ=255)
                                                                                                                           gTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGtgtt                            <  1:65626/150‑1 (MQ=255)
                                                                                                                            tCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGtgttg                           <  1:669806/150‑1 (MQ=255)
                                                                                                                              cgtgcccacTTCTCGGGTGTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTc                         >  2:1421519/7‑150 (MQ=255)
                                                                                                                              aaCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTc                         >  1:1075697/1‑150 (MQ=255)
                                                                                                                                   acacGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACTCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTcc                                                                              <  1:868912/92‑1 (MQ=255)
                                                                                                                                   acacGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACTCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTcc                                                                              >  2:868912/1‑92 (MQ=255)
                                                                                                                                    cacGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGAc                   <  2:293957/150‑1 (MQ=255)
                                                                                                                                       gTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                            <  2:635568/58‑1 (MQ=255)
                                                                                                                                       gTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                            >  1:635568/1‑58 (MQ=255)
                                                                                                                                       gTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAt                <  2:829072/150‑1 (MQ=255)
                                                                                                                                          gcgGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGa                    <  1:1368341/143‑1 (MQ=255)
                                                                                                                                          gcgGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGa                    >  2:1368341/1‑143 (MQ=255)
                                                                                                                                           cgGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtcttc            <  2:5855/150‑1 (MQ=255)
                                                                                                                                                gttGAACATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGa                    >  2:611524/2‑137 (MQ=255)
                                                                                                                                                gttGAACATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGa                    <  1:611524/136‑1 (MQ=255)
                                                                                                                                                 ttGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGgcgc      >  2:1364313/1‑150 (MQ=255)
                                                                                                                                                    aGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGc   <  2:138580/150‑1 (MQ=255)
                                                                                                                                                     gCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCt  >  1:48651/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCT  >  NC_002947/3942905‑3943203

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