Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,348,653 (T)6→7 coding (379/1077 nt) glpQ ← glycerophosphodiester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,348,6471.T86.9% 847.4 / 114.2 251coding (385/1077 nt)glpQglycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base . (16/17);  new base T (120/98);  total (136/115)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

Only 100 of 285 total aligned reads displayed.

ATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGACCGTCTT  >  NC_000913/2348555‑2348739
                                                                                             |                                                                                            
 tttAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACt                                                                                             >  1:5439966/1‑93 (MQ=255)
 tttAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACt                                                                                             >  2:3136738/1‑93 (MQ=255)
    aaCCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACtttt                                                                                          <  1:6438011/93‑1 (MQ=255)
      ccccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACtttttt                                                                                        <  2:7451146/93‑1 (MQ=255)
      ccccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACtttttt                                                                                        >  1:931724/1‑93 (MQ=255)
        ccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACttttttta                                                                                      >  2:2896403/1‑91 (MQ=255)
        ccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTAc                                                                                      <  1:6588766/93‑1 (MQ=255)
          tGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACtttttttacc                                                                                    >  1:5277394/1‑89 (MQ=255)
          tGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACtttttttacc                                                                                    >  1:2762309/1‑89 (MQ=255)
          tGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCg                                                                                    >  1:3326436/1‑93 (MQ=255)
            aaCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGtt                                                                                  <  2:3039995/93‑1 (MQ=255)
            aaCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGtt                                                                                  <  2:4493088/93‑1 (MQ=255)
             aCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGtt                                                                                  <  2:783109/93‑1 (MQ=255)
             aCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGttt                                                                                 >  1:5712091/1‑93 (MQ=255)
                  ttCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCaa                                                                             <  2:671480/93‑1 (MQ=255)
                  ttCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCaa                                                                             >  1:255599/1‑93 (MQ=255)
                   tCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAat                                                                            <  2:3288415/93‑1 (MQ=255)
                    cAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAata                                                                           >  2:2663646/1‑93 (MQ=255)
                     aaTCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAatat                                                                          <  2:4799978/93‑1 (MQ=255)
                      aTCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATc                                                                         <  2:1004577/93‑1 (MQ=255)
                       tctcTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCg                                                                        >  2:3836211/1‑93 (MQ=255)
                       tctcTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGa                                                                       >  2:559573/1‑93 (MQ=255)
                         tcttctTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGaa                                                                      >  1:7187165/1‑93 (MQ=255)
                         tcttctTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGaa                                                                      >  1:4927156/1‑93 (MQ=255)
                          cttcttCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGaaa                                                                     >  2:5816997/1‑93 (MQ=255)
                           ttcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAAc                                                                    <  1:6297602/93‑1 (MQ=255)
                           ttcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAAc                                                                    >  1:5888583/1‑93 (MQ=255)
                            tcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAAcc                                                                   >  2:3750618/1‑93 (MQ=255)
                            tcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAAcc                                                                   >  1:2717093/1‑93 (MQ=255)
                            tcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCt                                                                  <  2:2530975/93‑1 (MQ=255)
                             cttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCtt                                                                 <  1:4757146/93‑1 (MQ=255)
                              ttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCtt                                                                 >  1:3779715/1‑93 (MQ=255)
                              ttcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCtt                                                                 >  2:1650764/1‑93 (MQ=255)
                               tcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTc                                                                <  2:5114908/93‑1 (MQ=255)
                               tcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTc                                                                <  1:6364774/93‑1 (MQ=255)
                               tcAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCg                                                               <  1:1000637/93‑1 (MQ=255)
                                cAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCg                                                               >  1:2243077/1‑93 (MQ=255)
                                cAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCg                                                               >  2:3677917/1‑93 (MQ=255)
                                 aaaGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCgg                                                              >  2:236276/1‑93 (MQ=255)
                                 aaaGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCgg                                                              >  2:1488095/1‑93 (MQ=255)
                                 aaaGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCgg                                                              >  1:5514024/1‑93 (MQ=255)
                                  aaGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGt                                                             <  2:5355996/93‑1 (MQ=255)
                                   aGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTa                                                            >  1:1703742/1‑93 (MQ=255)
                                     gtgtGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTaaa                                                          <  1:1197946/93‑1 (MQ=255)
                                     gtgtGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTaaa                                                          <  1:4171682/93‑1 (MQ=255)
                                      tgtgCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAt                                                         <  2:6784653/93‑1 (MQ=255)
                                       gtgCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAtt                                                        <  1:7375939/93‑1 (MQ=255)
                                         gCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTc                                                      >  2:2111261/1‑93 (MQ=255)
                                           aCCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCaa                                                    >  2:5849764/1‑93 (MQ=255)
                                               ggAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGAc                                                >  1:7020382/1‑93 (MQ=255)
                                               ggAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGAc                                                <  1:1529798/93‑1 (MQ=255)
                                               ggAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGAc                                                >  2:6326821/1‑93 (MQ=255)
                                               ggAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGAc                                                <  2:1814174/93‑1 (MQ=255)
                                               ggAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACt                                               >  1:3656462/1‑93 (MQ=255)
                                                gAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACt                                               >  2:5937713/1‑93 (MQ=255)
                                                   gTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTaa                                            >  2:1583889/1‑93 (MQ=255)
                                                   gTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTaa                                            >  1:1397701/1‑93 (MQ=255)
                                                      ggACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAttt                                         <  1:4484489/93‑1 (MQ=255)
                                                      ggACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAttt                                         <  2:6904207/93‑1 (MQ=255)
                                                        aCTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCa                                       >  1:2937706/1‑93 (MQ=255)
                                                         cTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCAt                                      >  1:3271384/1‑93 (MQ=255)
                                                         cTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCAt                                      >  2:3005925/1‑93 (MQ=255)
                                                          ttACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATc                                     >  2:4700290/1‑93 (MQ=255)
                                                          ttACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATcc                                    <  2:6517575/93‑1 (MQ=255)
                                                            gcccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAGCCTTCGGTAAATTTCAACGACTTAATTTCATCCa                                   <  1:1119211/92‑1 (MQ=255)
                                                            gcccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCa                                   <  1:5122370/92‑1 (MQ=255)
                                                             cccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAg                                  >  2:4225444/1‑93 (MQ=255)
                                                             cccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTCCATCCAg                                  >  1:4328070/1‑93 (MQ=255)
                                                              ccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGc                                 <  1:1071952/93‑1 (MQ=255)
                                                              ccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGc                                 >  1:3950570/1‑93 (MQ=255)
                                                              ccATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGc                                 >  1:4539567/1‑93 (MQ=255)
                                                               gaTTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCg                                >  2:7321943/2‑93 (MQ=255)
                                                               cATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGt                               >  2:1707245/1‑93 (MQ=255)
                                                                aTTGGGAAACGCCCAGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTg                              >  2:6446090/1‑93 (MQ=255)
                                                                 ttGGGAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTg                              >  1:201012/1‑93 (MQ=255)
                                                                    ggAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGaaa                           <  2:1339627/93‑1 (MQ=255)
                                                                    ggAAACGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGaaa                           <  2:5476985/93‑1 (MQ=255)
                                                                      aaaCGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATc                         >  1:1825024/1‑93 (MQ=255)
                                                                       aaCGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCt                        >  1:6683255/1‑93 (MQ=255)
                                                                        aCGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTa                       >  1:594008/1‑93 (MQ=255)
                                                                         cGCCCCGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTAt                      <  2:6217965/93‑1 (MQ=255)
                                                                           ccccGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATcg                    <  1:215009/93‑1 (MQ=255)
                                                                              cGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgt                 >  1:6842051/1‑93 (MQ=255)
                                                                              cGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgt                 >  1:7076284/1‑93 (MQ=255)
                                                                              cGGATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgt                 >  2:7062729/1‑93 (MQ=255)
                                                                               ggATAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgta                <  1:6041222/93‑1 (MQ=255)
                                                                                  tAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgtagta             >  1:461361/1‑93 (MQ=255)
                                                                                  tAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgtagta             >  2:7252123/1‑93 (MQ=255)
                                                                                  tAAGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCgtagta             >  1:5360459/1‑93 (MQ=255)
                                                                                  tAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTaa            <  1:7501997/93‑1 (MQ=255)
                                                                                   aaGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTaa            <  1:7319522/93‑1 (MQ=255)
                                                                                    aGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAAc           >  2:5997871/1‑93 (MQ=255)
                                                                                    aGTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAAc           <  2:3099495/93‑1 (MQ=255)
                                                                                     gTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACg          >  1:2469901/1‑93 (MQ=255)
                                                                                     gTCTGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACg          >  1:1570000/1‑93 (MQ=255)
                                                                                      nntGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGa         <  2:3495013/91‑1 (MQ=255)
                                                                                        tGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGAcc       >  2:5028128/1‑93 (MQ=255)
                                                                                        tGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGAcc       <  1:6170114/93‑1 (MQ=255)
                                                                                        tGCACTTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGAcc       >  1:131304/1‑93 (MQ=255)
                                                                                            cTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGACCGTCtt  <  1:3163126/93‑1 (MQ=255)
                                                                                             |                                                                                            
ATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAATCTATCGCGTAGTAACGACCGTCTT  >  NC_000913/2348555‑2348739

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