Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,808,042 G→A 10.4% intergenic (+121/+34) mrcA → / ← maeB penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,808,0420GA10.4% 168.1 / 7.6 68intergenic (+121/+34)mrcA/maeBpenicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base A (3/4);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGATCA  >  NC_002947/5807903‑5808191
                                                                                                                                            |                                                                                                                                                      
gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCgggg                                                                                                                                              <  2:201704/150‑1 (MQ=255)
 ccGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGc                                                                                                                                             <  2:965433/150‑1 (MQ=255)
    aTGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGGTTTGCGTGGGCGCTCGATCTAATAGGCGCCACAAGGGCCCTGACCAACACCTTGCAGCCCTGATACATTCCCAAGCCCCATAAAAAAGCCCCGGCGCATGCGCACCGGGGCttt                                                                                                                                          >  2:1002887/1‑150 (MQ=255)
           tGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACAGCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATgg                                                                                                                                   <  1:609763/150‑2 (MQ=38)
            gTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCACATAAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATggg                                                                                                                                  <  1:111571/150‑3 (MQ=38)
             ttCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCg                                                                                                                                 >  1:195308/1‑150 (MQ=255)
             ttCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCg                                                                                                                                 >  1:1443715/1‑150 (MQ=255)
                      aGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATgg                                                                                                                                   >  1:261731/1‑138 (MQ=38)
                      aGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATgg                                                                                                                                   <  2:261731/139‑2 (MQ=38)
                         gcgcCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATc                                                                                                                     <  2:1538697/150‑1 (MQ=255)
                           gcCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCGTGGCACTGTCCCAAGCCCCATAACACCAGTCACGATAATGGAAGTCGTGGGTTTTTTATGTCGCTTGCTGCAATCAg                                                                                                                   >  1:1867900/1‑150 (MQ=255)
                           gcCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAg                                                                                                                   >  1:1725665/1‑150 (MQ=255)
                                    cggTTGTGCGTGGGCGCTCACACAAACAGGCGCCACAAGCGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAac                                                                                                         <  2:441791/149‑1 (MQ=255)
                                                     tCGATCTGCTAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGt                                                                                        >  2:1168168/1‑150 (MQ=255)
                                                     tCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCAc                                                                                              >  1:1451084/1‑144 (MQ=255)
                                                     tCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCAc                                                                                              <  2:1451084/144‑1 (MQ=255)
                                                             aTAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCg                                                                                                                <  2:876374/119‑1 (MQ=255)
                                                             aTAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCg                                                                                                                >  1:876374/1‑119 (MQ=255)
                                                               aGGCGCCACAAGGGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTaa                                                                              <  1:1333738/150‑1 (MQ=255)
                                                                 gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCgg                                                                                  >  2:595827/1‑144 (MQ=255)
                                                                 gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCgg                                                                                  <  1:595827/144‑1 (MQ=255)
                                                                 gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATg                                                                            <  2:1867900/150‑1 (MQ=255)
                                                                    ccACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                              >  2:317407/1‑82 (MQ=255)
                                                                    ccACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                              <  1:317407/82‑1 (MQ=255)
                                                                    ccACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCtt                                                                         >  1:836322/1‑150 (MQ=255)
                                                                      acaAGTGCCCTGACGAACACCTTGCAGCGCTGACACATTCCCAAGCACCATAAAAAAGCCCCGGCTCATGCGAGACGGGGCTTTTTTATGTCGCTTGCTGCAGTCAGCCGTTGAACACTTCCTCCACGCTGGTCAGCGGGTAATGCTTCg                                                                       >  1:1681488/1‑150 (MQ=255)
                                                                      acaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCg                                                                       >  1:1293625/1‑150 (MQ=255)
                                                                       caAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCgg                                                                      <  2:195308/150‑1 (MQ=255)
                                                                               cTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCAGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCa                                                              >  2:412660/1‑150 (MQ=255)
                                                                               cTGACCAACACCTAGCAACCCTGACACATTCCCCAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGCGCAACTTTATCTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGCTAAGGcc                                                              >  2:914123/1‑149 (MQ=255)
                                                                                  aCCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAggg                                                           <  1:1168168/150‑1 (MQ=255)
                                                                                   ccAACACCTTGCAGCCCTGCCACATGCCCAAGCCCCATAACAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCCGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGt                                                          <  1:1002887/150‑1 (MQ=255)
                                                                                    cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTa                                                         >  1:1027433/1‑150 (MQ=255)
                                                                                    cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTa                                                         <  1:1627431/150‑1 (MQ=255)
                                                                                     aaCACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAg                                                        >  2:1427884/1‑150 (MQ=255)
                                                                                      acacCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGc                                                       <  2:1412639/150‑1 (MQ=255)
                                                                                      acacCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTACATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGc                                                       >  1:681264/1‑150 (MQ=255)
                                                                                         ccTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCGTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCAc                                                    >  1:47421/1‑150 (MQ=255)
                                                                                          cTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACg                                                   >  1:560239/1‑150 (MQ=255)
                                                                                                 cccTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt                                            >  2:1227250/1‑150 (MQ=255)
                                                                                                 cccTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt                                            <  2:1837127/150‑1 (MQ=255)
                                                                                                   cTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCg                                          >  2:1671378/1‑150 (MQ=255)
                                                                                                    tGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGa                                         >  1:1285082/1‑150 (MQ=255)
                                                                                                      acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATg                                       <  2:1293625/150‑1 (MQ=255)
                                                                                                          aTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAg                                   >  1:1331948/1‑150 (MQ=255)
                                                                                                              ccAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTGCGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCtt                               >  1:959280/1‑150 (MQ=255)
                                                                                                               cAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTg                              <  2:1443715/150‑1 (MQ=255)
                                                                                                               cAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTg                              >  1:726673/1‑150 (MQ=255)
                                                                                                                    gacATAAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCcac                         >  2:716677/3‑150 (MQ=38)
                                                                                                                     ccAGAAAAAGGCCCCGGGTCATCCGAGCCGGGGATTTTCGATGCCGGGTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCcaca                        <  2:959280/150‑1 (MQ=255)
                                                                                                                      cATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAg                       >  1:1098365/1‑150 (MQ=255)
                                                                                                                         aaaaaaGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTTGATGGCAGCCTTGGCCACAGCAt                    >  2:1391646/1‑150 (MQ=39)
                                                                                                                          aaaaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATc                   <  2:1331948/150‑1 (MQ=255)
                                                                                                                             aaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGAATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCGCAGCATAgg                 >  2:517076/1‑150 (MQ=255)
                                                                                                                                 ccccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac             >  2:1166430/1‑150 (MQ=255)
                                                                                                                                  cccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTTGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGc                                                       <  1:1467797/107‑1 (MQ=255)
                                                                                                                                  cccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTTGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGc                                                       >  2:1467797/1‑107 (MQ=255)
                                                                                                                                  cccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                        >  1:1631842/1‑90 (MQ=255)
                                                                                                                                  cccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                        <  2:1631842/90‑1 (MQ=255)
                                                                                                                                   ccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGc                                                                                                                  >  1:1342610/1‑48 (MQ=255)
                                                                                                                                   ccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGc                                                                                                                  <  2:1342610/48‑1 (MQ=255)
                                                                                                                                   ccGACTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGATGCAATCGCTGCTTGATCACTTCATCCACGCTTTGCAGCGGGTAATGCTTACGGTAAGGCGGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacga           >  1:976236/1‑150 (MQ=255)
                                                                                                                                     ggCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacgacg         <  2:560239/150‑1 (MQ=255)
                                                                                                                                       cTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGGAACGACGGt       <  2:270807/150‑1 (MQ=255)
                                                                                                                                       cTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGt       <  2:836322/150‑1 (MQ=255)
                                                                                                                                          aTGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAgg                                                            >  1:1335623/1‑94 (MQ=255)
                                                                                                                                          aTGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAgg                                                            <  2:1335623/94‑1 (MQ=255)
                                                                                                                                          aTGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGAt    <  1:1166430/150‑1 (MQ=255)
                                                                                                                                            gCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAg                                                        <  1:847585/96‑1 (MQ=255)
                                                                                                                                            gCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAg                                                        >  2:847585/1‑96 (MQ=255)
                                                                                                                                            gCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGATCa  >  2:1831075/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
GCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGATCA  >  NC_002947/5807903‑5808191

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