Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,493,645 G→A intergenic (‑391/+122) frc ← / ← yfdX formyl‑CoA transferase, NAD(P)‑binding/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,493,6450GA95.6% 296.5 / 0.7 90intergenic (‑391/+122)frc/yfdXformyl‑CoA transferase, NAD(P)‑binding/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (45/41);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCTT  >  NC_000913/2493548‑2493741
                                                                                                 |                                                                                                
tATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACg                                                                                                <  1:1797548/100‑1 (MQ=255)
 aTGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGa                                                                                               >  1:406915/1‑100 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAg                                                                                              <  2:686519/100‑1 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                              >  1:1549536/1‑100 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                              >  1:1578611/1‑100 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                              <  2:70646/100‑1 (MQ=255)
   gATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGt                                                                                             >  1:1218971/1‑100 (MQ=255)
      tACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTta                                                                                          <  2:974627/100‑1 (MQ=255)
       aCGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtat                                                                                         <  2:1661777/100‑1 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                        <  1:414595/100‑1 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                        >  2:615620/1‑100 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                        <  2:851421/100‑1 (MQ=255)
         gCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                        <  1:1301943/99‑1 (MQ=255)
         gCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAg                                                                                       <  1:834660/100‑1 (MQ=255)
          cTTTTTAAAGCGTACCATAACACATAACATAATAAATGATTTGTGATGACACGGGATTCTTGCGGTTTACTTTGACGAATATAAAATACGAGTTTATAGc                                                                                      >  1:453628/1‑100 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                      >  2:1034280/1‑100 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                      >  1:657284/1‑100 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                      <  2:967361/100‑1 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                      >  2:1391422/1‑100 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                      >  2:188776/1‑100 (MQ=255)
               tAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACcgc                                                                                 <  1:781599/100‑1 (MQ=255)
                aaaGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACcgcg                                                                                >  1:1230087/1‑100 (MQ=255)
                  aGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAg                                                                              >  1:1753344/1‑100 (MQ=255)
                   gCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAgc                                                                             >  1:1272254/1‑100 (MQ=255)
                    cGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAgcg                                                                            <  1:1701339/100‑1 (MQ=255)
                     gTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAgcgc                                                                           >  1:859616/1‑100 (MQ=255)
                      tACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGcc                                                                          <  1:382041/100‑1 (MQ=255)
                      tACCATAAAAAATAACATAATAAATGATTTGTGATGACACGGGATTATTGCTGTTTACTTTGATGAATATAAAATACGAGTTTATAGCACCGCGAgcg                                                                            <  2:1520341/98‑1 (MQ=255)
                       aCCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAgcgc                                                                           <  1:1615255/98‑1 (MQ=255)
                       aCCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGccn                                                                         <  1:699172/100‑2 (MQ=255)
                        ccATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCaa                                                                        >  2:1476837/1‑100 (MQ=255)
                         cATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCaaa                                                                       >  1:249815/1‑100 (MQ=255)
                         cATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCaaa                                                                       >  2:1251266/1‑100 (MQ=255)
                          aTAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaat                                                                      >  1:1625148/1‑100 (MQ=255)
                             aaaCATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaataat                                                                   >  1:198070/1‑100 (MQ=255)
                              aacataacatAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAAtt                                                                  >  1:375338/1‑100 (MQ=255)
                                cataacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAAtttt                                                                <  2:1778328/100‑1 (MQ=255)
                                 ataacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAAtttt                                                                <  2:175014/99‑1 (MQ=255)
                                 ataacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTa                                                               >  2:1424563/1‑100 (MQ=255)
                                  taacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTAt                                                              <  1:349045/100‑1 (MQ=255)
                                   aacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                             <  1:261265/100‑1 (MQ=255)
                                   aacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                             <  1:96986/100‑1 (MQ=255)
                                    acataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGa                                                            >  1:1206203/1‑100 (MQ=255)
                                     cataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGAc                                                           >  1:1259354/1‑100 (MQ=255)
                                     cataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGAc                                                           <  2:1125142/100‑1 (MQ=255)
                                       taataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACca                                                         <  1:1423151/100‑1 (MQ=255)
                                       taataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACca                                                         >  2:1022131/1‑100 (MQ=255)
                                        aataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACcac                                                        >  1:750653/1‑100 (MQ=255)
                                         ataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACcaca                                                       <  2:1705693/100‑1 (MQ=255)
                                            aaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                    <  2:1596677/100‑1 (MQ=255)
                                            aaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                    <  2:1768569/100‑1 (MQ=255)
                                            aaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                    <  1:188775/100‑1 (MQ=255)
                                               gATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGt                                                 >  1:91461/1‑100 (MQ=255)
                                                ntttGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTc                                                >  1:955054/2‑100 (MQ=255)
                                                aTTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTc                                                >  2:1009639/1‑100 (MQ=255)
                                                      gatgaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtg                                          >  2:680008/1‑100 (MQ=255)
                                                      gatgaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtg                                          <  1:471361/100‑1 (MQ=255)
                                                      gatgaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtg                                          >  1:1251663/1‑100 (MQ=255)
                                                        tgaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgtg                                        >  1:609684/1‑100 (MQ=255)
                                                         gaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTgg                                       <  2:1031312/100‑1 (MQ=255)
                                                         gaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGAGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgt                                         >  1:882481/1‑98 (MQ=255)
                                                          acacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt                                      <  2:901661/100‑1 (MQ=255)
                                                           cacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGATCACAACCTGTCGTATTGTGGt                                      >  1:707868/1‑99 (MQ=255)
                                                           cacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTc                                     <  2:406917/100‑1 (MQ=255)
                                                           cacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTc                                     >  2:1615109/1‑100 (MQ=255)
                                                             cGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCAt                                   <  1:1142111/100‑1 (MQ=255)
                                                              gggATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATa                                  <  2:859617/100‑1 (MQ=255)
                                                                 attattGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATaaaa                               <  2:1195586/100‑1 (MQ=255)
                                                                  ttattGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGATAACAACCTGTCGTATTGTGGTCATaaaaa                              >  1:468354/1‑100 (MQ=255)
                                                                    attGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAAt                            >  2:153361/1‑100 (MQ=255)
                                                                    attGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGCCATAAAAAAt                            <  1:1423891/100‑1 (MQ=255)
                                                                    attGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATGGTGGTCATAAAAAAt                            >  2:575719/1‑100 (MQ=255)
                                                                       gCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGt                         >  2:937540/1‑100 (MQ=255)
                                                                        cGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTACTGTGGTCATAAAAAATTGTc                        <  1:1839323/100‑1 (MQ=255)
                                                                            ttACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCt                    >  2:722107/1‑100 (MQ=255)
                                                                             tACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTg                   <  2:778800/100‑1 (MQ=255)
                                                                              aCTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgg                  >  2:496776/1‑100 (MQ=255)
                                                                               cTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgga                 >  1:119823/1‑100 (MQ=255)
                                                                                tttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggag                <  1:740453/100‑1 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggg              <  1:1424555/100‑1 (MQ=255)
                                                                                    aGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagggag            >  1:1374762/1‑100 (MQ=255)
                                                                                       aaTATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTa         >  1:44395/1‑100 (MQ=255)
                                                                                       aaTATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTa         <  2:1730946/100‑1 (MQ=255)
                                                                                           taAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATc     >  1:1787155/1‑100 (MQ=255)
                                                                                           taAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATc     >  2:234281/1‑100 (MQ=255)
                                                                                            aaaaTGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTAGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc    <  1:815359/100‑1 (MQ=255)
                                                                                            aaaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc    <  2:1578622/100‑1 (MQ=255)
                                                                                            aaaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc    <  2:1549546/100‑1 (MQ=255)
                                                                                            aaaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc    <  1:680007/100‑1 (MQ=255)
                                                                                              aaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCtt  >  2:726762/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCTT  >  NC_000913/2493548‑2493741

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