Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,184,484 A→G F112L (TTT→CTT)  potA ← ATP‑binding component of spermidine/putrescine transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,184,4840AG99.3% 460.1 / ‑4.8 143F112L (TTT→CTT) potAATP‑binding component of spermidine/putrescine transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (58/84);  total (58/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 148 total aligned reads displayed.

GCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGG  >  NC_000913/1184396‑1184573
                                                                                        |                                                                                         
gCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGcc                                                                                       <  1:4646762/93‑1 (MQ=255)
  gAAGGTTTCCAACTGCACCATCCGCAGGGCGTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCcac                                                                                     >  2:2194695/1‑93 (MQ=255)
  gAAGGTTTCCAACTGCACCATCCGCAGGGAATCCATCACGCGGGGCGTAATTTCAGCAACGGGGGTTTTTTGCATGCGCCACCCAAGGGGcac                                                                                     >  1:2283045/1‑93 (MQ=255)
   aaGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGTTTTTTTGCATGCGCAACCCAAGGGCcaca                                                                                    <  1:2397055/93‑1 (MQ=255)
     ggTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACAtt                                                                                  <  1:1280175/93‑1 (MQ=255)
      gTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACAttt                                                                                 <  1:3206390/93‑1 (MQ=255)
        ttCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTc                                                                               <  2:3761108/93‑1 (MQ=255)
        ttCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTc                                                                               <  2:2692427/93‑1 (MQ=255)
         tCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCg                                                                              <  2:4392469/93‑1 (MQ=255)
          ccACCTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGa                                                                             <  2:2992836/93‑1 (MQ=255)
           cAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGaa                                                                            >  2:2167767/1‑93 (MQ=255)
           cAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGaa                                                                            >  2:161163/1‑93 (MQ=255)
            aaCTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAac                                                                           >  1:4103624/1‑93 (MQ=255)
              cTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAacac                                                                         <  2:2634339/93‑1 (MQ=255)
              cTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAacac                                                                         <  1:2806442/93‑1 (MQ=255)
               tGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACg                                                                        >  1:1818755/1‑93 (MQ=255)
               tGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACg                                                                        <  1:2120833/93‑1 (MQ=255)
                gCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACgg                                                                       <  1:1367267/93‑1 (MQ=255)
                 cACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACggc                                                                      >  2:1863129/1‑92 (MQ=255)
                  aCCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTTGAACACGGTc                                                                     >  1:3037006/1‑93 (MQ=255)
                  aCCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTc                                                                     >  1:4319362/1‑93 (MQ=255)
                      tCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCAtgg                                                                 >  1:42701/1‑92 (MQ=255)
                      tCCGCAGCGATTCCATCACGCGAGGCGTAATTTCATCAGCGGTGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCAtgt                                                                 <  1:195194/93‑1 (MQ=255)
                       ccGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCAtgtg                                                                <  1:883332/93‑1 (MQ=255)
                           aGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTggggg                                                            <  1:1172353/93‑1 (MQ=255)
                           aGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTggggg                                                            <  2:3533423/93‑1 (MQ=255)
                            gggCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGa                                                           >  2:3838550/1‑93 (MQ=255)
                            gggCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGa                                                           >  2:3838551/1‑93 (MQ=255)
                                ttCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAg                                                       >  2:2794225/1‑93 (MQ=255)
                                 tCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGt                                                      >  2:3378539/1‑93 (MQ=255)
                                  ccATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTg                                                     >  2:3300082/1‑93 (MQ=255)
                                  ccATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTg                                                     <  1:2436541/93‑1 (MQ=255)
                                  ccATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTg                                                     <  2:1787519/93‑1 (MQ=255)
                                  ccATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTg                                                     <  2:4412084/93‑1 (MQ=255)
                                   cATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGc                                                    >  2:504894/1‑93 (MQ=255)
                                    aTCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCg                                                   >  1:3321370/1‑93 (MQ=255)
                                      cACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTa                                                 <  2:3414458/93‑1 (MQ=255)
                                      cACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTa                                                 <  1:2556663/93‑1 (MQ=255)
                                      cACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTa                                                 <  2:709439/93‑1 (MQ=255)
                                        cgcgGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGc                                               <  1:4131353/93‑1 (MQ=255)
                                         gcgGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCt                                              >  1:1200566/1‑93 (MQ=255)
                                          cgGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCtt                                             <  2:2457573/93‑1 (MQ=255)
                                           ggggCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCttt                                            <  1:2624017/93‑1 (MQ=255)
                                           ggggCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGCCATGTGGGGGAAAAGGGCTTAGtttg                                            >  1:1102355/1‑92 (MQ=255)
                                            gggCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTg                                           <  1:4649637/93‑1 (MQ=255)
                                            gggCGTAATTTCAGCAGCGGGGGCTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTg                                           >  2:4296151/1‑93 (MQ=255)
                                              gCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGa                                         <  1:4658432/93‑1 (MQ=255)
                                                gTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaaa                                       >  1:897772/1‑93 (MQ=255)
                                                 tAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaaaa                                      <  1:3981230/93‑1 (MQ=255)
                                                 tAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGaaaa                                      <  1:2213012/93‑1 (MQ=255)
                                                 tAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGCCATGTGGGGGAAAAGGGCGTACCTTTGGaaaa                                      >  1:1436367/1‑93 (MQ=255)
                                                  aaTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGTAAAAc                                     >  1:725191/1‑93 (MQ=255)
                                                  aaTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAAc                                     <  2:916313/93‑1 (MQ=255)
                                                  aaTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAAc                                     <  1:3627319/93‑1 (MQ=255)
                                                  aaTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAAc                                     <  1:1915804/93‑1 (MQ=255)
                                                   aTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACa                                    <  1:4430984/93‑1 (MQ=255)
                                                   aTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACa                                    <  1:4672415/93‑1 (MQ=255)
                                                   aTTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACa                                    <  2:1214477/93‑1 (MQ=255)
                                                     ttCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgt                                  <  2:3150669/93‑1 (MQ=255)
                                                     ttCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgt                                  <  2:133699/93‑1 (MQ=255)
                                                     ttCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgt                                  <  1:1549260/93‑1 (MQ=255)
                                                     ttCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgt                                  <  1:512576/93‑1 (MQ=255)
                                                      tCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgtg                                 <  2:1988740/93‑1 (MQ=255)
                                                       cagcagCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgtgt                                <  2:1788421/93‑1 (MQ=255)
                                                       cagcagCGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAgtgt                                <  1:2170713/93‑1 (MQ=255)
                                                        agcagcGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGtt                               <  2:4643034/93‑1 (MQ=255)
                                                        agcagcGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGtt                               >  2:3206253/1‑93 (MQ=255)
                                                         gcagcGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTc                              <  1:152259/93‑1 (MQ=255)
                                                         gcagcGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTc                              <  1:1803147/93‑1 (MQ=255)
                                                         gcagcGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTc                              >  1:1259686/1‑93 (MQ=255)
                                                             cGGGGGTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACAt                          <  1:2003600/93‑1 (MQ=255)
                                                              gggggTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATa                         <  1:43879/93‑1 (MQ=255)
                                                              gggggTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATa                         >  1:3670923/1‑93 (MQ=255)
                                                               ggggTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAg                        >  1:883862/1‑93 (MQ=255)
                                                               ggggTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAg                        >  1:2049993/1‑93 (MQ=255)
                                                               ggggTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAg                        >  2:2921767/1‑93 (MQ=255)
                                                                  gTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCgg                     <  2:546980/93‑1 (MQ=255)
                                                                  gTTTTTTGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCgg                     >  1:3200384/1‑93 (MQ=255)
                                                                   ttttttGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGt                    <  2:4744536/93‑1 (MQ=255)
                                                                   tttgttGCATGCGCAGCCCGAGGGCCACATTTTCGAACAGGGTCATGTGGGGGACAAGTGCGTAGATTTGGAAAACAGTGTTCACATAggggg                    <  2:1668362/93‑2 (MQ=255)
                                                                     ttttGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGttt                  >  2:1083395/1‑93 (MQ=255)
                                                                     ttttGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGttt                  <  1:4228643/93‑1 (MQ=255)
                                                                      tttGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGtttt                 <  1:2599163/93‑1 (MQ=255)
                                                                      tttGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATACCGGttgt                 >  1:904154/1‑91 (MQ=255)
                                                                       ttGCATGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTc                <  1:3765199/93‑1 (MQ=255)
                                                                           aTGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTccgcc            <  2:86599/93‑1 (MQ=255)
                                                                           aTGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTccgcc            >  1:288757/1‑93 (MQ=255)
                                                                            tGCGCAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTccgccg           <  2:4497682/93‑1 (MQ=255)
                                                                              cgcAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGa         <  1:3151936/93‑1 (MQ=255)
                                                                               gcAAGCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGaa        >  1:4698753/1‑93 (MQ=255)
                                                                               gcAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGaa        >  2:1102956/1‑93 (MQ=255)
                                                                               gcAACCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGaa        >  2:352368/1‑93 (MQ=255)
                                                                                 aaCCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGTAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACg      >  1:1637322/1‑93 (MQ=255)
                                                                                 aaCCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGTAAAAGAGCGTAGCTTTGGAAACCAGTGTACGCCTAGAGAGGTTCCGCCGGCACg      >  1:1637321/1‑93 (MQ=255)
                                                                                 aaCCCAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACg      >  2:4141396/1‑93 (MQ=255)
                                                                                   cccAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTg    >  1:3391133/1‑93 (MQ=255)
                                                                                   cccAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTg    <  1:755992/93‑1 (MQ=255)
                                                                                    ccAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTgg   <  2:3674955/93‑1 (MQ=255)
                                                                                    ccAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTgg   >  1:51811/1‑93 (MQ=255)
                                                                                     cAAGGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTggg  >  1:4292720/1‑93 (MQ=255)
                                                                                        |                                                                                         
GCGAAGGTTTCCAACTGCACCATCCGCAGGGCTTCCATCACGCGGGGCGTAATTTCAGCAGCGGGGGTTTTTTGCATGCGCAACCCAAAGGCCACATTTTCGAACACGGTCATGTGGGGGAAAAGTGCGTAGCTTTGGAAAACAGTGTTCACATAGCGGTTTTCCGCCGGAACGTGGG  >  NC_000913/1184396‑1184573

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