Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,875,275 G→A 5.5% H49Y (CAC→TAC)  cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,875,2750GA 5.5% 244.1 / 3.7 73H49Y (CAC→TAC) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base G (27/42);  new base A (2/2);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-02

CCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTT  >  NC_000913/2875128‑2875422
                                                                                                                                                   |                                                                                                                                                   
ccACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgc                                                                                                                                                    >  2:789940/1‑149 (MQ=255)
 cACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgca                                                                                                                                                   >  1:103428/1‑149 (MQ=255)
 cACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgca                                                                                                                                                   >  4:530075/1‑149 (MQ=255)
 cACGTCAACGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgca                                                                                                                                                   >  1:25808/1‑149 (MQ=255)
    gTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgcagca                                                                                                                                                >  1:312783/1‑149 (MQ=255)
      cAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgca                                                                                                                                                   >  1:53910/1‑144 (MQ=255)
      cAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgca                                                                                                                                                   <  2:53910/144‑1 (MQ=255)
          ggTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAgacgac                                                                                                                                          <  3:530075/149‑1 (MQ=255)
          ggTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAgacgac                                                                                                                                          <  1:363841/149‑1 (MQ=255)
             gATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCaa                                                                                                                                       <  1:800520/149‑1 (MQ=255)
             gATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCaa                                                                                                                                       <  2:25808/149‑1 (MQ=255)
               tGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATc                                                                                                                                     <  3:711266/149‑1 (MQ=255)
               tGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATc                                                                                                                                     <  2:261160/149‑1 (MQ=255)
                  ccTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga                                                                                                                                  <  4:382772/149‑1 (MQ=255)
                  ccTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga                                                                                                                                  <  1:249946/149‑1 (MQ=255)
                  ccTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga                                                                                                                                  >  4:105139/1‑149 (MQ=255)
                  ccTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga                                                                                                                                  <  4:653273/149‑1 (MQ=255)
                            cTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGc                                                                                                                        <  2:372560/149‑1 (MQ=255)
                                 cTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCa                                                                                                                   >  3:462125/1‑149 (MQ=255)
                                      aGGCCGTCCACCAGCAGAGTCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGAc                                                                                                              <  1:440507/149‑1 (MQ=255)
                                      aGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGAc                                                                                                              >  1:58453/1‑149 (MQ=255)
                                        gCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAga                                                                                                                                              >  4:480368/1‑115 (MQ=255)
                                        gCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAga                                                                                                                                              <  3:480368/115‑1 (MQ=255)
                                        gCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGc                                                                                                            <  4:35982/149‑1 (MQ=255)
                                           gTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                                                                                         <  3:105139/149‑1 (MQ=255)
                                             ccaccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACca                                                                                                       <  2:150916/149‑1 (MQ=255)
                                              caccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACcac                                                                                                      >  1:328499/1‑149 (MQ=255)
                                               accaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACcaca                                                                                                     >  3:436191/1‑149 (MQ=255)
                                                      gagCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGa                                                                                              >  3:275798/1‑149 (MQ=255)
                                                         ccAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTACAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                                                                                             >  1:699671/1‑131 (MQ=255)
                                                         ccAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTACAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                                                                                             <  2:699671/131‑1 (MQ=255)
                                                         ccAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTACAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                                                                                             >  1:806211/1‑131 (MQ=255)
                                                         ccAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTACAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACg                                                                                                             <  2:806211/131‑1 (MQ=255)
                                                            gATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga                                                                                                                                  >  2:304065/1‑107 (MQ=255)
                                                            gATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga                                                                                                                                  <  1:304065/107‑1 (MQ=255)
                                                             aTCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgca                                                                                       >  2:562638/1‑149 (MQ=255)
                                                                 aGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTTATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGcagcag                                                                                   >  1:710616/1‑149 (MQ=255)
                                                                 aGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGCGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCAt                                                                                                                  <  4:675833/118‑1 (MQ=255)
                                                                 aGCTTTTCGCCCTGGGTGCCGTGAAGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCAt                                                                                                                  >  3:675833/1‑118 (MQ=255)
                                                                     tttCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCtt                                                                               <  4:436191/149‑1 (MQ=255)
                                                                         gCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTAt                                                                           <  4:787904/149‑1 (MQ=255)
                                                                            cTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                                                                                         <  4:771455/116‑1 (MQ=255)
                                                                            cTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                                                                                         >  3:771455/1‑116 (MQ=255)
                                                                               ggTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATgttgtt                                                                     <  4:275798/149‑1 (MQ=255)
                                                                                 tGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTgc                                                                   <  1:816572/149‑1 (MQ=255)
                                                                                    cGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                                                                                         >  3:175540/1‑108 (MQ=255)
                                                                                    cGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAc                                                                                                         <  4:175540/108‑1 (MQ=255)
                                                                                    cGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCa                                                                <  2:205067/149‑1 (MQ=255)
                                                                                     gTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCaa                                                               <  3:703195/149‑1 (MQ=255)
                                                                                     gTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCaa                                                               <  1:786430/149‑1 (MQ=255)
                                                                                            ttACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                                                                                               <  3:107509/94‑1 (MQ=255)
                                                                                            ttACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGa                                                                                                               >  4:107509/1‑94 (MQ=255)
                                                                                                tGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTa                                                                                                          <  2:348852/95‑1 (MQ=255)
                                                                                                tGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTa                                                                                                          >  1:348852/1‑95 (MQ=255)
                                                                                                tGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgc                                                                                        >  3:732603/1‑113 (MQ=255)
                                                                                                tGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgc                                                                                        <  4:732603/113‑1 (MQ=255)
                                                                                                tGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCt                                                    >  4:323309/1‑149 (MQ=255)
                                                                                                  tNGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTc                                                  >  3:204704/1‑149 (MQ=255)
                                                                                                    gTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGa                                                >  2:248612/1‑149 (MQ=255)
                                                                                                       aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGAcgc                                             <  2:312783/149‑1 (MQ=255)
                                                                                                       aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGAcgc                                             <  2:35045/149‑1 (MQ=255)
                                                                                                            gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGTCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCt                                        <  1:202869/149‑1 (MQ=255)
                                                                                                              gATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTc                                      >  3:673268/1‑149 (MQ=255)
                                                                                                              gATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTc                                      >  1:131763/1‑149 (MQ=255)
                                                                                                                  agagCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTg                                  <  2:568086/149‑1 (MQ=255)
                                                                                                                  agagCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTg                                  <  2:522200/149‑1 (MQ=255)
                                                                                                                   gagCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTgg                                 >  4:281009/1‑149 (MQ=255)
                                                                                                                        gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGAtt                            >  2:418005/1‑149 (MQ=255)
                                                                                                                        gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGAtt                            <  2:131763/149‑1 (MQ=255)
                                                                                                                        gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGAtt                            <  1:562638/149‑1 (MQ=255)
                                                                                                                         aTCTTCGTAGAGTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGAttt                           >  3:693497/1‑149 (MQ=255)
                                                                                                                                   aTTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAg                 >  1:101146/1‑149 (MQ=255)
                                                                                                                                    tttGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGt                <  2:213932/149‑1 (MQ=255)
                                                                                                                                            gTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGtt        >  2:9438/1‑149 (MQ=255)
                                                                                                                                             tATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCGCCAGGCTTTTAGTTGGTTGGGCAATCATCCAGCCTTGACCGCCGCCTCCTTTGCCGTTTGGTGTTGCAATTGCgggtttc       >  1:438245/1‑144 (MQ=255)
                                                                                                                                             tATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTCGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTc       <  4:462125/149‑1 (MQ=255)
                                                                                                                                               tCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCAt     <  1:437233/149‑1 (MQ=255)
                                                                                                                                                cGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATc    <  4:794411/149‑1 (MQ=255)
                                                                                                                                                 gTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCt   <  2:58453/149‑1 (MQ=255)
                                                                                                                                                  tGCAGCAAACGAACCATCAGAGTTATTTTGCCGTCATCGACCGTTACACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCtt  <  2:438245/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTT  >  NC_000913/2875128‑2875422

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