Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,785,382 G→A T306T (ACG→ACA b3613 → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,785,3820GA97.3% 631.2 / 0.6 184T306T (ACG→ACAb3613putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (80/99);  total (82/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 189 total aligned reads displayed.

AACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAT  >  NC_000913/3785291‑3785471
                                                                                           |                                                                                         
 accgaccgAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACt                                                                                         <  2:683139/93‑1 (MQ=255)
 accgaccgAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACt                                                                                         >  2:3018208/1‑93 (MQ=255)
  ccgaccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTg                                                                                        >  2:723930/1‑93 (MQ=255)
    gaccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCa                                                                                      >  1:4058958/1‑93 (MQ=255)
    gaccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCa                                                                                      >  1:951680/1‑93 (MQ=255)
     accgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCAt                                                                                     <  1:3994252/93‑1 (MQ=255)
      ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCCTc                                                                                    >  2:2688812/1‑93 (MQ=255)
      ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATc                                                                                    >  1:3310356/1‑93 (MQ=255)
      ccgaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATc                                                                                    <  2:3175279/93‑1 (MQ=255)
        gaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGc                                                                                  <  1:225802/93‑1 (MQ=255)
        gaAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGc                                                                                  <  1:711249/93‑1 (MQ=255)
         aaaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCt                                                                                 >  1:3248422/1‑93 (MQ=255)
          aaGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTa                                                                                >  2:1612634/1‑93 (MQ=255)
           aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTTGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               >  2:1069595/1‑93 (MQ=255)
           aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               <  2:2057395/93‑1 (MQ=255)
           aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTAt                                                                               <  2:2689541/93‑1 (MQ=255)
           aGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTAGCGGGCCGACACTGCATCGCTAt                                                                               <  1:1824155/93‑1 (MQ=255)
            gCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATg                                                                              <  2:1799750/93‑1 (MQ=255)
             ggAAAAATCGCTGATGTCCCGTACGGGTGGACTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATgg                                                                             <  1:3893612/92‑1 (MQ=255)
              gAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGc                                                                            >  1:2827717/1‑93 (MQ=255)
              gAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGc                                                                            >  1:1459352/1‑93 (MQ=255)
              gAAAAATCGCTGATGTCCCGTACCGGCGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGc                                                                            >  2:3377901/1‑93 (MQ=255)
                 aaaTCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGaa                                                                         <  1:1773442/93‑1 (MQ=255)
                   aTCGCTGATGTCCCGTACCGGTGCTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACa                                                                       <  2:2609121/93‑1 (MQ=255)
                     cGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGc                                                                     <  1:1030212/93‑1 (MQ=255)
                      gCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCGGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCt                                                                    <  1:1027444/93‑1 (MQ=255)
                       cTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTa                                                                   <  1:3765309/93‑1 (MQ=255)
                        cgATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTAc                                                                  <  2:4078557/92‑1 (MQ=255)
                         gATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACa                                                                 >  1:2341490/1‑93 (MQ=255)
                           tGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAgg                                                               <  1:1438198/93‑1 (MQ=255)
                            gTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAggg                                                              <  1:3013058/93‑1 (MQ=255)
                            gTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAggg                                                              <  1:453690/93‑1 (MQ=255)
                             tCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGt                                                             <  2:4432082/93‑1 (MQ=255)
                              cccGTCCCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTg                                                            <  2:3778649/93‑1 (MQ=255)
                               ccGTACCGGTGGTCTGGGCGCTCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGGGa                                                           >  1:836115/1‑93 (MQ=255)
                                 gTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAt                                                         <  2:4617559/93‑1 (MQ=255)
                                 gTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAt                                                         <  1:67435/93‑1 (MQ=255)
                                  tACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAtt                                                        >  2:965985/1‑93 (MQ=255)
                                  tACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAAtt                                                        <  2:3357205/93‑1 (MQ=255)
                                   aCCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCGGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTa                                                       <  2:3591838/93‑1 (MQ=255)
                                    ccGGTGGTCTGGGCGCGCCGCGTGGTAAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTAc                                                      <  1:3047025/93‑1 (MQ=255)
                                      ggtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGc                                                    <  1:1175314/93‑1 (MQ=255)
                                       gtggtCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCt                                                   <  1:3991184/93‑1 (MQ=255)
                                            cTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGaaa                                              >  1:4574443/1‑93 (MQ=255)
                                             tGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAg                                             >  1:473034/1‑93 (MQ=255)
                                             tGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAg                                             <  2:254635/93‑1 (MQ=255)
                                             tGGGCAGGCCGTGTGGTCAGGCATTCTGGCCGGTTCGGGGGCCGACACTGCATCTCTATGGCGAAAAGCTACAGGGTGAATTACGCTGGAAAg                                             <  1:4248451/93‑1 (MQ=255)
                                              gggCGCGCCGTGTGTACAGGGATTCTGGCAGGTTGTATGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTCAATTACGCTGGAAAgg                                            <  2:3132386/93‑1 (MQ=255)
                                               ggggcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGt                                           <  1:3115342/90‑1 (MQ=255)
                                               ggCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGt                                           <  1:4315016/93‑1 (MQ=255)
                                                gcgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTa                                          <  1:182628/93‑1 (MQ=255)
                                                gcgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTa                                          <  1:2811709/93‑1 (MQ=255)
                                                gcgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTa                                          <  2:4508303/93‑1 (MQ=255)
                                                 cgcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTAt                                         <  1:2281019/93‑1 (MQ=255)
                                                 cacgccgcGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTAt                                         <  2:2078997/91‑1 (MQ=255)
                                                  gcgcCGCGTGGTCAGGCATTCTGGCCGGTTTGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGTATTACGCTGGAAAGGTATg                                        <  1:1386273/93‑1 (MQ=255)
                                                  gcgcCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATg                                        <  2:293716/93‑1 (MQ=255)
                                                    gccgcGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGt                                      <  1:41476/93‑1 (MQ=255)
                                                      cgcgTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTa                                    <  1:2730732/93‑1 (MQ=255)
                                                        ggTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATc                                  <  2:1051224/92‑1 (MQ=255)
                                                          tGGTGAGGCATTCTGGCCGGTTCGCGGGCCGGCACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                <  2:1452134/93‑1 (MQ=255)
                                                          tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                >  1:2520769/1‑93 (MQ=255)
                                                          tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                >  2:1635985/1‑93 (MQ=255)
                                                          tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                >  2:497266/1‑93 (MQ=255)
                                                          tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCgg                                >  1:790491/1‑93 (MQ=255)
                                                          tGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGAGAACAGCTACAGGGCGAATTACGCTGGAAAGGAATGGTTATCgg                                >  2:4025217/1‑93 (MQ=255)
                                                           ggTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGt                               >  1:2442558/1‑93 (MQ=255)
                                                            gTCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTg                              >  1:3392543/1‑93 (MQ=255)
                                                             tCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGc                             >  1:1561496/1‑93 (MQ=255)
                                                             tCAGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGc                             <  1:1120060/93‑1 (MQ=255)
                                                              cAGGCATTCTGGCCGTTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCt                            <  2:1857072/93‑1 (MQ=255)
                                                               aGGCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCtt                           >  2:240643/1‑93 (MQ=255)
                                                                ggCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTc                          <  2:1955647/93‑1 (MQ=255)
                                                                ggCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTc                          >  1:2918319/1‑93 (MQ=255)
                                                                ggCATTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTc                          <  1:602155/93‑1 (MQ=255)
                                                                   aTTCTTGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGa                       <  1:3866565/93‑1 (MQ=255)
                                                                   aTTCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGa                       <  2:424232/93‑1 (MQ=255)
                                                                   aTTCTGGACGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGa                       >  1:4623105/1‑93 (MQ=255)
                                                                     tCTGGCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAg                     <  2:278175/93‑1 (MQ=255)
                                                                      cTGGCCGGTTCGGGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAgg                    <  1:2168609/93‑1 (MQ=255)
                                                                        ggCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGCTATCGGCGCTTCTGAAGGTa                  >  2:2678269/1‑93 (MQ=255)
                                                                         gCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTAc                 >  2:3063009/1‑93 (MQ=255)
                                                                         gCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTAc                 <  2:3520552/93‑1 (MQ=255)
                                                                         gCCGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGCATGGTTATCGGTGCTTCTGAAGGTAc                 >  2:3714824/1‑93 (MQ=255)
                                                                          ccGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACt                >  2:869012/1‑93 (MQ=255)
                                                                           cGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTg               >  2:794045/1‑93 (MQ=255)
                                                                           cGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTg               >  2:1339567/1‑93 (MQ=255)
                                                                           cGGTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTg               <  1:2856873/93‑1 (MQ=255)
                                                                            ggTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGa              >  1:4372465/1‑93 (MQ=255)
                                                                            ggTTCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGa              >  1:3692679/1‑93 (MQ=255)
                                                                              ttCGCGGGCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAg            >  1:3322219/1‑93 (MQ=255)
                                                                                   gggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTaaa       <  1:3524247/93‑1 (MQ=255)
                                                                                    ggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAg      >  1:3420310/1‑93 (MQ=255)
                                                                                    ggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAg      <  2:4569042/93‑1 (MQ=255)
                                                                                    ggCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGGATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAg      >  2:3374167/1‑93 (MQ=255)
                                                                                     gCCGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGc     >  1:4402706/1‑93 (MQ=255)
                                                                                      ccGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCg    <  1:2922832/93‑1 (MQ=255)
                                                                                       cGACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGa   <  2:1203916/93‑1 (MQ=255)
                                                                                        gACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCTAt  >  2:1637854/1‑93 (MQ=255)
                                                                                        gACACTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAt  <  2:1496405/93‑1 (MQ=255)
                                                                                           |                                                                                         
AACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGCGAACAGCTACAGGGTGAATTACGCTGGAAAGGTATGGTTATCGGTGCTTCTGAAGGTACTGAAGTTAAAGCGAT  >  NC_000913/3785291‑3785471

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