Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,382,545 A→C 5.6% D254E (GAT→GAG)  PP_1204 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,382,5450AC 5.6% 199.8 / 3.9 74D254E (GAT→GAG) PP_1204lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base C (2/2);  total (35/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 = 4.86e-01

AGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAACGGGCT  >  NC_002947/1382402‑1382694
                                                                                                                                               |                                                                                                                                                       
aGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                                      <  2:28868/147‑1 (MQ=255)
aGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                                      >  1:28868/1‑147 (MQ=255)
aGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAg                                                                                                                                                   >  2:1182597/1‑150 (MQ=255)
 ggAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGa                                                                                                                                                  >  1:1325259/1‑150 (MQ=255)
    aCTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGTCAGTCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGATGACCTCGAAGAtga                                                                                                                                               <  2:1412294/150‑3 (MQ=255)
       ccTGCTGCAGGTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCt                                                                                                                                            <  2:299468/150‑1 (MQ=255)
              cAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGc                                                                                                                                    >  2:779146/1‑150 (MQ=255)
                 cTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTc                                                                                                                                 <  1:1138364/150‑1 (MQ=255)
                   ccGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGa                                                                                                                               <  1:1036263/150‑1 (MQ=255)
                   ccGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGa                                                                                                                               <  2:1772489/150‑1 (MQ=255)
                    cGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTTCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGag                                                                                                                              >  1:151827/1‑150 (MQ=255)
                    cGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGag                                                                                                                              >  2:972771/1‑150 (MQ=255)
                      gagcGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGagca                                                                                                                            >  1:1084590/4‑150 (MQ=255)
                           ttGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑Ta                                                                                                                      >  1:581563/1‑150 (MQ=255)
                                      cTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGGGGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAgggg                                                                                                           >  1:1450855/1‑149 (MQ=255)
                                       tCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATAGAAAGTGTGTGTATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTg                                                                                                          >  2:1096500/1‑150 (MQ=255)
                                       tCGATCAGGCTCTGTTCCATTTCACTGGCTACCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCAACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTCCGCGGGGTg                                                                                                          >  2:244926/1‑150 (MQ=255)
                                         gATCAGGCTCTGTTCCATTTCACTGGCTTCCTATTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGgc                                                                                                        >  1:1400000/1‑150 (MQ=255)
                                              ggCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACGTCATCATCATCGAGGCCGACAAAAATCACCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGa                                                                                                   <  1:1470653/150‑1 (MQ=255)
                                                    gTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCCAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACa                                                                                             <  1:1096500/150‑1 (MQ=255)
                                                     ttCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAg                                                                                            >  2:1332953/1‑150 (MQ=255)
                                                        cATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATg                                                                                                                                        >  2:1871314/1‑105 (MQ=255)
                                                        cATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATg                                                                                                                                        <  1:1871314/105‑1 (MQ=255)
                                                        cATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCt                                                                                         <  2:1103111/150‑1 (MQ=255)
                                                         aTTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                                        <  1:943859/150‑1 (MQ=255)
                                                             cACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCAt                                                                                    <  1:1907684/150‑1 (MQ=255)
                                                              aCTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATc                                                                                   <  2:581563/150‑1 (MQ=255)
                                                                tGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑CGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGCCTGGCATCGc                                                                                 <  2:1325259/150‑1 (MQ=255)
                                                                 ggCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGcc                                                                                <  2:1920857/150‑1 (MQ=255)
                                                                 ggCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTCGCCCCGGTGCACACACTGGCCTACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGcc                                                                                <  2:151827/150‑1 (MQ=255)
                                                                   cTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGNCTc                                                                              <  2:1200610/150‑1 (MQ=255)
                                                                    ttCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑AATTGCGCAGGGAGGACCTGAACAACAGGCTAGCATCGACtag                                                                             >  2:781663/1‑148 (MQ=255)
                                                                     tCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGt                                                                            >  2:280676/1‑150 (MQ=255)
                                                                      ccTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTa                                                                           >  1:981298/1‑150 (MQ=255)
                                                                        tctTCAAGGCCCAGTTCGGAGAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTAcc                                                                         >  1:965039/1‑150 (MQ=255)
                                                                        tctTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTAcc                                                                         >  2:999249/1‑150 (MQ=255)
                                                                        tctTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTAcc                                                                         >  1:939750/1‑150 (MQ=255)
                                                                          ttCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGagcag                                                                                                                           >  2:810171/1‑99 (MQ=255)
                                                                          ttCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGagcag                                                                                                                           <  1:810171/99‑1 (MQ=255)
                                                                             aaGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCAtt                                                                    <  1:972771/150‑1 (MQ=255)
                                                                             aaGGACCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCAtt                                                                    <  2:608324/150‑1 (MQ=255)
                                                                             aaGCCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCTAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCAtt                                                                    <  1:7522/150‑1 (MQ=255)
                                                                              aGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTc                                                                                                                                 <  2:314415/89‑1 (MQ=255)
                                                                              aGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTc                                                                                                                                 >  1:314415/1‑89 (MQ=255)
                                                                               ggCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTCGCATCGCCTCGTACCGCATTGt                                                                  <  2:1226962/150‑1 (MQ=255)
                                                                                gCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTc                                                                 >  2:153361/1‑150 (MQ=255)
                                                                                   cAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTc                                                              <  1:453179/150‑1 (MQ=255)
                                                                                      ttCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTc                                                           >  1:655320/1‑150 (MQ=255)
                                                                                        cGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCaa                                                         <  1:1798580/150‑1 (MQ=255)
                                                                                            cAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATc                                                                                   >  2:1889810/1‑120 (MQ=255)
                                                                                            cAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATc                                                                                   <  1:1889810/120‑1 (MQ=255)
                                                                                                 gcTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGaa                                                                                                                                                    <  2:1623442/52‑1 (MQ=255)
                                                                                                 gcTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGaa                                                                                                                                                    >  1:1623442/1‑52 (MQ=255)
                                                                                                   ttGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGC‑CAGGAAGCTTCCAGCATCG‑ACCTGAGCGGGGTGGGGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcagg                                              >  1:1177701/1‑150 (MQ=255)
                                                                                                    tGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggt                                             <  2:263125/150‑1 (MQ=255)
                                                                                                     taCCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGGCATCGAAGATCAGCCCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggag                                            >  2:556144/2‑148 (MQ=255)
                                                                                                      aCCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCAtt                                                                    >  2:1509883/1‑125 (MQ=255)
                                                                                                      aCCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCAtt                                                                    <  1:1509883/125‑1 (MQ=255)
                                                                                                            tGCACACACTGGCCGACATCATCATCTTCGAGGTCGTCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTTCTTGCAGGGGCAGGGTg                                     <  2:192845/150‑1 (MQ=255)
                                                                                                                cccACCCGCCGACATCATCATCTTCGAGGTCATAGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGa                                 <  2:1400000/150‑1 (MQ=255)
                                                                                                                    cTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTg                             <  1:999249/150‑1 (MQ=255)
                                                                                                                         cGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGa                        >  2:991507/1‑150 (MQ=255)
                                                                                                                           aCATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAgg                      <  2:939750/150‑1 (MQ=255)
                                                                                                                                 tgatcTTCGATGACCTCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGc                >  1:216236/3‑150 (MQ=255)
                                                                                                                                 tcatcTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGc                >  2:1294356/1‑150 (MQ=255)
                                                                                                                                 tcatcTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGc                <  1:1294356/150‑1 (MQ=255)
                                                                                                                                  catcTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCt               >  1:107434/1‑150 (MQ=255)
                                                                                                                                   atcTTCGATGACCTCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTc                                                                 >  2:560679/1‑99 (MQ=255)
                                                                                                                                   atcTTCGATGACCTCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTc                                                                 <  1:560679/99‑1 (MQ=255)
                                                                                                                                   atcTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTc              >  1:634855/1‑150 (MQ=255)
                                                                                                                                      ttCGAGGTCATCGAAGATCAGCTCATG‑AGAATAGAG‑ACCTGTAATGCGCAGGGTGGCGCTGACCAACAGGTCCTGATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGga           <  1:888088/150‑1 (MQ=255)
                                                                                                                                           ggTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAACg      <  1:991507/150‑1 (MQ=255)
                                                                                                                                           ggTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAACg      <  1:153361/150‑1 (MQ=255)
                                                                                                                                             tCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAACggg    <  2:634855/150‑1 (MQ=255)
                                                                                                                                               aTCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAACGGGCt  >  1:193045/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
AGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATG‑GGCTTCGAGCAGCT‑TATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAACGGGCT  >  NC_002947/1382402‑1382694

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