Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,446,647 T→C T54A (ACC→GCC)  rplP ← 50S ribosomal subunit protein L16

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,446,6470TC99.6% 838.0 / NA 240T54A (ACC→GCC) rplP50S ribosomal subunit protein L16
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (132/107);  total (133/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

TACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCG  >  NC_000913/3446557‑3446736
                                                                                          |                                                                                         
taccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCa                                                                                         >  2:1835120/1‑93 (MQ=255)
  cccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATa                                                                                       <  2:5258373/93‑1 (MQ=255)
   ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAg                                                                                      >  1:274512/1‑93 (MQ=255)
   ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAg                                                                                      <  2:4230198/93‑1 (MQ=255)
   ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAg                                                                                      <  1:5884644/93‑1 (MQ=255)
   ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAg                                                                                      <  1:2304501/93‑1 (MQ=255)
    cATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGc                                                                                     <  1:4622702/93‑1 (MQ=255)
    cATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGc                                                                                     >  1:5582202/1‑93 (MQ=255)
    cATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGc                                                                                     >  1:182227/1‑93 (MQ=255)
      tACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCac                                                                                   >  1:1465806/1‑93 (MQ=255)
      tACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCac                                                                                   <  2:1573102/93‑1 (MQ=255)
        cGCCCTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacga                                                                                 <  1:5268996/93‑1 (MQ=255)
        cGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacga                                                                                 >  2:272326/1‑93 (MQ=255)
        cGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacga                                                                                 <  2:4829100/93‑1 (MQ=255)
         gCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacgac                                                                                <  1:1069659/93‑1 (MQ=255)
         gCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacgac                                                                                <  1:5053705/93‑1 (MQ=255)
         gCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacgac                                                                                <  2:4318935/93‑1 (MQ=255)
         gCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacgac                                                                                <  1:1612716/93‑1 (MQ=255)
          cACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCacgacg                                                                               <  2:2345650/93‑1 (MQ=255)
           aCTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGt                                                                              >  1:923592/1‑93 (MQ=255)
             tGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGtgc                                                                            >  1:3072598/1‑93 (MQ=255)
             tGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGtgc                                                                            <  2:3358066/93‑1 (MQ=255)
              gCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGtgct                                                                           <  2:2651326/93‑1 (MQ=255)
              gCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGtgct                                                                           <  2:4046349/93‑1 (MQ=255)
                cAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGtgctgc                                                                         >  1:587533/1‑93 (MQ=255)
                 aGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTgctgct                                                                        >  1:326771/1‑93 (MQ=255)
                 aGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTgctgct                                                                        <  1:3364193/93‑1 (MQ=255)
                  gCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCtt                                                                       >  1:6414956/1‑93 (MQ=255)
                   cGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTc                                                                      <  1:6002374/93‑1 (MQ=255)
                     gCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGa                                                                    >  2:5079272/1‑93 (MQ=255)
                       ttttCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATc                                                                  >  2:402282/1‑93 (MQ=255)
                        tttCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCt                                                                 >  1:4993275/1‑93 (MQ=255)
                        tttCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCt                                                                 >  2:1315295/1‑93 (MQ=255)
                         ttCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTg                                                                <  2:628161/93‑1 (MQ=255)
                         ttCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTg                                                                <  1:5792728/93‑1 (MQ=255)
                          tCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGa                                                               <  2:1161695/93‑1 (MQ=255)
                           cAGTGATCGTTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGAc                                                              <  1:3607951/93‑1 (MQ=255)
                           cAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGAc                                                              >  1:5220591/1‑93 (MQ=255)
                           cAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGAc                                                              >  1:6079422/1‑93 (MQ=255)
                            aGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACg                                                             >  1:1116691/1‑93 (MQ=255)
                             gTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACgg                                                            <  2:5340835/93‑1 (MQ=255)
                               gATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGc                                                          <  2:547280/93‑1 (MQ=255)
                               gATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGc                                                          <  1:2211446/93‑1 (MQ=255)
                                 tCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAg                                                        >  1:874427/1‑93 (MQ=255)
                                      ttGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAga                                                   <  2:295474/93‑1 (MQ=255)
                                         tCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAgacga                                                >  2:571723/1‑93 (MQ=255)
                                         tCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAgacga                                                >  1:4408443/1‑93 (MQ=255)
                                             ggAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCAc                                            >  2:4490443/1‑93 (MQ=255)
                                              gAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACg                                           <  1:679985/93‑1 (MQ=255)
                                              gAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACg                                           >  2:4802264/1‑93 (MQ=255)
                                               aaCACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACgg                                          >  2:1465086/1‑93 (MQ=255)
                                               aaCACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACgg                                          <  2:4378457/93‑1 (MQ=255)
                                               aaCACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACgc                                          >  2:4209230/1‑92 (MQ=255)
                                                acacacGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGc                                         >  2:1923008/1‑93 (MQ=255)
                                                acacacGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGc                                         >  1:1503580/1‑93 (MQ=255)
                                                  acacGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCa                                       <  2:3733216/93‑1 (MQ=255)
                                                  acacGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCa                                       >  2:5243336/1‑93 (MQ=255)
                                                   cacGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGTCaa                                      >  1:4769111/1‑93 (MQ=255)
                                                    acGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAAc                                     >  2:4643274/1‑93 (MQ=255)
                                                       gATCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATGTGACGGGCAGTCAGACGACCACGGCCAACAGc                                  <  1:5717445/93‑1 (MQ=255)
                                                        aTCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCt                                 >  2:65639/1‑93 (MQ=255)
                                                        aTCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCt                                 >  2:6277770/1‑93 (MQ=255)
                                                         tCCAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCtt                                >  2:5863210/1‑93 (MQ=255)
                                                           cAGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTc                              >  2:2179484/1‑93 (MQ=255)
                                                            aGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCa                             <  2:3004696/93‑1 (MQ=255)
                                                            aGATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCa                             <  2:3806396/93‑1 (MQ=255)
                                                             gATCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAg                            >  2:3635624/1‑93 (MQ=255)
                                                              aTCTTACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGa                           <  2:2076978/93‑1 (MQ=255)
                                                                  tACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGa                       <  2:5029532/93‑1 (MQ=255)
                                                                  tACCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGa                       <  1:1192985/93‑1 (MQ=255)
                                                                   aCCTTGACGCTTAACTGCACGGGCCATAGCACGTCGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGaa                      >  1:1962517/1‑93 (MQ=255)
                                                                   aCCTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGaa                      >  1:5705854/1‑93 (MQ=255)
                                                                    ccTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAg                     >  2:1084943/1‑93 (MQ=255)
                                                                     cTTGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGc                    >  1:4209255/1‑93 (MQ=255)
                                                                      ttGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCt                   >  2:732711/1‑93 (MQ=255)
                                                                      ttGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCt                   >  2:2313320/1‑93 (MQ=255)
                                                                       tGACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTg                  >  1:1747433/1‑93 (MQ=255)
                                                                        gACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGc                 <  1:3781925/93‑1 (MQ=255)
                                                                        gACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGc                 >  1:4152719/1‑93 (MQ=255)
                                                                        gACGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGc                 >  1:4462305/1‑93 (MQ=255)
                                                                         aCGCTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGcc                <  2:3063580/93‑1 (MQ=255)
                                                                            cTTAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGaa             <  1:2905347/93‑1 (MQ=255)
                                                                             ttAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAg            >  1:5106463/1‑93 (MQ=255)
                                                                              tAACTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGc           <  2:936568/93‑1 (MQ=255)
                                                                               aaCTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCt          >  1:2518623/1‑93 (MQ=255)
                                                                               aaCTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCt          <  2:802745/93‑1 (MQ=255)
                                                                                aCTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTa         >  2:2222000/1‑93 (MQ=255)
                                                                                aCTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTa         >  1:3090360/1‑93 (MQ=255)
                                                                                aCTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTa         >  2:3150913/1‑93 (MQ=255)
                                                                                 cTGCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTaa        <  2:3988412/93‑1 (MQ=255)
                                                                                   gCACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACa      <  1:2403196/93‑1 (MQ=255)
                                                                                    cACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACAt     <  1:4558163/93‑1 (MQ=255)
                                                                                    cACGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACAt     >  2:5163148/1‑93 (MQ=255)
                                                                                     aCGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATc    <  2:3834102/93‑1 (MQ=255)
                                                                                     aCGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATc    <  2:3930402/93‑1 (MQ=255)
                                                                                     aCGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGCCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATc    >  1:15746/1‑93 (MQ=255)
                                                                                      cGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATcc   >  1:5352580/1‑93 (MQ=255)
                                                                                      cGGGCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATcc   >  1:3525111/1‑93 (MQ=255)
                                                                                       gggCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCg  >  1:3758859/1‑93 (MQ=255)
                                                                                       gggCCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCg  <  1:3384776/93‑1 (MQ=255)
                                                                                          |                                                                                         
TACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCG  >  NC_000913/3446557‑3446736

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