Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,696,460 G→A 6.0% A220A (GCG→GCA)  mltF → membrane‑bound lytic transglycosylase F, murein hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,696,4600GA 6.0% 208.8 / 2.5 66A220A (GCG→GCA) mltFmembrane‑bound lytic transglycosylase F, murein hydrolase
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

CCCGGAATTAAGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATACCTGGGGCATGG  >  NC_000913/2696316‑2696606
                                                                                                                                                |                                                                                                                                                  
cccGGAATTAAGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTc                                                                                                                                                 <  6:182941/148‑1 (MQ=255)
 ccGGAATTAAGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGCTTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGa                                                                                                                                               <  2:214615/149‑1 (MQ=255)
 ccGGAATTAAGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGa                                                                                                                                               <  8:318820/149‑1 (MQ=255)
      aTTAAGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCa                                                                                                                                          <  5:196381/149‑1 (MQ=255)
         aaGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTCTTTCAGCGCGTTCACCCCGAGCTCGCCGTAGCGCTCGATATCAcc                                                                                                                                        >  1:33590/1‑148 (MQ=255)
                  ggTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGaacaac                                                                                                                              <  5:272090/149‑1 (MQ=255)
                    tAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTTCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCg                                                                                                                            >  2:66685/1‑149 (MQ=255)
                      gacgacAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGt                                                                                                                          <  8:270466/149‑1 (MQ=255)
                        cgacAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGa                                                                                                                        <  7:162009/149‑1 (MQ=255)
                         gacAAAAAATTCTCTGAGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCAGCAGCCTGTTTCAGCGCGTTCACACGGAGCTCGCCGTAGCGCTCGATATAAACGATGAACAACCGGTGAc                                                                                                                       >  3:245432/1‑149 (MQ=255)
                         gacAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCCGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGaa                                                                                                                       >  1:85063/1‑148 (MQ=255)
                              aaaaGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACACCCGGTGACTTGGt                                                                                                                  >  1:116237/1‑149 (MQ=255)
                                aaGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGttt                                                                                                                <  7:272354/149‑1 (MQ=255)
                                   gCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGc                                                                                                             >  4:119768/1‑149 (MQ=255)
                                   gCTCTGCGGAATTAATGGAAGATGTCATCGAAAGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGc                                                                                                             <  3:104373/149‑1 (MQ=255)
                                     tctGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGccc                                                                                                           >  3:127165/1‑149 (MQ=255)
                                     tctGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGccc                                                                                                           >  5:325593/1‑149 (MQ=255)
                                        gCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGt                                                                                                         >  8:312874/1‑148 (MQ=255)
                                              ttAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGAt                                                                                                                                     <  3:279744/114‑1 (MQ=255)
                                              ttAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGAt                                                                                                                                     >  4:279744/1‑114 (MQ=255)
                                                 aTGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGa                                                                                               <  8:224193/149‑1 (MQ=255)
                                                 aTGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTAGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCGCGCCGTAGCGCTCGATATCACCGATGAACATCCGGTGACTTGGTTTAGCCCGTTAGATGGCGa                                                                                               <  6:65690/149‑1 (MQ=255)
                                                    gAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGataa                                                                                            >  2:315107/1‑149 (MQ=255)
                                                    gAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGataa                                                                                            >  1:224326/1‑149 (MQ=255)
                                                       gATGTCATCGAAGGAAACCTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATAc                                                                                         >  5:108246/1‑149 (MQ=255)
                                                             aTCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCATCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCCACGATGAACACCCGGTGACTTGGTTGAGCCCGTTAGATGTCGATAATACGCTTTc                                                                                   >  4:187618/1‑149 (MQ=255)
                                                             aTCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTc                                                                                   >  8:52115/1‑149 (MQ=255)
                                                                  aGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTccgccg                                                                              >  8:221844/1‑149 (MQ=255)
                                                                      aaaCTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCt                                                                          <  5:126550/149‑1 (MQ=255)
                                                                        aCTCGAATACACCATTGCTGATTCTGTCGCCATCAGCCAGTATCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGAAATCACCGATGAACAACCGGTGACTTGGGTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGc                                                                        <  8:245594/149‑1 (MQ=255)
                                                                            gcttACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTc                                                                      <  5:314221/145‑1 (MQ=255)
                                                                            gATTACCCCATTGCTGATTCTGTCGACATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGGGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTc                                                                      <  5:207013/147‑1 (MQ=255)
                                                                                 cacCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCTTAGCGCTCGATATCACCGATGAACTACCGGTGACTTGGTTTAGCCCGTTCGATGGCGATAATACCCTTTCCGCCGCCCTTCTCGActtc                                                                <  4:62657/148‑1 (MQ=255)
                                                                                 cacCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGActtct                                                               >  3:276230/1‑149 (MQ=255)
                                                                                 cacCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGActtc                                                                <  8:88684/148‑1 (MQ=255)
                                                                                    cATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACttcttc                                                             >  6:119398/1‑148 (MQ=255)
                                                                                      ttGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCACTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATAc                                                                                         <  7:93689/118‑1 (MQ=255)
                                                                                      ttGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCACTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATAc                                                                                         <  7:98281/118‑1 (MQ=255)
                                                                                      ttGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCACTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATAc                                                                                         >  8:93689/1‑118 (MQ=255)
                                                                                      ttGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCACTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATAc                                                                                         >  8:98281/1‑118 (MQ=255)
                                                                                           gATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGActtctt                                                              <  4:61898/140‑1 (MQ=255)
                                                                                           gATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGActtctt                                                              >  3:61898/1‑140 (MQ=255)
                                                                                           gATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaat                                                     >  4:155523/1‑149 (MQ=255)
                                                                                            aTTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCGGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaat                                                     <  7:221844/148‑1 (MQ=255)
                                                                                               cTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCTTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatgaa                                                  <  3:119768/148‑1 (MQ=255)
                                                                                                tGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGGCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatgaatg                                                <  2:33590/149‑1 (MQ=255)
                                                                                                  tCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatgaat                                                 <  1:29582/146‑1 (MQ=255)
                                                                                                    gCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatgaa                                                  >  4:136946/1‑143 (MQ=255)
                                                                                                    gCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatgaa                                                  <  3:136946/143‑1 (MQ=255)
                                                                                                    gCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGa                                            >  5:217992/1‑149 (MQ=255)
                                                                                                    gCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGa                                            <  1:66685/149‑1 (MQ=255)
                                                                                                     ccATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGAc                                           >  1:253464/1‑149 (MQ=255)
                                                                                                      cATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatga                                                   <  7:85451/140‑1 (MQ=255)
                                                                                                      cATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAaatga                                                   >  8:85451/1‑140 (MQ=255)
                                                                                                            ccTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGc                                    <  2:253464/149‑1 (MQ=255)
                                                                                                                      cgcgTTCACACGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGCTACGCTGGCACGCAt                          >  2:233533/1‑149 (MQ=255)
                                                                                                                       gcgTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCAtt                         <  4:127165/149‑1 (MQ=255)
                                                                                                                          ttCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTg                        >  8:90486/1‑147 (MQ=255)
                                                                                                                           tCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTg                        >  1:39982/1‑146 (MQ=255)
                                                                                                                             aCCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATTAATGAAGACGGTACGCTGGCACGCATTGaa                      >  8:321353/1‑146 (MQ=255)
                                                                                                                                cGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGaaa                <  7:312874/149‑1 (MQ=255)
                                                                                                                                cGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGaaa                <  4:276230/149‑1 (MQ=255)
                                                                                                                                   aGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATAc             <  7:52115/149‑1 (MQ=255)
                                                                                                                                        gCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATACCTggg        <  8:59658/149‑1 (MQ=255)
                                                                                                                                        gCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATACCTg          <  5:119398/147‑1 (MQ=255)
                                                                                                                                           gTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATACCTGGGGCa     >  5:74622/1‑149 (MQ=255)
                                                                                                                                              gcgcTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATACCTGGGGCATgg  <  2:39982/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCCGGAATTAAGCTGGAAGGTAGACGACAAAAAAGGCTCTGCGGAATTAATGGAAGATGTCATCGAAGGAAAACTCGATTACACCATTGCTGATTCTGTCGCCATCAGCCTGTTTCAGCGCGTTCACCCGGAGCTCGCCGTAGCGCTCGATATCACCGATGAACAACCGGTGACTTGGTTTAGCCCGTTAGATGGCGATAATACCCTTTCCGCCGCCCTGCTCGACTTCTTCAACGAAATGAATGAAGACGGTACGCTGGCACGCATTGAAGAGAAATACCTGGGGCATGG  >  NC_000913/2696316‑2696606

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