Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,493,645 G→A 100% 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,6450GA89.3% 287.4 / 22.5 112intergenic (‑391/+122)frc/yfdXformyl‑CoA transferase, NAD(P)‑binding/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (9/3);  new base A (48/52);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 112 total aligned reads displayed.

TATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCC  >  NC_000913/2493548‑2493739
                                                                                                 |                                                                                              
ttatggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTATGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGc                                                                                               >  2:23156/7‑99 (MQ=255)
 aTGAATACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGa                                                                                             <  2:413672/100‑1 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAg                                                                                            >  2:1052088/1‑100 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                            >  2:1125223/1‑100 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                            <  1:641184/100‑1 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                            <  1:1127202/100‑1 (MQ=255)
  atggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                            >  1:1778125/5‑100 (MQ=255)
   tggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGt                                                                                           >  1:1973268/4‑100 (MQ=255)
   tggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGt                                                                                           >  2:2106585/4‑100 (MQ=255)
   tggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGt                                                                                           >  2:898081/4‑100 (MQ=255)
    aTTACGCTTTTTAGAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGtt                                                                                          >  1:1743316/1‑100 (MQ=255)
    aTTACGCTTTTTAAATCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACTGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGtt                                                                                          >  2:303481/1‑100 (MQ=255)
      tACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTta                                                                                        <  1:1250345/100‑1 (MQ=255)
       aCGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtat                                                                                       <  2:590490/100‑1 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                      >  1:2275772/1‑100 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                      >  1:2366105/1‑100 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                      >  1:2176914/1‑100 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                      <  2:939682/100‑1 (MQ=255)
         gatttttAAAGCGTACGAGAAAACATATCATAATAAAAGATTTGGGAAGACTCGGGATTAGTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                      <  2:271936/97‑1 (MQ=255)
             tttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCAcc                                                                                 <  2:1955256/100‑1 (MQ=255)
             tttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCAcc                                                                                 >  1:310643/1‑100 (MQ=255)
              ttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACcg                                                                                >  2:1620007/1‑100 (MQ=255)
               tAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACcgc                                                                               <  2:146451/100‑1 (MQ=255)
                aaaGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACcgcg                                                                              >  1:2307074/1‑100 (MQ=255)
                 aaGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGa                                                                             >  2:707409/1‑100 (MQ=255)
                    cGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAgcg                                                                          >  1:399438/1‑100 (MQ=255)
                       aCCATAAAACATAACATAATAAATGATTTGTGATGACACAGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAgcgac                                                                        >  2:192097/1‑97 (MQ=255)
                        ccATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCaa                                                                      <  2:2451989/100‑1 (MQ=255)
                        ccATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCaa                                                                      >  2:1037778/1‑100 (MQ=255)
                         cATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCaaa                                                                     <  1:1930444/100‑1 (MQ=255)
                          aTAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaat                                                                    <  2:2438052/100‑1 (MQ=255)
                          aTAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTATACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaat                                                                    >  2:1006799/1‑100 (MQ=255)
                            aaaaCATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGATTAATTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaataa                                                                  >  2:527185/1‑100 (MQ=255)
                             aaaCATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaataat                                                                 >  2:102891/1‑100 (MQ=255)
                               acataacataATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGAGAGCGCCAAATAAttt                                                               <  2:117394/100‑1 (MQ=255)
                                cataacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAAtttt                                                              >  2:1979695/1‑100 (MQ=255)
                                  taacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTAt                                                            <  2:671819/100‑1 (MQ=255)
                                   aacataaTAAATGATTTGTGATGACACGGTATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTAt                                                            <  2:1200360/99‑1 (MQ=255)
                                   aacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                           <  2:763973/100‑1 (MQ=255)
                                   aacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                           <  1:1787804/100‑1 (MQ=255)
                                   aacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                           <  2:610098/100‑1 (MQ=255)
                                    acataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCACCAAATAATTTTATGa                                                          <  1:2309829/100‑1 (MQ=255)
                                      ataataAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGAcc                                                        <  2:100288/100‑1 (MQ=255)
                                         ataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACcaca                                                     >  1:1322681/1‑100 (MQ=255)
                                           aaaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACGGCGAGCGCCAAATAATTTTATGAcaacaa                                                    >  1:2017994/1‑99 (MQ=255)
                                           aaaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAc                                                   <  2:1742068/100‑1 (MQ=255)
                                            aaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                  <  1:1037230/100‑1 (MQ=255)
                                                aTTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTAGAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTc                                              >  2:10298/1‑100 (MQ=255)
                                                 tttGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCg                                             <  2:1708498/100‑1 (MQ=255)
                                                 tttGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCg                                             <  1:2020789/100‑1 (MQ=255)
                                                   tgtgATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTa                                           >  1:1603167/1‑100 (MQ=255)
                                                          acacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt                                    <  2:1322672/100‑1 (MQ=255)
                                                           cacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTc                                   >  1:2325859/1‑100 (MQ=255)
                                                            actggATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTc                                   <  2:951688/96‑1 (MQ=255)
                                                            acGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa                                  <  2:845820/100‑1 (MQ=255)
                                                            acGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa                                  >  2:579850/1‑100 (MQ=255)
                                                            acGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa                                  <  1:1512360/100‑1 (MQ=255)
                                                              gggATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATa                                <  1:1937414/100‑1 (MQ=255)
                                                               ggATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATaa                               >  2:1289121/1‑100 (MQ=255)
                                                               ggATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATaa                               >  1:2117114/1‑100 (MQ=255)
                                                               ggATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATa                                >  2:2077432/1‑99 (MQ=255)
                                                                  ttattGCGGTTGACTTTGAGAAATATAAAATACGAGTTTATAGCACCGCGAGCTCCAAATAATTTTATGACCACAATACGACAGGTTGTGGTCATaaaat                            <  2:2267106/100‑2 (MQ=18)
                                                                   tattgggatttACTTTGAGGAAGATAAAATGCGTGTTTATAGAACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATaaaaaa                           >  2:1812625/9‑100 (MQ=255)
                                                                    attGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAAt                          >  2:1744057/1‑100 (MQ=255)
                                                                    attGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAAt                          <  1:989179/100‑1 (MQ=255)
                                                                     ttGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAAtt                         <  1:1276257/100‑1 (MQ=255)
                                                                     ttGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAAtt                         >  2:1064473/1‑100 (MQ=255)
                                                                      tGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTg                        >  2:2560930/1‑100 (MQ=255)
                                                                       gCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAAtt                         >  2:577883/1‑98 (MQ=255)
                                                                       gCGGTTTACTTTGAGGAATATAAAATACGAGTTCATAGCACCTCGAGCGCCAAATAATTTTATGACCACAATACGACAGGTTGTGGTCATaaaattattt                       <  1:110150/100‑7 (MQ=17)
                                                                        cGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTc                      >  1:342113/1‑100 (MQ=255)
                                                                          gcttACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCTCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCg                     <  2:875388/97‑1 (MQ=255)
                                                                           tttACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGc                   <  1:1052096/100‑1 (MQ=255)
                                                                           tttACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGc                   <  1:1628621/100‑1 (MQ=255)
                                                                             tACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTg                 >  2:1919926/1‑100 (MQ=255)
                                                                              aCTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgg                <  1:2461436/100‑1 (MQ=255)
                                                                              aCTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgg                <  1:107299/100‑1 (MQ=255)
                                                                               cTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgga               >  1:69626/1‑100 (MQ=255)
                                                                               cTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgga               <  1:610380/100‑1 (MQ=255)
                                                                                tttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggag              <  1:1055587/100‑1 (MQ=255)
                                                                                tttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggag              <  2:2473284/100‑1 (MQ=255)
                                                                                tttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggag              >  1:1934215/1‑100 (MQ=255)
                                                                                 ttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagg             <  2:2524455/100‑1 (MQ=255)
                                                                                 ttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagg             <  1:898089/100‑1 (MQ=255)
                                                                                  tGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggg            >  1:1997682/1‑100 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggg            <  2:1130308/100‑1 (MQ=255)
                                                                                   gAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggga           >  1:1553301/1‑100 (MQ=255)
                                                                                    aGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagggag          >  1:2254164/1‑100 (MQ=255)
                                                                                     ggAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggga           >  2:849660/1‑98 (MQ=255)
                                                                                     ggAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGt         >  2:1753351/1‑100 (MQ=255)
                                                                                     ggAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGt         >  2:467507/1‑100 (MQ=255)
                                                                                     ggAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGt         >  2:302381/1‑100 (MQ=255)
                                                                                     ggAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGt         >  2:153459/1‑100 (MQ=255)
                                                                                       aaTATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTa       <  1:1765571/100‑1 (MQ=255)
                                                                                         tataAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTaaa     >  1:953724/1‑100 (MQ=255)
                                                                                          ataAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAAt    <  1:2226940/100‑1 (MQ=255)
                                                                                            aaaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc  <  2:730704/100‑1 (MQ=255)
                                                                                            aaaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc  <  1:579108/100‑1 (MQ=255)
                                                                                             aaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCGCAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc  <  2:1629546/99‑1 (MQ=255)
                                                                                             aaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc  <  1:229941/99‑1 (MQ=255)
                                                                                                 |                                                                                              
TATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCC  >  NC_000913/2493548‑2493739

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