Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,151,874 T→A 7.2% intergenic (+23/‑65) acpP → / → fabF acyl carrier protein (ACP)/3‑oxoacyl‑[acyl‑carrier‑protein] synthase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,8740TA 7.2% 280.3 / 9.6 83intergenic (+23/‑65)acpP/fabFacyl carrier protein (ACP)/3‑oxoacyl‑[acyl‑carrier‑protein] synthase II
Reads supporting (aligned to +/- strand):  ref base T (35/42);  new base A (3/3);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

CTTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGCTCTG  >  NC_000913/1151726‑1152022
                                                                                                                                                    |                                                                                                                                                    
ctTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGtt                                                                                                                                                      <  1:586380/149‑1 (MQ=255)
ctTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGtt                                                                                                                                                      <  4:354101/149‑1 (MQ=255)
 ttGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                     <  1:8791/149‑1 (MQ=255)
 ttGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                     <  3:80259/149‑1 (MQ=255)
 ttGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                     <  3:470011/149‑1 (MQ=255)
    acacCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGAc                                                                                                                                                  <  1:741967/149‑1 (MQ=255)
      acCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCg                                                                                                                                                >  3:846967/1‑149 (MQ=255)
      acCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCg                                                                                                                                                >  2:239691/1‑149 (MQ=255)
        cGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGcc                                                                                                                                              <  3:713999/149‑1 (MQ=255)
          ttGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTg                                                                                                                                            >  2:797454/1‑149 (MQ=255)
          ttGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCTg                                                                                                                                            <  2:815310/149‑1 (MQ=255)
            gAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCTgga                                                                                                                                          <  2:774974/149‑3 (MQ=255)
               cTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGttt                                                                                                                                       >  4:329726/1‑149 (MQ=255)
                  gTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTAt                                                                                                                                    <  3:275579/149‑1 (MQ=255)
                    aaTGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCt                                                                                                                                  <  3:296269/149‑1 (MQ=255)
                     aTGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCtt                                                                                                                                 >  3:585431/1‑149 (MQ=255)
                      tGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCttt                                                                                                                                >  2:531274/1‑149 (MQ=255)
                      tGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCttt                                                                                                                                <  3:718590/149‑1 (MQ=255)
                        gCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCttttt                                                                                                                              <  1:794227/149‑1 (MQ=255)
                              gaagaaGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCGAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCa                                                                                                                        <  3:773400/149‑1 (MQ=255)
                                 gaagaGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTa                                                                                                                     >  3:618344/1‑149 (MQ=255)
                                  aagaGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAg                                                                                                                    >  1:12705/1‑149 (MQ=255)
                                    gagTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCt                                                                                                                                             >  1:357430/1‑122 (MQ=255)
                                    gagTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGAACGACCGCCt                                                                                                                                             <  2:357430/122‑1 (MQ=255)
                                       tttGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCa                                                                                                               >  1:825080/1‑149 (MQ=255)
                                            tctcgagaTTCCGGACGAAGAAGCTGACAAAATCCCCACCGTTCAGGCTGCCATTGATTACATCAACGGCCATCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCAGTAGAATCAttttt                                                                                                          <  4:763831/145‑1 (MQ=255)
                                             aCTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAtttttt                                                                                                         <  1:84373/149‑1 (MQ=255)
                                              cTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtc                                                                                                        <  1:134833/149‑1 (MQ=255)
                                              cTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtc                                                                                                        >  2:372472/1‑149 (MQ=255)
                                              cTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtc                                                                                                        <  4:618344/149‑1 (MQ=255)
                                               tGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtcc                                                                                                       >  1:295308/1‑149 (MQ=255)
                                               tGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtcc                                                                                                       >  4:661402/1‑149 (MQ=255)
                                               tGAGATTCCGGACGAAGAAGCCGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtcc                                                                                                       >  1:256740/1‑149 (MQ=255)
                                                gagaTTCCGGAGGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAttttttaga                                                                                                      >  3:763831/1‑146 (MQ=255)
                                                     tCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTccctccct                                                                                                 <  1:835652/149‑1 (MQ=255)
                                                      ccGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTg                                                                                                <  3:565744/149‑1 (MQ=255)
                                                         gACGAAGAACCTGATTCAACCACCCCCGTACAGGCTGCCATTGATTACATCACCGGCCACCAGGCGCAAGTAAACATCTCCAGGCGGTCGCTCGACCGCCTGAGTTCTATCTTTTTGTCCCACTAGAATCATTTTTTGCCTCCCTggag                                                                                             <  3:476591/149‑1 (MQ=255)
                                                                  gCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgt                                                                                    <  2:880812/149‑1 (MQ=255)
                                                                    tGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgtgt                                                                                  >  3:316961/1‑149 (MQ=255)
                                                                    tGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgtgt                                                                                  >  2:877024/1‑149 (MQ=255)
                                                                    tGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgtgt                                                                                  >  1:584083/1‑149 (MQ=255)
                                                                    tGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgt                                                                                    <  3:81282/147‑1 (MQ=255)
                                                                    tGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCGTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgt                                                                                    >  4:81282/1‑147 (MQ=255)
                                                                    tGAGAAAATCACCACCGTTCAGGATGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgtgt                                                                                  >  2:634970/1‑149 (MQ=255)
                                                                      agaAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCt                                                                                <  3:661402/149‑1 (MQ=255)
                                                                        aaaaTCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTaa                                                                              >  2:406003/1‑149 (MQ=255)
                                                                               ccaccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTccctccct                                                                                                 <  4:787886/123‑1 (MQ=255)
                                                                               ccaccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTccctccct                                                                                                 >  3:787886/1‑123 (MQ=255)
                                                                                 accGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCgtgt                                                                     >  2:652316/1‑149 (MQ=255)
                                                                                 accGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCgtgt                                                                     >  1:608173/1‑149 (MQ=255)
                                                                                  ccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtccctccc                                                                                                  >  3:422067/1‑119 (MQ=255)
                                                                                  ccGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTtccctccc                                                                                                  <  4:422067/119‑1 (MQ=255)
                                                                                     ttCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGtt                                                                 <  1:372472/149‑1 (MQ=255)
                                                                                       cAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTtgt                                                               <  1:366097/149‑1 (MQ=255)
                                                                                        aGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTtgtg                                                              >  4:172860/1‑149 (MQ=255)
                                                                                              tcATTGATTACATGAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATGTTATTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGa                                                        >  4:572372/2‑149 (MQ=255)
                                                                                              ccATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAtt                                                                                                             <  2:242354/96‑1 (MQ=255)
                                                                                              ccATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAtt                                                                                                             >  1:242354/1‑96 (MQ=255)
                                                                                                    aTTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGc                                                  >  4:385392/1‑149 (MQ=255)
                                                                                                       aCATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtg                                               >  2:122215/1‑149 (MQ=255)
                                                                                                         aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgtt                                             <  2:608173/149‑1 (MQ=255)
                                                                                                         aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgtt                                             <  1:531274/149‑1 (MQ=255)
                                                                                                         aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgtt                                             <  1:406003/149‑1 (MQ=255)
                                                                                                           cAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgttgt                                           >  4:111303/1‑149 (MQ=255)
                                                                                                           cAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgttgt                                           <  3:145667/149‑1 (MQ=255)
                                                                                                            aaCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtc                                          <  3:572372/149‑1 (MQ=255)
                                                                                                             aCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtct                                         <  2:281150/149‑1 (MQ=255)
                                                                                                             aCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtct                                         >  2:504956/1‑149 (MQ=255)
                                                                                                             aCGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGtct                                         >  2:602308/1‑149 (MQ=255)
                                                                                                                    ccaGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCg                                  <  4:444539/149‑1 (MQ=255)
                                                                                                                        gCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCaa                              <  1:795276/149‑1 (MQ=255)
                                                                                                                        gCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCaa                              >  2:749663/1‑149 (MQ=255)
                                                                                                                               tGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTa                       <  1:122215/149‑1 (MQ=255)
                                                                                                                                 aaCATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTaga                     <  4:311202/149‑1 (MQ=255)
                                                                                                                                 aaCATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTaga                     <  4:74420/149‑1 (MQ=255)
                                                                                                                                 aaCATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTaga                     <  2:797412/149‑1 (MQ=255)
                                                                                                                                     aactcAGGCGGTCGAACGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCt                 >  2:728688/5‑149 (MQ=255)
                                                                                                                                      actcAGGCGGTCGAACGACCGCCTGGGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTa                >  3:767389/4‑149 (MQ=255)
                                                                                                                                       tcCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTAc               >  2:133948/1‑149 (MQ=255)
                                                                                                                                         cAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCt             <  3:629043/149‑1 (MQ=255)
                                                                                                                                         cAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCt             <  4:585431/149‑1 (MQ=255)
                                                                                                                                               gTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAg       <  2:151122/149‑1 (MQ=255)
                                                                                                                                                tCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGc      <  1:239691/149‑1 (MQ=255)
                                                                                                                                                  gTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGaa                                                                                                                  <  2:160359/39‑1 (MQ=255)
                                                                                                                                                  gTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGaa                                                                                                                  >  1:160359/1‑39 (MQ=255)
                                                                                                                                                    tCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGCTCTg  <  2:584083/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGCTCTG  >  NC_000913/1151726‑1152022

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