Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,338,141 (C)5→6 coding (321/552 nt) mlaD ← intermembrane phospholipid transport system, substrate binding protein MlaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,338,1361.C99.6% 1048.2 / ‑5.0 249E109G (GAA→GGA) mlaDintermembrane phospholipid transport system, substrate binding protein MlaD
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (107/141);  minor base A (1/0);  total (108/141)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 282 total aligned reads displayed.

CACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAG  >  NC_000913/3338036‑3338231
                                                                                                     |                                                                                               
cACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGtt                                                                                                  <  1:4891065/101‑1 (MQ=255)
cACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGtt                                                                                                  <  2:4890864/101‑1 (MQ=255)
   cATCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTccc                                                                                              <  2:4891100/101‑1 (MQ=255)
    aTCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTcccc                                                                                             <  1:4891333/101‑1 (MQ=255)
    aTCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTcccc                                                                                             >  1:4891627/1‑101 (MQ=255)
      cGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTcccccc                                                                                           >  1:4891630/1‑100 (MQ=255)
       gCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTcccccca                                                                                          >  2:4891632/1‑99 (MQ=255)
        cAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTccccccag                                                                                         <  2:4891036/101‑4 (MQ=255)
          gACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTccccccagca                                                                                       >  2:4891636/1‑96 (MQ=255)
          gACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTccccccagca                                                                                       >  2:4891639/1‑96 (MQ=255)
           aCTTATTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCcagcag                                                                                       <  2:4890999/101‑1 (MQ=255)
           aCTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCcagcag                                                                                       <  2:4891262/101‑1 (MQ=255)
           aCTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCcagcag                                                                                       <  1:4891172/101‑1 (MQ=255)
           aCTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCcagcag                                                                                       >  1:4891640/1‑101 (MQ=255)
            cTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAgg                                                                                      <  1:4891505/101‑1 (MQ=255)
              tAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGcc                                                                                    >  1:4891644/1‑101 (MQ=255)
               aGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGAAGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCg                                                                                   >  2:4891647/1‑101 (MQ=255)
                 tgtCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGa                                                                                 >  1:4891649/1‑101 (MQ=255)
                   tCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAg                                                                               <  1:4891087/101‑1 (MQ=255)
                    ccTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGt                                                                              >  1:4891652/1‑101 (MQ=255)
                        aaTTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGa                                                                          >  2:4891656/1‑101 (MQ=255)
                        aaTTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGa                                                                          <  1:4891204/101‑1 (MQ=255)
                         aTTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGaa                                                                         <  2:4891252/101‑1 (MQ=255)
                          ttGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAAt                                                                        <  2:4891083/101‑1 (MQ=255)
                            gTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCgg                                                                                  >  2:4891660/1‑89 (MQ=255)
                              aTCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTc                                                                    >  1:4891664/1‑101 (MQ=255)
                               tCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCa                                                                   <  2:4891090/101‑1 (MQ=255)
                                 gCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGc                                                                 <  2:4891113/101‑1 (MQ=255)
                                    aTCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGaa                                                              >  1:4891668/1‑101 (MQ=255)
                                    aTCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGaa                                                              >  2:4891667/1‑101 (MQ=255)
                                     tCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAAc                                                             >  2:4891671/1‑101 (MQ=255)
                                       cTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTg                                                           <  2:4891295/101‑1 (MQ=255)
                                         tCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGt                                                         <  2:4891240/101‑1 (MQ=255)
                                          cAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTa                                                        >  1:4891674/1‑101 (MQ=255)
                                           aGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTAt                                                       <  1:4890972/101‑1 (MQ=255)
                                             gATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCt                                                     <  2:4891281/101‑1 (MQ=255)
                                             gATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCt                                                     <  1:4891155/101‑1 (MQ=255)
                                              aTAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTg                                                    <  1:4891314/101‑1 (MQ=255)
                                              aTAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTg                                                    <  2:4891276/101‑1 (MQ=255)
                                               tAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTgg                                                   <  2:4891492/101‑1 (MQ=255)
                                                aGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGa                                                  <  1:4891223/101‑1 (MQ=255)
                                                  cAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAAt                                                <  2:4891177/101‑1 (MQ=255)
                                                  cAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTACCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAAt                                                >  1:4891678/1‑101 (MQ=255)
                                                   aGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAAtg                                               >  2:4891680/1‑101 (MQ=255)
                                                    gTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAAtgt                                              >  1:4891683/1‑101 (MQ=255)
                                                        ccAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCg                                                                >  1:4891688/1‑79 (MQ=255)
                                                         cAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTta                                         >  1:4891690/1‑101 (MQ=255)
                                                         cAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTta                                         <  1:4891146/101‑1 (MQ=255)
                                                          aGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTtat                                        >  2:4891692/1‑101 (MQ=255)
                                                          aGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTtat                                        <  1:4891621/101‑1 (MQ=255)
                                                          aGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTtat                                        >  2:4891695/1‑101 (MQ=255)
                                                           gTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTtata                                       >  1:4891697/1‑101 (MQ=255)
                                                           gTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTtata                                       <  2:4891145/101‑1 (MQ=255)
                                                            ttCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATaa                                      >  1:4891700/1‑101 (MQ=255)
                                                               cGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACgt                                   >  2:4891703/1‑101 (MQ=255)
                                                                gggTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACgtt                                  <  1:4891260/101‑1 (MQ=255)
                                                                 ggTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACgttg                                 <  1:4891289/101‑1 (MQ=255)
                                                                 ggTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACgttg                                 >  2:4891708/1‑101 (MQ=255)
                                                                  gTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACgttgt                                <  1:4891336/101‑1 (MQ=255)
                                                                  gTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACgttgt                                >  1:4891709/1‑101 (MQ=255)
                                                                    cTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTc                              <  2:4891415/101‑1 (MQ=255)
                                                                     ttCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAAc                                     <  1:4891447/93‑1 (MQ=255)
                                                                     ttCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCa                             <  1:4891414/101‑1 (MQ=255)
                                                                      tCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCaa                            <  1:4891400/101‑1 (MQ=255)
                                                                      tCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCaa                            >  1:4891712/1‑101 (MQ=255)
                                                                       cAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAAt                           >  2:4891714/1‑101 (MQ=255)
                                                                       cAAAACCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAAt                           >  1:4891717/1‑101 (MQ=255)
                                                                        aaaaCCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAAtt                          >  1:4891719/1‑101 (MQ=255)
                                                                          aaCCGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTc                        <  2:4891190/101‑1 (MQ=255)
                                                                            ccGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCa                      <  2:4891360/101‑1 (MQ=255)
                                                                            ccGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCa                      >  1:4891724/1‑101 (MQ=255)
                                                                            ccGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCa                      <  2:4891224/101‑1 (MQ=255)
                                                                            ccGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCa                      <  1:4891365/101‑1 (MQ=255)
                                                                             cGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCGGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAg                     <  1:4891008/101‑1 (MQ=255)
                                                                             cGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAg                     <  1:4891249/101‑1 (MQ=255)
                                                                             cGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAg                     <  1:4891378/101‑1 (MQ=255)
                                                                             cGACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAg                     <  1:4891467/101‑1 (MQ=255)
                                                                              gACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGc                    <  1:4891234/101‑1 (MQ=255)
                                                                              gACGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGCATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGc                    >  2:4891725/1‑101 (MQ=255)
                                                                               aCGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCg                   <  1:4891161/101‑1 (MQ=255)
                                                                                cGTTTAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGt                  >  2:4891728/1‑101 (MQ=255)
                                                                                  tttAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTa                >  1:4891733/1‑101 (MQ=255)
                                                                                  tttAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTa                <  1:4891424/101‑1 (MQ=255)
                                                                                    tAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcg              <  2:4891474/101‑1 (MQ=255)
                                                                                    tAATGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcg              <  2:4891122/101‑1 (MQ=255)
                                                                                     aaTGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgc             >  1:4891736/1‑101 (MQ=255)
                                                                                     aaTGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgc             >  2:4891738/1‑101 (MQ=255)
                                                                                     aaTGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgc             <  2:4891304/101‑1 (MQ=255)
                                                                                      aTGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgcg            >  1:4891741/1‑101 (MQ=255)
                                                                                      aTGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgcg            <  2:4891339/101‑1 (MQ=255)
                                                                                      aTGCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgcg            <  2:4891344/101‑1 (MQ=255)
                                                                                        gCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgcgcg          >  1:4891744/1‑101 (MQ=255)
                                                                                        gCCAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTAcgcgcg          <  2:4891388/101‑1 (MQ=255)
                                                                                         ccAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCgg         <  2:4891214/101‑1 (MQ=255)
                                                                                          cAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGc        >  1:4891747/1‑101 (MQ=255)
                                                                                          cAGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGc        <  1:4891273/101‑1 (MQ=255)
                                                                                           aGATATTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCa       >  2:4891749/1‑101 (MQ=255)
                                                                                             atatTGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGa     <  2:4891312/101‑1 (MQ=255)
                                                                                                ttGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAg  <  2:4891540/101‑1 (MQ=255)
                                                                                                ttGTTCCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAg  <  2:4891350/101‑1 (MQ=255)
                                                                                                     |                                                                                               
CACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCATCCTTCAGGATAGCAGTCCCCAGTTCCGGGTCTTCAAAACCGACGTTTAATGCCAGATATTGTTCCCCCAGCAGGCCGGAAGTACGAATGCTCAGCGAACTGGTATCTGGAATGTGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAG  >  NC_000913/3338036‑3338231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: