Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,000,237 T→A 9.4% intergenic (+92/+55) uvrD → / ← yigE DNA‑dependent ATPase I and helicase II/DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,000,2370TA 9.4% 269.2 / 14.2 95intergenic (+92/+55)uvrD/yigEDNA‑dependent ATPase I and helicase II/DUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base T (42/44);  new base A (6/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCA  >  NC_000913/4000089‑4000385
                                                                                                                                                    |                                                                                                                                                     
catggCCTGGGTATTAAATGGCTGGTTGCGGCATTCGCCCGGCTGGAGTCTGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGttt                                                                                                                                                      <  2:271041/147‑1 (MQ=255)
cAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCTGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGttt                                                                                                                                                      <  1:1133009/150‑1 (MQ=255)
cAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCGCCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGttt                                                                                                                                                      >  2:1784726/1‑150 (MQ=255)
 aGGGCCAGGGTATTTAATGTCTGGTGGCGGCATACGCCCTGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTa                                                                                                                                                     <  1:835125/150‑1 (MQ=255)
  gggCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAc                                                                                                                                                    >  2:1741157/1‑150 (MQ=255)
   ggCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCTGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcg                                                                                                                                                   <  1:1625379/150‑1 (MQ=255)
      cAGGGTATTAACTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccg                                                                                                                                                <  2:713672/150‑1 (MQ=255)
      cAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccg                                                                                                                                                <  1:879520/150‑1 (MQ=255)
       aGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccgc                                                                                                                                               <  1:911364/150‑1 (MQ=255)
         ggTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCAt                                                                                                                                             <  1:1325272/150‑1 (MQ=255)
         ggTATTAAATGGCTGGTGGCGGAATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCAt                                                                                                                                             <  1:1797317/150‑1 (MQ=255)
           tATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATcc                                                                                                                                           <  1:1620933/150‑1 (MQ=255)
           tATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATcc                                                                                                                                           <  1:1810609/150‑1 (MQ=255)
           tATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATcc                                                                                                                                           <  1:390126/150‑1 (MQ=255)
            aTTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCCGATAGGGCGTTTACGCCGCATCCg                                                                                                                                          >  1:54729/1‑150 (MQ=255)
             ttAAATGGCTGGGGCCGGCATACGCCCGGCTGGAGTCGGGGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCgg                                                                                                                                         <  1:185504/150‑1 (MQ=255)
               aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCCCCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCa                                                                                                                                       >  2:668367/1‑150 (MQ=255)
               aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCa                                                                                                                                       <  1:72250/150‑1 (MQ=255)
               aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGGCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCa                                                                                                                                       <  2:1798426/150‑1 (MQ=255)
                  tGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCGGcctac                                                                                                                                    <  2:150374/150‑5 (MQ=255)
                    gCTGGTGGCGGCATACGCCCGGCTCGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTCCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGAATCTgc                                                                                                                                  >  2:777698/1‑150 (MQ=255)
                     cTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATGCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcg                                                                                                                                 >  2:942810/1‑150 (MQ=255)
                        gtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCa                                                                                                                              <  1:495181/150‑1 (MQ=255)
                             ggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTc                                                                                                                         <  2:1047794/150‑1 (MQ=255)
                             ggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTc                                                                                                                         <  1:1547613/150‑1 (MQ=255)
                               cATACGCCCGTCTGGAGTCGTTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAg                                                                                                                       <  1:1433559/150‑1 (MQ=255)
                                aTACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCAGTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGt                                                                                                                      >  2:218905/1‑150 (MQ=255)
                                   cGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATc                                                                                                                   >  2:1088381/1‑150 (MQ=255)
                                     cccGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTg                                                                                                                 <  1:1660240/150‑1 (MQ=255)
                                     cccGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTg                                                                                                                 <  2:1495710/150‑1 (MQ=255)
                                          cTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGAcaca                                                                                                            >  1:36320/1‑150 (MQ=255)
                                                 cGGTGTAACCATTCAGGATGCGGTGCACAGCAACTTATTTGGCCTAAAAAATCATTCAGATTCAATAAAATGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGAATATGACACAAAACTag                                                                                                     >  2:96143/1‑149 (MQ=38)
                                                  ggTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATATGACACAAAACTATc                                                                                                    >  1:1343357/1‑150 (MQ=255)
                                                   gtgtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCg                                                                                                   <  2:1235000/150‑1 (MQ=255)
                                                     gtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtt                                                                                                 <  2:870772/150‑1 (MQ=255)
                                                     gtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtt                                                                                                 <  2:862924/150‑1 (MQ=255)
                                                       aaCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtttt                                                                                               <  2:1092282/150‑1 (MQ=255)
                                                        aCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTa                                                                                              <  2:1234796/150‑1 (MQ=255)
                                                        aCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTa                                                                                              <  2:81586/150‑1 (MQ=255)
                                                         cGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTaa                                                                                             >  2:231734/1‑150 (MQ=255)
                                                          gTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAc                                                                                            >  2:1757953/1‑150 (MQ=255)
                                                          gTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAc                                                                                            <  2:333879/150‑1 (MQ=255)
                                                          gTTGCCAGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAc                                                                                            <  2:1209139/150‑1 (MQ=255)
                                                            tGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCt                                                                                          >  1:46898/1‑150 (MQ=255)
                                                                ggATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGTCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTcc                                                                                      <  2:1399786/150‑1 (MQ=255)
                                                                ggATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTcc                                                                                      <  2:1490865/150‑1 (MQ=255)
                                                                     cGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                 <  2:926136/150‑1 (MQ=255)
                                                                     cGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                 <  2:1086619/150‑1 (MQ=255)
                                                                      ggTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGa                                                                                >  2:1579951/1‑150 (MQ=255)
                                                                               cACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCa                                                                      >  1:639273/1‑150 (MQ=255)
                                                                                aCCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGTAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCG‑TCGAGCGGaaaggtt                                                                     <  1:1380672/150‑5 (MQ=39)
                                                                                 ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                    >  2:1485875/1‑150 (MQ=255)
                                                                                 ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATg                                                                    >  2:1695851/1‑150 (MQ=255)
                                                                                 ccTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGAGTTTACGCCGCATCCGCCATCTGCGGCATCTTCAGTATCTGGCACAAAACTATCGATTTAACCTTTCCGATCGA‑CGGAAATCATg                                                                    >  2:1775609/1‑150 (MQ=255)
                                                                                          gCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACaaaa                                                           >  1:1830221/1‑150 (MQ=255)
                                                                                           ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAg                                                          >  1:452164/1‑150 (MQ=255)
                                                                                           ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAg                                                          >  2:1068013/1‑150 (MQ=255)
                                                                                           ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAg                                                          >  1:1536249/1‑150 (MQ=255)
                                                                                           ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAg                                                          >  1:1367937/1‑150 (MQ=255)
                                                                                             tAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATAATGGTGACAAAAGGa                                                        >  1:1813994/1‑150 (MQ=255)
                                                                                             tAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAAATATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAggt                                                        >  1:180992/1‑149 (MQ=255)
                                                                                              aaaaaaTCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGAt                                                       <  2:143442/150‑1 (MQ=255)
                                                                                                       ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGcc                                              >  1:1637666/1‑150 (MQ=255)
                                                                                                        tCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTTACAAAAGGATAGCGTTGccc                                             >  2:73074/1‑149 (MQ=255)
                                                                                                            aTTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACggg                                         <  1:616160/150‑1 (MQ=255)
                                                                                                                aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAACACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCCCGGGATTg                                     >  2:847535/1‑150 (MQ=255)
                                                                                                                aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                     >  1:518435/1‑150 (MQ=255)
                                                                                                                aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                     >  2:1098304/1‑150 (MQ=255)
                                                                                                                 aTAAATTGCAACGTCATGTAGGACGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc                                    <  1:208434/150‑1 (MQ=255)
                                                                                                                   aaaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGccc                                  >  2:420289/1‑150 (MQ=255)
                                                                                                                   aaaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGacc                                  >  1:551854/1‑150 (MQ=255)
                                                                                                                    aaTTTCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                 <  2:778136/150‑1 (MQ=255)
                                                                                                                    aaTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGcccc                                 <  2:778133/150‑1 (MQ=255)
                                                                                                                         cAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCtt                            >  2:867472/1‑150 (MQ=255)
                                                                                                                           aCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                          <  1:1251158/150‑1 (MQ=255)
                                                                                                                           aCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                          <  2:1536249/150‑1 (MQ=255)
                                                                                                                            cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                         <  2:1418048/150‑1 (MQ=255)
                                                                                                                            cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCAGCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                         <  1:305744/150‑1 (MQ=255)
                                                                                                                             gTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCat                        >  2:848769/1‑150 (MQ=255)
                                                                                                                              tCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCata                       <  2:645109/150‑1 (MQ=255)
                                                                                                                               cATGTAGGCCGGATAGGGCGTTTACGCCGCATCAGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAGTCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCatat                      >  1:548624/1‑150 (MQ=255)
                                                                                                                                gcagatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCAATCTTCAGTATCTGACACAAAACTATCGTTTTATCCTTTCCGCTCGA‑CGGAATTCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCAtata                     >  1:1098981/7‑150 (MQ=38)
                                                                                                                                 cagatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATAc                    >  1:116344/6‑150 (MQ=255)
                                                                                                                                   tAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACat                  >  1:455874/1‑150 (MQ=255)
                                                                                                                                   gatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACat                  >  1:1628669/4‑150 (MQ=255)
                                                                                                                                    cggCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCATCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACata                 <  2:1577908/149‑1 (MQ=255)
                                                                                                                                    aGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACata                 >  2:669876/1‑150 (MQ=255)
                                                                                                                                     tgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGCCAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCGTTCATATACatat                >  2:1737029/2‑150 (MQ=255)
                                                                                                                                     tgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTATCATATACattt                >  2:1324546/2‑148 (MQ=255)
                                                                                                                                     tgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACGATCGTGTTACCCTTTCCGCTCGA‑CGGAAATCATGGGGACAACAGGATAGCGTTGCCACGGGATTGCCCCGCCTTGCATATACatat                >  2:1318213/2‑150 (MQ=255)
                                                                                                                                      gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtg               >  1:1488200/1‑150 (MQ=255)
                                                                                                                                      gCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtg               <  1:1741157/150‑1 (MQ=255)
                                                                                                                                       ccGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGC‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgt              >  1:1588615/1‑150 (MQ=255)
                                                                                                                                       ccGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgt              >  1:1749430/1‑150 (MQ=255)
                                                                                                                                            tAGGGCGGTTACGCCGCATGCGGCATCTGCGCCATCTTCAGTATGTGACACAAAACTATCGTTTTAACCTTTGCGCTCGA‑CGGAAATCATTGTGACAAAAGGATAGGGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGaatt         <  2:1478030/150‑3 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGt      >  1:1606757/1‑150 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGt      >  1:387532/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCAGCTGCGCCATCTTCAGTATCTGACACAATACTATCGTTTTAACCTTTCCGATCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGc    <  1:319502/150‑1 (MQ=255)
                                                                                                                                                   tttACGCCGCATCCGGCATCTGCGACATCTTCAGTATATGACACAAAACTATCGTTTTAACATTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGTATTGCCCCGCATTTCAAATACATATGTGAAATCGTGCCa  <  1:1225649/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGA‑CGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCA  >  NC_000913/4000089‑4000385

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