Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 436,707:1 +C 5.9% coding (119/519 nt) pgpA → phosphatidylglycerophosphatase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,7071.C 5.9% 276.1 / 12.6 68coding (119/519 nt)pgpAphosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base . (35/29);  new base C (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

AACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGGCAGC  >  NC_000913/436562‑436855
                                                                                                                                                  |                                                                                                                                                    
aaCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGAT‑cgggc                                                                                                                                                 <  1:953/150‑3 (MQ=255)
 aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCt                                                                                                                                                 <  2:1365614/150‑1 (MQ=255)
 aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGATCGGGCt                                                                                                                                                 <  2:858171/150‑1 (MQ=255)
  ccTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGATCGGGCtt                                                                                                                                                <  1:650620/150‑2 (MQ=255)
     gTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCt                                                                                                                                             <  2:670942/150‑1 (MQ=255)
       tACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTgg                                                                                                                                           >  2:723645/1‑150 (MQ=255)
        aCATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGc                                                                                                                                          >  2:567089/1‑150 (MQ=255)
          aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAg                                                                                                                                        <  2:14317/150‑1 (MQ=255)
          aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAg                                                                                                                                        <  2:1426574/150‑1 (MQ=255)
             aGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGa                                                                                                                                     >  1:205639/1‑150 (MQ=255)
                cTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTc                                                                                                                                  >  1:1325229/1‑150 (MQ=255)
                        gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑g                                                                                                                                                     >  1:1012577/1‑123 (MQ=255)
                        gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑g                                                                                                                                                     <  2:1012577/123‑1 (MQ=255)
                        gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTgg                                                                                                                          <  2:223072/150‑1 (MQ=255)
                            tGACCTTTTTGCGACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGc                                                                                                                                       <  1:955562/133‑1 (MQ=255)
                            tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGc                                                                                                                                       >  2:955562/1‑133 (MQ=255)
                            tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGt                                                                                                                               >  2:659949/1‑141 (MQ=255)
                            tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGt                                                                                                                               <  1:659949/141‑1 (MQ=255)
                                  ttttGCCACGCCATAATGATGTCGCGAAATGTCGCCTGAAGATGAGTAATCGGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCGCGTTCGTTCCTGGGTCGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgatga                                                                                                                <  2:991045/150‑1 (MQ=255)
                                    ttGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGAcc                                                                                                              >  2:958991/1‑150 (MQ=255)
                                       ccACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttt                                                                                                           >  1:510815/1‑150 (MQ=255)
                                        cACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCtttt                                                                                                          >  2:1396252/1‑150 (MQ=255)
                                           gCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGc                                                                                                       <  1:723645/150‑1 (MQ=255)
                                              aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                                    >  2:461715/1‑150 (MQ=255)
                                              aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                                    >  1:259525/1‑150 (MQ=255)
                                                 aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                                 >  1:650273/1‑150 (MQ=255)
                                                 aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGGCGGATTCGGAAGTGGATTAAGCCCGAGCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                                 <  1:1252892/150‑1 (MQ=255)
                                                  aGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGGTCTGGTATCTGATGACCTTTTTGCCCTGGCa                                                                                                <  2:1289428/150‑1 (MQ=255)
                                                      gTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctc                                                                                            >  1:1087486/1‑150 (MQ=255)
                                                           gAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTc                                                                                       <  1:93280/150‑1 (MQ=255)
                                                                gTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                                  >  1:953538/1‑150 (MQ=255)
                                                                   gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                                               >  2:114157/1‑150 (MQ=255)
                                                                   gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                                               >  1:669892/1‑150 (MQ=255)
                                                                    ccTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTATTTGCCCTGGCAGCTCTACTCGCTggtggt                                                                              >  1:1356610/1‑150 (MQ=255)
                                                                      tGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGa                                                                            >  1:776567/1‑150 (MQ=255)
                                                                          gATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                        <  2:510815/150‑1 (MQ=255)
                                                                          gATGAGTAATCCGTGGCATCTACTTGCTGTCAGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                        <  1:1385744/150‑1 (MQ=255)
                                                                            tGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                                           >  1:206344/1‑145 (MQ=255)
                                                                            tGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                                           <  2:206344/145‑1 (MQ=255)
                                                                               gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                                   >  1:787611/1‑150 (MQ=255)
                                                                               gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                                   >  1:115409/1‑150 (MQ=255)
                                                                                 aaTCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATc                                                                 <  1:114157/150‑1 (MQ=255)
                                                                                   tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTg                                                               <  1:567089/150‑1 (MQ=255)
                                                                                   tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTg                                                               <  2:865042/150‑1 (MQ=255)
                                                                                           aTCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCg                                                       <  1:1382370/150‑1 (MQ=255)
                                                                                              tACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCt                                                    >  1:1144659/1‑150 (MQ=255)
                                                                                                  tGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCt                                                >  1:274089/1‑150 (MQ=255)
                                                                                                    cTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCttt                                              >  1:162978/1‑150 (MQ=255)
                                                                                                      gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                            <  1:958991/150‑1 (MQ=255)
                                                                                                      gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                            >  1:350698/1‑150 (MQ=255)
                                                                                                                gAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTg                                                                                                        >  2:794575/1‑80 (MQ=255)
                                                                                                                gAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTg                                                                                                        <  1:794575/80‑1 (MQ=255)
                                                                                                                  aGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCg                                >  2:156866/1‑150 (MQ=255)
                                                                                                                     ggATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                             <  2:274089/150‑1 (MQ=255)
                                                                                                                      gATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAg                            <  2:259525/150‑1 (MQ=255)
                                                                                                                            gCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                      >  2:433100/1‑150 (MQ=255)
                                                                                                                             gcccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                      >  1:555406/5‑150 (MQ=255)
                                                                                                                             cccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg                     <  2:1144659/150‑1 (MQ=255)
                                                                                                                             cccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTAATCAAACGGCGAAAGACATggg                     >  1:454955/1‑150 (MQ=255)
                                                                                                                              cccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg                     >  2:363304/4‑150 (MQ=255)
                                                                                                                                gATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTTTCATCAAACGGCGAAAGACATGGgtgt                  >  1:413306/1‑150 (MQ=255)
                                                                                                                                gATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGgtgt                  >  1:188198/1‑150 (MQ=255)
                                                                                                                                          gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                               >  2:16084/1‑127 (MQ=255)
                                                                                                                                          gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                               <  1:16084/127‑1 (MQ=255)
                                                                                                                                            gACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATgg      <  2:162978/150‑1 (MQ=255)
                                                                                                                                            gACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATgg      <  1:1413827/150‑1 (MQ=255)
                                                                                                                                               gAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                               <  2:837860/122‑1 (MQ=255)
                                                                                                                                               gAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                               >  1:837860/1‑122 (MQ=255)
                                                                                                                                               gAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGa                           <  1:1396973/126‑1 (MQ=255)
                                                                                                                                               gAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGa                           >  2:1396973/1‑126 (MQ=255)
                                                                                                                                                aT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGgcagc  >  2:842185/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGGCAGC  >  NC_000913/436562‑436855

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