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,6450GA96.9% 397.4 / 1.8 129intergenic (‑391/+122)frc/yfdXformyl‑CoA transferase, NAD(P)‑binding/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (64/60);  total (66/63)
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.79e-01

Only 100 of 130 total aligned reads displayed.

TATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCTTTAC  >  NC_000913/2493548‑2493744
                                                                                                 |                                                                                                   
tATGATTACGCTTTTTAAAGCGTACCATAAAACATATCATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACg                                                                                                   <  1:895327/100‑1 (MQ=255)
tATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACg                                                                                                   <  1:2716287/100‑1 (MQ=255)
  tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGa                                                                                                  >  2:1806571/1‑99 (MQ=255)
  atggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAg                                                                                                 >  2:1185816/5‑100 (MQ=255)
      tACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTta                                                                                             <  2:2082364/100‑1 (MQ=255)
       aCGCTTTTTAAATCGTACCATAAAAAATAATATAATAAATAATTTGTTTTGACTCTGGATTATTGCGGTTTACTTTTACGAATATAAAATACGATTatat                                                                                            >  2:960561/1‑100 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATACCATAATAATTGATTTGTGGTGACACGGGAATATTGCGGTTTAATTTGAGGAATATAAAATACAAGTTtata                                                                                           >  1:1203345/1‑100 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                           <  1:19881/100‑1 (MQ=255)
        cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTtata                                                                                           >  1:1923859/1‑100 (MQ=255)
         gCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAg                                                                                          <  1:1014444/100‑1 (MQ=255)
          cTTTTTAAAGCGTATCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAAGATAAAATACGAGTTtata                                                                                           <  1:1231858/98‑1 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                         >  1:1347403/1‑100 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                         <  2:854224/100‑1 (MQ=255)
          cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACAAGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGc                                                                                         >  1:941919/1‑100 (MQ=255)
           tttttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCa                                                                                        <  2:2681697/100‑1 (MQ=255)
           ccgatctaGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCa                                                                                        <  1:1250574/93‑1 (MQ=255)
             tttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCAcc                                                                                      >  1:2492668/1‑100 (MQ=255)
             taaaaaGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCAcc                                                                                      >  2:395211/4‑100 (MQ=255)
              ttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACcg                                                                                     >  2:520102/1‑100 (MQ=255)
               tAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACcgc                                                                                    >  1:451253/1‑100 (MQ=255)
                aaaGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACcgcg                                                                                   >  1:2431388/1‑100 (MQ=255)
                      tACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGcc                                                                             >  1:683107/1‑100 (MQ=255)
                          aTAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaat                                                                         >  2:2568809/1‑100 (MQ=255)
                            aaaaCATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaataa                                                                       <  1:2301878/100‑1 (MQ=255)
                             aaaCATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAaataat                                                                      >  2:1553201/1‑100 (MQ=255)
                              aacataacatAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAAtt                                                                     >  1:2380267/1‑100 (MQ=255)
                                cataacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAAtttt                                                                   <  1:1322038/100‑1 (MQ=255)
                                  taacataaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTAt                                                                 >  2:1308421/1‑100 (MQ=255)
                                   aacattaTAAATGATTTGTGATGACATGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                                <  1:373836/100‑1 (MQ=255)
                                   aacacaaTAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATg                                                                <  1:2105873/100‑1 (MQ=255)
                                       taataaATGATTTGTGATGATACGGGATTATTGCGGTTTACTTTAAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACca                                                            <  1:520105/100‑1 (MQ=255)
                                       taataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACca                                                            >  2:2261595/1‑100 (MQ=255)
                                       taataaATGATTTGTGATGACACGGGAATATTGCCGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACca                                                            >  2:39802/1‑100 (MQ=255)
                                         ataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACcaca                                                          <  1:1349655/100‑1 (MQ=255)
                                         ataaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACcaca                                                          <  2:696072/100‑1 (MQ=255)
                                          taaATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACaa                                                         >  1:2019973/1‑100 (MQ=255)
                                            aaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                       >  1:188434/1‑100 (MQ=255)
                                            aaTGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                       <  2:1778590/100‑1 (MQ=255)
                                             aTGATTTGTGATTACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                       <  1:504434/99‑1 (MQ=255)
                                                aTTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                                                       >  2:1232033/1‑96 (MQ=255)
                                                  ttGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGt                                                 >  1:1288036/1‑100 (MQ=255)
                                                   tgtgATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTa                                                >  2:1089338/1‑100 (MQ=255)
                                                    gtgATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTAt                                               >  2:1404886/1‑100 (MQ=255)
                                                     tgatgaCACGGGAGTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTAtt                                              >  1:2063143/1‑100 (MQ=255)
                                                       atgaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATAACCACAACCTGTCGTAtt                                              >  1:1633495/1‑98 (MQ=255)
                                                       atgaCACGGGACTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgt                                            <  2:1210602/100‑1 (MQ=255)
                                                        tgaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgtg                                           <  1:272273/100‑1 (MQ=255)
                                                         gaCACGGGATTATTGCGGTTTACTTTGAGTAATATAATATACGAGTTTATAGCACCGCGAGCGCCAATTAATTTTATGACCACAACCTGTCGTATtgtg                                           <  2:183151/99‑1 (MQ=255)
                                                         gaCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTgg                                          <  1:166023/100‑1 (MQ=255)
                                                          acacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt                                         <  2:1030288/100‑1 (MQ=255)
                                                          acacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt                                         <  2:1299443/100‑1 (MQ=255)
                                                          acacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt                                         >  1:1030295/1‑100 (MQ=255)
                                                          acacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt                                         >  1:2395684/1‑100 (MQ=255)
                                                           cacGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgt                                            >  2:1666062/1‑96 (MQ=255)
                                                           cacGGGATTATTGCGGTTTACTTAGAGGAATATAAAATACGAGTTTATAACACCGCGAGAGAAAAATAATTTTATGACCACAACCTGTCGTATTGTGGTc                                        >  2:708538/1‑100 (MQ=255)
                                                            acGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa                                       <  1:111814/100‑1 (MQ=255)
                                                             cGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCAt                                      >  2:2464007/1‑100 (MQ=255)
                                                             cGGGATTATTGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCAt                                      >  1:2478613/1‑100 (MQ=255)
                                                                    attGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAATACGACAGGTTGTGGTCATaaaat                                 <  2:2433965/98‑2 (MQ=21)
                                                                      tGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTg                             >  2:2425839/1‑100 (MQ=255)
                                                                      tGCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTg                             <  2:532524/100‑1 (MQ=255)
                                                                       gCGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGt                            >  2:416847/1‑100 (MQ=255)
                                                                        cGGTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTc                           <  1:2777014/100‑1 (MQ=255)
                                                                          gTTTACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCg                          >  1:2259375/1‑99 (MQ=255)
                                                                           tttACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGc                        <  1:1013067/100‑1 (MQ=255)
                                                                            ttACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGca                       >  2:1013061/1‑99 (MQ=255)
                                                                             tACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTg                      >  1:246955/1‑100 (MQ=255)
                                                                             tACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTg                      <  1:80154/100‑1 (MQ=255)
                                                                             tACTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCATATTGTGAACATAAAAAATTGTCAGc                        >  1:2965667/1‑98 (MQ=255)
                                                                              aCTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgg                     >  2:1213980/1‑100 (MQ=255)
                                                                              aCTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgg                     <  1:1228560/100‑1 (MQ=255)
                                                                              aCTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgg                     <  2:684649/100‑1 (MQ=255)
                                                                               cTTTTAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCTAAATAATTTTATGAACACAACCTGTCGTAATGTGGTCATAAAAAATTGTCGGCCgga                    >  2:57278/1‑100 (MQ=255)
                                                                               cTTTGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTgga                    >  1:2914707/1‑100 (MQ=255)
                                                                                tttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggag                   <  2:592461/100‑1 (MQ=255)
                                                                                tttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggag                   <  2:933514/100‑1 (MQ=255)
                                                                                 ttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagg                  >  2:2249616/1‑100 (MQ=255)
                                                                                 ttGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagg                  >  2:2515618/1‑100 (MQ=255)
                                                                                  taagtaaTATAAATTACCATTTTCTTGCACCTCTATTGTCAAATAATTATATGACTACAACCTTTCTTATTGTGTTCATAAAAAATTGTAGGCTggag                   <  2:2681335/93‑1 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAACAAATTGTCGGCTggacgg                 <  2:2852113/100‑1 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggg                 <  1:805531/100‑1 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggg                 >  1:2748540/1‑100 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaggg                 <  2:2456407/100‑1 (MQ=255)
                                                                                  tGAGGAATATAAAATACGAGATTATAGCACCGCGAGCGCCAAATAATTTTATGACAACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggaag                  >  2:1117234/1‑97 (MQ=255)
                                                                                    aGGAATATAAAATACGAGTTTATAGCACCGCGAGTGCCTAATAATTTTATGACCACAAGCTGTTGTATTGTGGTCATAAAAAATTGTCGGCTggagggag               <  1:2507664/100‑1 (MQ=255)
                                                                                    aGGAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTggagggag               >  1:1597058/1‑100 (MQ=255)
                                                                                     ggAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGt              <  1:416850/100‑1 (MQ=255)
                                                                                      gAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGtt             <  2:1890023/100‑1 (MQ=255)
                                                                                      gAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGtt             <  2:1063939/100‑1 (MQ=255)
                                                                                      gAATATAAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGt              <  1:2261620/99‑1 (MQ=255)
                                                                                         tataAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTaaa          <  2:2439788/100‑1 (MQ=255)
                                                                                         tataAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTAGAGGGAGTTaaa          >  1:2811541/1‑100 (MQ=255)
                                                                                          ataAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGGTAAAt         >  2:1633032/1‑100 (MQ=255)
                                                                                           taAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCAGCTGGAGGGAGTTAAATc        >  1:604179/1‑100 (MQ=255)
                                                                                           taAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTATAGGCTGGAAGGAGTTAAATc        >  1:1106550/1‑100 (MQ=255)
                                                                                           taAAATACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCAGTCGTATTGTGGTCATAAAAAATAGTCGGATGGAGGGAGTTaaa          >  1:1907817/1‑98 (MQ=255)
                                                                                            aaaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATcc       >  1:1049033/1‑100 (MQ=255)
                                                                                             aaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCt      <  1:2978726/100‑1 (MQ=255)
                                                                                             aaaTACGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCt      <  1:2978826/100‑1 (MQ=255)
                                                                                                 gCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCTTTAc  <  1:2629033/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAAAAATTGTCGGCTGGAGGGAGTTAAATCCTTTAC  >  NC_000913/2493548‑2493744

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: