Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,493,563 T→A 5.3% intergenic (‑309/+204) frc ← / ← yfdX formyl‑CoA transferase/protein YfdX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,493,5630TA 5.3% 451.1 / 9.3 142intergenic (‑309/+204)frc/yfdXformyl‑CoA transferase/protein YfdX
Reads supporting (aligned to +/- strand):  ref base T (64/70);  new base A (4/4);  total (68/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 144 total aligned reads displayed.

CTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAA  >  NC_000913/2493418‑2493711
                                                                                                                                                 |                                                                                                                                                    
  gCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCg                                                                                                                                                >  2:1002254/1‑150 (MQ=255)
   cATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGt                                                                                                                                               >  2:449564/1‑150 (MQ=255)
     tGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTCTTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTAc                                                                                                                                             <  2:1433280/150‑1 (MQ=255)
     tGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTAc                                                                                                                                             <  1:429489/150‑1 (MQ=255)
      ggTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCTCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCGTaat                                                                                                                                            <  1:1527240/150‑3 (MQ=255)
      ggTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCACGGGTATGATTACGCTTTAAAACGCGTaat                                                                                                                                            <  1:1106654/150‑3 (MQ=255)
      ggTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCGTaat                                                                                                                                            <  2:779946/150‑3 (MQ=255)
       gTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTAAAAAGCGTaatc                                                                                                                                           <  1:1007079/150‑4 (MQ=255)
        tAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCAt                                                                                                                                          >  2:758040/1‑150 (MQ=255)
        tAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCAt                                                                                                                                          >  1:998436/1‑150 (MQ=255)
         aaGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATa                                                                                                                                         <  2:264674/150‑1 (MQ=255)
         aaGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATa                                                                                                                                         >  2:798562/1‑150 (MQ=255)
          agtattGAGTAATAAAAGAAAAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATaa                                                                                                                                        >  2:320445/1‑150 (MQ=255)
             agtGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaac                                                                                                                                     <  1:506152/150‑1 (MQ=255)
              gtGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaaca                                                                                                                                    <  2:224102/150‑1 (MQ=255)
              gtGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaaca                                                                                                                                    >  2:885455/1‑150 (MQ=255)
               tGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATaaaa                                                                                                                                      >  1:304365/1‑147 (MQ=255)
               tGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaacat                                                                                                                                   >  2:723846/1‑150 (MQ=255)
                gAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaacata                                                                                                                                  >  1:965684/1‑150 (MQ=255)
                gAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaacata                                                                                                                                  >  2:1092425/1‑150 (MQ=255)
                gAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaacata                                                                                                                                  <  1:24575/150‑1 (MQ=255)
                  ggACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaacataac                                                                                                                                <  1:471119/150‑1 (MQ=255)
                  ggACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAaacataac                                                                                                                                <  1:1459892/150‑1 (MQ=255)
                       aaaaGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACataat                                                                                                                           <  2:1600387/150‑1 (MQ=255)
                       aaaaGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACataat                                                                                                                           >  2:1391827/1‑150 (MQ=255)
                        aaaGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACataata                                                                                                                          <  2:1202965/150‑1 (MQ=255)
                        aaaGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACataata                                                                                                                          <  1:1467485/150‑1 (MQ=255)
                                atatGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGAttt                                                                                                                  >  2:247330/1‑150 (MQ=255)
                                   tGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTtgtg                                                                                                               <  2:1059540/150‑1 (MQ=255)
                                     aTAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGtgat                                                                                                             <  1:247330/150‑1 (MQ=255)
                                          tatTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGacac                                                                                                        >  2:1516593/1‑150 (MQ=255)
                                             tGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATTATTACGCTTTTTAAAGCGTACCATAAAACATAACACAATAAATGATTTGTGATGACACggg                                                                                                     <  2:487482/150‑1 (MQ=255)
                                             tGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACggg                                                                                                     <  1:1391827/150‑1 (MQ=255)
                                               aTTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTCCAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGat                                                                                                   >  2:1171386/1‑150 (MQ=255)
                                                    tGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTg                                                                                              <  1:660498/150‑1 (MQ=255)
                                                       tAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCgg                                                                                           <  2:565980/150‑1 (MQ=255)
                                                         aTCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGtt                                                                                         >  1:196792/1‑150 (MQ=255)
                                                            agagagagTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTAc                                                                                      >  2:1403277/1‑150 (MQ=255)
                                                             gagagagTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACt                                                                                     >  1:348000/1‑150 (MQ=255)
                                                              agagagTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACtt                                                                                    <  2:1490397/150‑1 (MQ=255)
                                                                  agTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAg                                                                                <  1:1504577/150‑1 (MQ=255)
                                                                   gTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAgg                                                                               >  1:115266/1‑150 (MQ=255)
                                                                   gTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAgg                                                                               <  1:798562/150‑1 (MQ=255)
                                                                     cgcgACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGaa                                                                             <  1:566450/150‑1 (MQ=255)
                                                                       cgACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAata                                                                           <  1:1002254/150‑1 (MQ=255)
                                                                          cATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATaa                                                                        <  2:16945/150‑1 (MQ=255)
                                                                          cATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATATACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATaa                                                                        >  2:1182126/1‑150 (MQ=255)
                                                                              tgtgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATg                                                                    >  2:200512/1‑150 (MQ=255)
                                                                               gtgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGc                                                                   >  2:1341165/1‑150 (MQ=255)
                                                                               gtgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGc                                                                   >  2:1202500/1‑150 (MQ=255)
                                                                               gtgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGAATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGc                                                                   >  1:233103/1‑150 (MQ=255)
                                                                                tgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCg                                                                  >  2:934361/1‑150 (MQ=255)
                                                                                tgATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCg                                                                  <  2:115266/150‑1 (MQ=255)
                                                                                 taTGTTATTCTTACAACTTTACCAAGCGCCGGGCATTTAAAACCAAGGGTATGATTACGCTTTTTAATGCGTACCATAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAg                                                                <  1:714472/149‑1 (MQ=255)
                                                                                  aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAg                                                                <  1:311149/150‑1 (MQ=255)
                                                                                  aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAg                                                                <  1:320445/150‑1 (MQ=255)
                                                                                  aTGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAg                                                                <  1:723846/150‑1 (MQ=255)
                                                                                   tGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGt                                                               >  1:380059/1‑150 (MQ=255)
                                                                                     ttATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGttt                                                             <  1:1408877/150‑1 (MQ=255)
                                                                                     cgcttCTTACAGTTGTAAGAATCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATaaaa                                                                      >  2:488352/4‑141 (MQ=255)
                                                                                     cgcttCTTACAGTTGTAAGAATCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATaaaa                                                                      <  1:488352/138‑1 (MQ=255)
                                                                                       aTTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTtat                                                           <  1:934361/150‑1 (MQ=255)
                                                                                        ttCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTtata                                                          <  2:443020/150‑1 (MQ=255)
                                                                                        ttCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTtata                                                          <  2:965684/150‑1 (MQ=255)
                                                                                               aaatttAAGAATCGGCTGTGATTTACAATCAATGGTATGCTTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACcgc                                                   <  1:296550/145‑1 (MQ=255)
                                                                                                aCTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACcgcg                                                  <  1:1268704/150‑1 (MQ=255)
                                                                                                 cTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGa                                                 >  2:336176/1‑150 (MQ=255)
                                                                                                  tGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAg                                                <  1:958488/150‑1 (MQ=255)
                                                                                                    tAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAgcg                                              >  1:935515/1‑150 (MQ=255)
                                                                                                     aagaagCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAgcgc                                             <  1:1192594/150‑1 (MQ=255)
                                                                                                      agaagCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGcc                                            >  1:858731/1‑150 (MQ=255)
                                                                                                       gaagCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTTTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCa                                           >  2:1280264/1‑150 (MQ=255)
                                                                                                       gaagCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCa                                           <  1:885455/150‑1 (MQ=255)
                                                                                                          gcgcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAaat                                        <  1:200512/150‑1 (MQ=255)
                                                                                                          gcgcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAaat                                        >  1:586185/1‑150 (MQ=255)
                                                                                                           cgcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAaata                                       <  2:714445/150‑1 (MQ=255)
                                                                                                            gcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATATTAAAGGAGTTGTGAGGCAACGGGATTATTGCGGTGTACTTTGTGGAATATTAATTGCGAGTTTATAGCACCGCGAGCGCCAaatat                                      >  2:1092591/1‑149 (MQ=255)
                                                                                                            gcCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAaataa                                      >  2:1076700/1‑150 (MQ=255)
                                                                                                                  gATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGGCACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTAt                                >  1:709474/1‑150 (MQ=255)
                                                                                                                    tttACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGa                              >  1:1006613/1‑150 (MQ=255)
                                                                                                                      tACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGAcc                            >  1:164780/1‑150 (MQ=255)
                                                                                                                            caatcgTATGATTACGCTTTTAAAAGCGTACCAAAAAACATAACATAATAAATGATTTGTGATGACACCGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAAcc                      >  2:808448/6‑150 (MQ=255)
                                                                                                                               gggTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGt                   <  1:758040/150‑1 (MQ=255)
                                                                                                                                ggTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTg                    <  2:397767/148‑1 (MQ=255)
                                                                                                                                ggTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTc                  >  1:1545613/1‑150 (MQ=255)
                                                                                                                                 gTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCg                 <  2:61111/150‑1 (MQ=255)
                                                                                                                                  ttatggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGt                >  1:805472/7‑150 (MQ=255)
                                                                                                                                    tGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTAt              <  2:815257/150‑1 (MQ=255)
                                                                                                                                      ggtACGCTTTAAAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGAGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtg            >  2:1431458/3‑150 (MQ=255)
                                                                                                                                       ttACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgt           <  2:887423/150‑1 (MQ=255)
                                                                                                                                        tACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATtgtg          >  2:1608310/1‑150 (MQ=255)
                                                                                                                                         aCGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATg                                                                    >  1:1254487/1‑91 (MQ=255)
                                                                                                                                         aCGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATg                                                                    <  2:1254487/91‑1 (MQ=255)
                                                                                                                                          cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt        <  2:1006613/150‑1 (MQ=255)
                                                                                                                                          cGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGt        <  1:1092425/150‑1 (MQ=255)
                                                                                                                                           gCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGta       >  2:1593809/1‑149 (MQ=255)
                                                                                                                                           gCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTc       <  2:230264/150‑1 (MQ=255)
                                                                                                                                            cTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCa      <  1:1182126/150‑1 (MQ=255)
                                                                                                                                              ttttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATa    >  1:435068/1‑150 (MQ=255)
                                                                                                                                                ttAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATaaa  >  1:162033/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGCATGGTAAGTAGTGAGGACGAAAAGAAGAATATGATAAATATTGATTCGTGGTAATCAGAGAGAGTCGCGACATTTGTGATGTTATTCTTACAACTGTAAGAAGCGCCGGTGATTTACAAGCAAGGGTATGATTACGCTTTTTAAAGCGTACCATAAAACATAACATAATAAATGATTTGTGATGACACGGGATTATTGCGGTTTACTTTGAGGAATATAAAATGCGAGTTTATAGCACCGCGAGCGCCAAATAATTTTATGACCACAACCTGTCGTATTGTGGTCATAAA  >  NC_000913/2493418‑2493711

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