Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 436,699 T→A 7.7% P37P (CCT→CCA pgpA → phosphatidylglycerophosphatase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,6990TA 7.7% 198.4 / 5.9 65P37P (CCT→CCApgpAphosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base T (34/26);  new base A (2/3);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

TTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGT  >  NC_000913/436550‑436839
                                                                                                                                                     |                                                                                                                                             
ttCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCt                                                                                                                                               >  1:623045/1‑150 (MQ=255)
 tCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTg                                                                                                                                              <  2:949433/150‑1 (MQ=255)
   gTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTggg                                                                                                                                            >  1:641665/1‑150 (MQ=255)
    tGACGGCGAACCTGTTACATTAGACTGGAAAGGAAATGACCATTTTGCCACGCCATAAAGATTTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGAGGCATCTACTTGCTGTCGGATTCGGAACTGGATTAAGGCCGATCGTTCCTCGGa                                                                                                                                           >  2:789979/1‑150 (MQ=255)
            aaCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATcgggc                                                                                                                                 <  1:953/150‑3 (MQ=255)
             aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCt                                                                                                                                 <  2:1365614/150‑1 (MQ=255)
             aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATCGGGCt                                                                                                                                 <  2:858171/150‑1 (MQ=255)
              ccTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATCGGGCtt                                                                                                                                <  1:650620/150‑2 (MQ=255)
                 gTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCt                                                                                                                             <  2:670942/150‑1 (MQ=255)
                   tACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTgg                                                                                                                           >  2:723645/1‑150 (MQ=255)
                    aCATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGc                                                                                                                          >  2:567089/1‑150 (MQ=255)
                      aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAg                                                                                                                        <  2:14317/150‑1 (MQ=255)
                      aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAg                                                                                                                        <  2:1426574/150‑1 (MQ=255)
                         aGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGa                                                                                                                     >  1:205639/1‑150 (MQ=255)
                            cTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTc                                                                                                                  >  1:1325229/1‑150 (MQ=255)
                                    gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATg                                                                                                                                     >  1:1012577/1‑123 (MQ=255)
                                    gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATg                                                                                                                                     <  2:1012577/123‑1 (MQ=255)
                                    gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTgg                                                                                                          <  2:223072/150‑1 (MQ=255)
                                        tGACCTTTTTGCGACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGc                                                                                                                       <  1:955562/133‑1 (MQ=255)
                                        tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGc                                                                                                                       >  2:955562/1‑133 (MQ=255)
                                        tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGt                                                                                                               <  1:659949/141‑1 (MQ=255)
                                        tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGt                                                                                                               >  2:659949/1‑141 (MQ=255)
                                              ttttGCCACGCCATAATGATGTCGCGAAATGTCGCCTGAAGATGAGTAATCGGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCGCGTTCGTTCCTGGGTCGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgatga                                                                                                <  2:991045/150‑1 (MQ=255)
                                                ttGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGAcc                                                                                              >  2:958991/1‑150 (MQ=255)
                                                   ccACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttt                                                                                           >  1:510815/1‑150 (MQ=255)
                                                    cACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCtttt                                                                                          >  2:1396252/1‑150 (MQ=255)
                                                       gCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGc                                                                                       <  1:723645/150‑1 (MQ=255)
                                                          aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                    >  2:461715/1‑150 (MQ=255)
                                                          aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                    >  1:259525/1‑150 (MQ=255)
                                                             aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                 >  1:650273/1‑150 (MQ=255)
                                                             aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGGCGGATTCGGAAGTGGATTAAGCCCGAGCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                 <  1:1252892/150‑1 (MQ=255)
                                                              aGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGGTCTGGTATCTGATGACCTTTTTGCCCTGGCa                                                                                <  2:1289428/150‑1 (MQ=255)
                                                                  gTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctc                                                                            >  1:1087486/1‑150 (MQ=255)
                                                                       gAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTc                                                                       <  1:93280/150‑1 (MQ=255)
                                                                            gTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                  >  1:953538/1‑150 (MQ=255)
                                                                               gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                               >  2:114157/1‑150 (MQ=255)
                                                                               gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                               >  1:669892/1‑150 (MQ=255)
                                                                                ccTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTATTTGCCCTGGCAGCTCTACTCGCTggtggt                                                              >  1:1356610/1‑150 (MQ=255)
                                                                                  tGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGa                                                            >  1:776567/1‑150 (MQ=255)
                                                                                      gATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                        <  2:510815/150‑1 (MQ=255)
                                                                                      gATGAGTAATCCGTGGCATCTACTTGCTGTCAGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                        <  1:1385744/150‑1 (MQ=255)
                                                                                        tGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                           >  1:206344/1‑145 (MQ=255)
                                                                                        tGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                           <  2:206344/145‑1 (MQ=255)
                                                                                           gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                   >  1:115409/1‑150 (MQ=255)
                                                                                           gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                   >  1:787611/1‑150 (MQ=255)
                                                                                             aaTCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATc                                                 <  1:114157/150‑1 (MQ=255)
                                                                                               tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTg                                               <  1:567089/150‑1 (MQ=255)
                                                                                               tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTg                                               <  2:865042/150‑1 (MQ=255)
                                                                                                       aTCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCg                                       <  1:1382370/150‑1 (MQ=255)
                                                                                                          tACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCt                                    >  1:1144659/1‑150 (MQ=255)
                                                                                                              tGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCt                                >  1:274089/1‑150 (MQ=255)
                                                                                                                cTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCttt                              >  1:162978/1‑150 (MQ=255)
                                                                                                                  gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                            <  1:958991/150‑1 (MQ=255)
                                                                                                                  gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                            >  1:350698/1‑150 (MQ=255)
                                                                                                                            gAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTg                                                                                        >  2:794575/1‑80 (MQ=255)
                                                                                                                            gAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTg                                                                                        <  1:794575/80‑1 (MQ=255)
                                                                                                                              aGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCg                >  2:156866/1‑150 (MQ=255)
                                                                                                                                 ggATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa             <  2:274089/150‑1 (MQ=255)
                                                                                                                                  gATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAg            <  2:259525/150‑1 (MQ=255)
                                                                                                                                        gCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg      >  2:433100/1‑150 (MQ=255)
                                                                                                                                         gcccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg      >  1:555406/5‑150 (MQ=255)
                                                                                                                                         cccGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg     <  2:1144659/150‑1 (MQ=255)
                                                                                                                                         cccGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTAATCAAACGGCGAAAGACATggg     >  1:454955/1‑150 (MQ=255)
                                                                                                                                          cccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg     >  2:363304/4‑150 (MQ=255)
                                                                                                                                            gATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTTTCATCAAACGGCGAAAGACATGGgtgt  >  1:413306/1‑150 (MQ=255)
                                                                                                                                            gATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGgtgt  >  1:188198/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGT  >  NC_000913/436550‑436839

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: