Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,991,579 T→A W143R (TGG→AGG)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,5790TA94.7% 372.9 / 8.2 113W143R (TGG→AGG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base A (53/54);  total (56/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

Only 100 of 113 total aligned reads displayed.

ATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGC  >  NC_000913/3991435‑3991726
                                                                                                                                                |                                                                                                                                                   
aTTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGgc                                                                                                                                                 <  6:59905/149‑1 (MQ=255)
   aCTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGgccgc                                                                                                                                              >  1:308478/1‑149 (MQ=255)
     tGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCt                                                                                                                                            >  6:104720/1‑149 (MQ=255)
       gtgtATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCg                                                                                                                                          >  4:156402/1‑149 (MQ=255)
       gtgtATACACCATGGGCAGCAACTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTTGATCGATAGAGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACATGGCCGCCTCg                                                                                                                                          >  5:33739/1‑149 (MQ=255)
         gtATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCt                                                                                                                                        >  8:507797/1‑149 (MQ=255)
          tataCACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCt                                                                                                                                        <  3:156402/148‑1 (MQ=255)
            taCACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTggg                                                                                                                                     >  2:395952/1‑149 (MQ=255)
            taCACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTAGTTCCTCTGACCGTGATATCTGGGTCTGTCATCAAGCCTAGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTgg                                                                                                                                      <  3:17271/148‑1 (MQ=255)
              cacCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGgtg                                                                                                                                   <  8:508304/149‑1 (MQ=255)
              cacCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGgt                                                                                                                                    <  2:453572/148‑1 (MQ=255)
               acCATTGGCAGCACCTCTTCCGTAGTGCAAAGTTGTTCCTCTTACCTGGATATCTTGGTCTGTCATCAATCCTGGCTCGATAGCTAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGgtgt                                                                                                                                  <  6:33739/149‑1 (MQ=255)
               acCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGgtgt                                                                                                                                  >  4:457818/1‑149 (MQ=255)
                ccATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGtgtg                                                                                                                                 <  7:507797/149‑1 (MQ=255)
                   tGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                              <  2:308478/149‑1 (MQ=255)
                    gggCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGTGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                              <  5:482516/148‑1 (MQ=255)
                    gggCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTgg                                                                                                                                >  6:482516/1‑146 (MQ=255)
                          caCCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATAAATCCTGGCTCCATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGct                                                                                                                       >  7:260597/1‑149 (MQ=255)
                             cTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGgt                                                                                                                                    >  8:190769/1‑133 (MQ=255)
                             cTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGgt                                                                                                                                    <  7:190769/133‑1 (MQ=255)
                             cTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                              <  2:18151/139‑1 (MQ=255)
                             cTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                              >  1:18151/1‑139 (MQ=255)
                             cTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGcttct                                                                                                                    >  8:43402/1‑149 (MQ=255)
                                gTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                              >  5:58402/1‑136 (MQ=255)
                                gTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGaa                                                                                                                              <  6:58402/136‑1 (MQ=255)
                                gTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCttcttc                                                                                                                  <  3:161074/148‑1 (MQ=255)
                                gTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTcc                                                                                                                 <  3:167773/149‑1 (MQ=255)
                                gTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTcc                                                                                                                 >  8:510709/1‑149 (MQ=255)
                                   cGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtga                                                                                                              <  7:43402/149‑1 (MQ=255)
                                   cGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtga                                                                                                              >  2:144162/1‑149 (MQ=255)
                                    gTAGGGCAAAGTTGTTCCTCTGAGCTGGATATCTGTGTCTGTCATCTATACTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAATCAGGGCCGCCTCGCTGTGTGTGGAAGTCAGCTTCTTCCtgat                                                                                                             <  3:455168/149‑1 (MQ=255)
                                        ggCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtgattgat                                                                                                         >  2:71194/1‑149 (MQ=255)
                                        ggCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtgattgat                                                                                                         <  1:144162/149‑1 (MQ=255)
                                           aaaGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTgg                                                                                                                                      <  6:226757/117‑1 (MQ=255)
                                           aaaGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTgg                                                                                                                                      >  5:226757/1‑117 (MQ=255)
                                                     ctctGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGAGGAAAACCGCTTcc                                                                                            >  4:20817/1‑149 (MQ=255)
                                                            cTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATaa                                                                                     >  2:184316/1‑149 (MQ=255)
                                                               gATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAggg                                                                                                                                                  >  7:268242/1‑85 (MQ=255)
                                                               gATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAggg                                                                                                                                                  <  8:268242/85‑1 (MQ=255)
                                                               gATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAAt                                                                                    <  4:392063/147‑1 (MQ=255)
                                                               gATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATg                                                                                   >  1:69304/1‑148 (MQ=255)
                                                               gATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATg                                                                                   >  1:219656/1‑148 (MQ=255)
                                                                 tatCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTcc                                                                                                                 >  5:227714/1‑116 (MQ=255)
                                                                 tatCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTcc                                                                                                                 <  6:227714/116‑1 (MQ=255)
                                                                  atCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAg                                                                               >  6:430222/1‑149 (MQ=255)
                                                                     tGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGgcc                                                                                                                                                >  1:226899/1‑81 (MQ=255)
                                                                     tGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGgcc                                                                                                                                                <  2:226899/81‑1 (MQ=255)
                                                                     tGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCgg                                                                            >  8:447021/1‑149 (MQ=255)
                                                                       ggTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTg                                                                                                                                       >  4:467660/1‑88 (MQ=255)
                                                                       ggTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTg                                                                                                                                       <  3:467660/88‑1 (MQ=255)
                                                                           tGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtg                                                                                                               <  7:406139/108‑1 (MQ=255)
                                                                           tGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtg                                                                                                               >  8:406139/1‑108 (MQ=255)
                                                                            gTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCtgat                                                                                                             <  3:195308/109‑1 (MQ=255)
                                                                            gTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTc                                                                                         >  5:265601/1‑129 (MQ=255)
                                                                            gTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTc                                                                                         <  6:265601/129‑1 (MQ=255)
                                                                             tcatcaATCCTGGCTCGATAGCGAAGAGCGCCAATTTCTACAACGTAAATGTAGCCTTCTGGAAAACAGGGCCGCCTCGCTGGGGGTGGAAGTCAGc                                                                                                                        <  8:433628/97‑1 (MQ=255)
                                                                             tcatcaATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTgg                                                                    <  1:395952/149‑1 (MQ=255)
                                                                               atcaATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTgggg                                                                  >  8:491428/1‑149 (MQ=255)
                                                                               atcaATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTgggg                                                                  >  1:484515/1‑149 (MQ=255)
                                                                               atcaATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTggg                                                                   <  2:69304/148‑1 (MQ=255)
                                                                                 caATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTgggggg                                                                >  2:242015/1‑149 (MQ=255)
                                                                                   aTCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGc                                                               >  1:31983/1‑148 (MQ=255)
                                                                                      cTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCa                                                                          <  2:471257/134‑1 (MQ=255)
                                                                                      cTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCa                                                                          >  1:471257/1‑134 (MQ=255)
                                                                                      cTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGaa                                                            <  7:510709/148‑1 (MQ=255)
                                                                                      cTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAg                                                           <  1:71194/149‑1 (MQ=255)
                                                                                         gCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGa                                                          <  8:260597/147‑1 (MQ=255)
                                                                                         gCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGAtt                                                        <  2:219656/149‑1 (MQ=255)
                                                                                           tCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATtgt                                                      <  2:484515/149‑1 (MQ=255)
                                                                                             gATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTgg                                                    <  4:114684/149‑1 (MQ=255)
                                                                                             gATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTgg                                                    <  3:20817/149‑1 (MQ=255)
                                                                                             gATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGGGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATtgt                                                      <  7:521239/147‑1 (MQ=255)
                                                                                              aTAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGGGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGc                                                   <  1:184316/149‑1 (MQ=255)
                                                                                                  cGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCa                                               >  6:556647/1‑149 (MQ=255)
                                                                                                            cAATTGCTACAACGTAAATGTAGCCTGCTTGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTggg                                                                   <  1:162148/119‑1 (MQ=255)
                                                                                                            cAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTggg                                                                   >  2:162148/1‑119 (MQ=255)
                                                                                                                   tACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATActgctg                              >  6:490367/1‑149 (MQ=255)
                                                                                                                   tACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGAtt                                                        <  8:273208/123‑1 (MQ=255)
                                                                                                                   tACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGAtt                                                        >  7:273208/1‑123 (MQ=255)
                                                                                                                   tACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCAc                                              <  3:501035/133‑1 (MQ=255)
                                                                                                                   tACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCAc                                              >  4:501035/1‑133 (MQ=255)
                                                                                                                   tACAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATActgctg                              >  2:301526/1‑149 (MQ=255)
                                                                                                                    aCAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTc                                                 <  4:175547/129‑1 (MQ=255)
                                                                                                                    aCAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACtgctgc                             >  7:13622/1‑149 (MQ=255)
                                                                                                                     cAACGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTgctgct                            <  2:31983/149‑1 (MQ=255)
                                                                                                                      aaCGTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGAGTGGAAGTCAGCTTCTTCCTGATTGATGAACACCGCTTCCCTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCtt                           >  5:248999/1‑149 (MQ=255)
                                                                                                                         gTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTTATTGATGAAAACCGCTTCCGTCATAATTAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCAcc                                             <  4:118175/128‑1 (MQ=255)
                                                                                                                         gTAAATGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGAc                        >  5:489331/1‑149 (MQ=255)
                                                                                                                         gTAAATGTAGCCTGCTGGAAAACAGGGCAGCCTAGCTGGGTTTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCAcc                                             >  3:118175/1‑128 (MQ=255)
                                                                                                                           aaaTGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGa                                                             >  3:545638/1‑110 (MQ=255)
                                                                                                                           aaaTGTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGa                                                             <  4:545638/110‑1 (MQ=255)
                                                                                                                               gTAGCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAAttt                  >  6:521922/1‑149 (MQ=255)
                                                                                                                                  gCCTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTa                <  1:176013/148‑1 (MQ=255)
                                                                                                                                   ccTGCTGGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATc              <  5:430222/149‑1 (MQ=255)
                                                                                                                                    ctgctgGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCg             >  7:226059/1‑149 (MQ=255)
                                                                                                                                    ctgctgGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCg             >  7:119046/1‑149 (MQ=255)
                                                                                                                                     tgctgGAAAACAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCg             >  6:137485/1‑148 (MQ=255)
                                                                                                                                        tgGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAAt                    <  3:501598/138‑1 (MQ=255)
                                                                                                                                        tgGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAAt                    >  4:501598/1‑138 (MQ=255)
                                                                                                                                             aaCAGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCat                                       >  2:23300/1‑114 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGC  >  NC_000913/3991435‑3991726

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