Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,609 T→C 54.8% K212K (AAA→AAG)  fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,070,6090TC54.8% 63.8 / 221.2 186K212K (AAA→AAG) fbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base T (41/43);  new base C (47/54);  total (89/97)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 188 total aligned reads displayed.

GAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCAC  >  NC_000913/3070511‑3070706
                                                                                                  |                                                                                                 
gaATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGGAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGAttt                                                                                                  <  2:732694/100‑1 (MQ=255)
 aaTCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGAtttt                                                                                                 <  1:2871927/100‑1 (MQ=255)
   tCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGAAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGc                                                                                               <  2:2900685/100‑1 (MQ=255)
    cACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCt                                                                                              >  1:1684009/1‑100 (MQ=255)
     aCGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTc                                                                                             <  1:795693/100‑1 (MQ=255)
     aCGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTc                                                                                             >  1:1841418/1‑100 (MQ=255)
     aCGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTc                                                                                             <  1:2324067/100‑1 (MQ=255)
     aCGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTc                                                                                             >  1:322924/1‑100 (MQ=255)
      cGCAGGCTGGTCGGGGTCAGACCCACGCTACCCGGCTTGTAAACTCCGTGTACGTTACCGATGGGCGGCGGGATGGTGAAACGCGGGCTGTTTTTGc                                                                                               <  2:2852179/97‑1 (MQ=255)
      cGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCa                                                                                            <  1:1440334/100‑1 (MQ=255)
         aGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGtt                                                                                         >  1:2414895/1‑100 (MQ=255)
          ggATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTc                                                                                        >  2:4033169/1‑100 (MQ=255)
          ggATGGTCGGAGTCAGAACCACGCTACCCGGCTTGTAAACCCCGTGTACGTTACCGAAGGACGCTGCGATGGTGCAACGCGGGCTGATTTTGCTCAGTTc                                                                                        >  2:943169/1‑100 (MQ=255)
          ggATGGTCGGAGGCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTc                                                                                        <  2:1858420/100‑1 (MQ=255)
            aTGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCgg                                                                                      <  1:4520/100‑1 (MQ=255)
              ggTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGgtg                                                                                    <  1:4028177/100‑1 (MQ=255)
              ggTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGgtg                                                                                    <  1:2012083/100‑1 (MQ=255)
              ggTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGgtg                                                                                    <  1:1316587/100‑1 (MQ=255)
               gTCGGGGTAAGAACCACGTTACCCGGCTTGTAAACACCGTGTGCGTTACCGAAGGACGCTCCGATGGTGAAGCGCGGGCTGACATTGATGAGTTCGgt                                                                                     >  2:2233076/1‑98 (MQ=255)
               gTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGgtgt                                                                                   >  2:695058/1‑100 (MQ=255)
                 cGGAGTCAGATCCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGgtg                                                                                    <  2:1519330/97‑1 (MQ=255)
                  ggAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATg                                                                                <  1:3819935/100‑1 (MQ=255)
                  ggAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATg                                                                                <  2:3044435/100‑1 (MQ=255)
                     gTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGt                                                                             >  1:3989638/1‑100 (MQ=255)
                      tCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGtn                                                                            <  2:4063269/100‑2 (MQ=255)
                       cAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAAAGCGGCCTGATTTTGCTGAGTTCGGTGTAGGCg                                                                              >  2:2669999/1‑97 (MQ=255)
                        aGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCCGCGCTGGTGAAACGAGGGCTGATTTTGCTCAGTTCGGTGTATGAGTa                                                                            >  2:1893708/1‑98 (MQ=255)
                         gAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAAt                                                                          >  1:2220689/1‑99 (MQ=255)
                         gAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATc                                                                         >  1:991735/1‑100 (MQ=255)
                             cACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACg                                                                     <  2:2125887/100‑1 (MQ=255)
                             cACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACg                                                                     >  2:3232836/1‑100 (MQ=255)
                              aCGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGt                                                                    >  2:2787910/1‑100 (MQ=255)
                              aCGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGt                                                                    <  2:1050374/100‑1 (MQ=255)
                              aCGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGt                                                                    <  2:229405/100‑1 (MQ=255)
                                 ttACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGAAGCTGCGATGGTGAAACGCGGGCTAATCTTGCTAAGCACGGCCTATACGTAATCTACGTCt                                                                  >  1:761747/1‑99 (MQ=255)
                                  tACCCGGCTTGTAAACACGGTGTACGTTACTGAAGGGCGCTGCGATGGCGAAAAGCGGGCTGACCTTGCTCAGGCCAGTGTATGCGTAATCGACGGCTTc                                                                >  2:3278807/1‑100 (MQ=255)
                                  tACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTc                                                                >  1:3159022/1‑100 (MQ=255)
                                   aCCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTcc                                                               >  1:595531/1‑100 (MQ=255)
                                    cccGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCg                                                              >  1:3440149/1‑100 (MQ=255)
                                     ccGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCCGTTCGGTGTATGCGTAATCAACGTCTTCCgg                                                             >  2:1627691/1‑100 (MQ=255)
                                      cGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGc                                                            <  1:4127512/100‑1 (MQ=255)
                                       ggCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCt                                                           >  2:3358704/1‑100 (MQ=255)
                                         cTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTgg                                                         <  2:3440129/100‑1 (MQ=255)
                                          ttGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGTGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTggg                                                        <  2:2157615/100‑1 (MQ=255)
                                          ttGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTggg                                                        <  1:1523743/100‑1 (MQ=255)
                                          ttGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTggg                                                        <  1:2933706/100‑1 (MQ=255)
                                            gTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCGCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGgtg                                                      <  2:1444679/100‑1 (MQ=255)
                                              aaaCACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTa                                                    >  2:2688488/1‑100 (MQ=255)
                                               aaCACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAGCGTCTTCCGGCTGGGTGTAc                                                   <  2:2524829/100‑1 (MQ=255)
                                               aaCACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGgtgt                                                     >  2:383335/1‑98 (MQ=255)
                                                acacCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACa                                                  <  2:3474205/100‑1 (MQ=255)
                                                  acCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGt                                                <  1:1862445/100‑1 (MQ=255)
                                                   ccGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTg                                               <  1:3488830/100‑1 (MQ=255)
                                                   ccGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGGTCGGTGTATGAGTAATCAACGTCTTCCGGCTGGGTGGACAGt                                                <  2:3419059/99‑1 (MQ=255)
                                                    cGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCTGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGc                                              <  1:191429/100‑1 (MQ=255)
                                                    cGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGc                                              <  1:3956240/100‑1 (MQ=255)
                                                    cGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGc                                              <  2:2763018/100‑1 (MQ=255)
                                                       gtACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGa                                           <  2:247878/100‑1 (MQ=255)
                                                        tACGGTACCGAAAGACGCTGCGAGGGTGAAACGCGGGCGGCTTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCTGCTAGGTGTACAGTg                                               <  2:1426036/95‑1 (MQ=255)
                                                          ctttACCGAAGGGCGCTGCGATGGTTAGACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGaa                                          <  2:2433603/96‑1 (MQ=255)
                                                          cGTTACCGCAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACATTGCAGAAGc                                        >  2:604293/1‑100 (MQ=255)
                                                           gTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCg                                       >  1:1122711/1‑100 (MQ=255)
                                                           gTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCg                                       >  1:2166206/1‑100 (MQ=255)
                                                              aCCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTcc                                    <  2:3387170/100‑1 (MQ=255)
                                                              aCCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCCCATTTCCTTGTACGCGCACCCCACGCCTTCACGCCGGGTGCACATTGCAGAAGcacc                                     >  2:3037234/1‑96 (MQ=255)
                                                               ccGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCa                                   >  1:2472265/1‑100 (MQ=255)
                                                               ccGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCa                                   >  1:275362/1‑100 (MQ=255)
                                                               ccGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCa                                   >  1:1289228/1‑100 (MQ=255)
                                                               ccGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCa                                   >  2:802556/1‑100 (MQ=255)
                                                               ccGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCa                                   <  2:2871078/100‑1 (MQ=255)
                                                                cGAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAt                                  >  1:1448421/1‑100 (MQ=255)
                                                                 gAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAtg                                 >  1:3100365/1‑100 (MQ=255)
                                                                 gAAGGACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAtg                                 >  2:593088/1‑100 (MQ=255)
                                                                  aaGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAtgt                                <  1:1564028/100‑1 (MQ=255)
                                                                     gACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGc                             <  1:2588049/100‑1 (MQ=255)
                                                                     gACGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGc                             <  1:2147706/100‑1 (MQ=255)
                                                                      aCGCTGCGATGGTGAAACGCGGGCTGATTCTGCTCAGTTCGGCGTATGCGTAAGCAACGTCTTCCGGCTGGGTGTACAGTGCAGAACCGTCCATGTTGCt                            >  2:3307383/1‑100 (MQ=255)
                                                                      aCGCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCt                            <  2:1334827/100‑1 (MQ=255)
                                                                        gCTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtgt                          <  2:843513/100‑1 (MQ=255)
                                                                         cTGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtgtt                         <  1:3523395/100‑1 (MQ=255)
                                                                          tGCGATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtgttg                        <  1:2956801/100‑1 (MQ=255)
                                                                             gATGGTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTGCATGTGGCTGTTGTcc                     >  2:960476/1‑100 (MQ=255)
                                                                                ggTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACg                  >  2:1727059/1‑100 (MQ=255)
                                                                                 gTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGc                 <  2:450554/100‑1 (MQ=255)
                                                                                 gTGAAACGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGc                 >  2:103005/1‑100 (MQ=255)
                                                                                   gAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCg               >  1:3939396/1‑100 (MQ=255)
                                                                                   gAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCg               >  2:2659351/1‑100 (MQ=255)
                                                                                    aaaCGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGt              >  1:2529023/1‑100 (MQ=255)
                                                                                     aaCGCGGGCTTATCTTGCGCATTTCGGGGTATGCGTAATCGACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGATGTTGTCCACGCCg               >  2:2408778/1‑98 (MQ=255)
                                                                                     aaCGCGGGCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGtc             >  2:2980654/1‑100 (MQ=255)
                                                                                           ggCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCttcttc       >  2:3111696/1‑100 (MQ=255)
                                                                                           ggCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCttcttc       <  1:1575054/100‑1 (MQ=255)
                                                                                           ggCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCttcttc       <  2:388614/100‑1 (MQ=255)
                                                                                           ggCTGATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCttcttc       <  1:121516/100‑1 (MQ=255)
                                                                                            gCTGTTCTTCCTCAGTACGGTGTATGCGTAATGAACGTCTTCAGGCTGAGGGTACAGTGCAGAAGCGTCCATGTGGCTGTCGTCCACGCCGTCttcttc       <  2:2674996/99‑1 (MQ=255)
                                                                                            gCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTca      <  1:3671762/100‑1 (MQ=255)
                                                                                            gCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTca      <  1:1335632/100‑1 (MQ=255)
                                                                                             cTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTcac     <  1:3497886/100‑1 (MQ=255)
                                                                                               gATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTcacca   >  2:923187/1‑100 (MQ=255)
                                                                                               gATCTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTcacca   >  1:1144570/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
GAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCAC  >  NC_000913/3070511‑3070706

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑