Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,763,451 C→T 7.8% V187M (GTG→ATG) ‡ sufB ← component of SufBCD Fe‑S cluster assembly scaffold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,763,4510CT 7.8% 283.4 / 10.9 91V187M (GTG→ATG) ‡sufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base T (3/4);  total (45/46)
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.93e-01

AAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTT  >  NC_000913/1763305‑1763598
                                                                                                                                                  |                                                                                                                                                   
aaaCTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACgg                                                                                                                                                   >  1:341692/1‑149 (MQ=255)
 aaCTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATCAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACgg                                                                                                                                                   <  3:737091/148‑1 (MQ=255)
 aaCTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACgg                                                                                                                                                   <  3:930405/148‑1 (MQ=255)
   cTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGc                                                                                                                                                <  2:247789/149‑1 (MQ=255)
    tGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGcc                                                                                                                                               >  1:955968/1‑149 (MQ=255)
     gCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGcc                                                                                                                                               >  3:168598/1‑148 (MQ=255)
      cccGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCg                                                                                                                                              <  3:152869/148‑1 (MQ=255)
          gTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGGACCACGGTGCCGAGat                                                                                                                                          >  4:873691/1‑148 (MQ=255)
          gTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCgggga                                                                                                                                          <  4:764446/148‑5 (MQ=255)
              tttCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATAttt                                                                                                                                      <  3:467144/148‑1 (MQ=255)
                   gCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCa                                                                                                                                 >  3:175259/1‑148 (MQ=255)
                   gCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCa                                                                                                                                 >  4:15441/1‑148 (MQ=255)
                     gTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCa                                                                                                                              >  2:54339/1‑149 (MQ=255)
                     gTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCa                                                                                                                              >  1:986091/1‑149 (MQ=255)
                     gTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCACATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACTGTGCCGAGATATTTACGCACCa                                                                                                                              <  1:652049/149‑1 (MQ=255)
                          tGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCACCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                                                                         >  2:405050/1‑149 (MQ=255)
                            cGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCgg                                                                                                                       <  1:477129/149‑1 (MQ=255)
                            cGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCgg                                                                                                                       <  1:872212/149‑1 (MQ=255)
                            cGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCgg                                                                                                                       <  1:980361/149‑1 (MQ=255)
                             gAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCggg                                                                                                                      >  2:120639/1‑149 (MQ=255)
                                aaTAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCA‑AGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                                                                                  <  2:347128/149‑1 (MQ=255)
                                 aTAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                                                                                  <  1:395648/149‑1 (MQ=255)
                                       ggAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTgg                                                                                                             <  3:133558/148‑1 (MQ=255)
                                        gAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGa                                                                                                            <  3:369260/148‑1 (MQ=255)
                                        gAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGa                                                                                                            <  3:866040/148‑1 (MQ=255)
                                         aaaGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATc                                                                                                          <  1:405050/149‑1 (MQ=255)
                                             ttCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCt                                                                                                      >  1:2471/1‑149 (MQ=255)
                                              tCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCt                                                                                                      >  3:67761/1‑148 (MQ=255)
                                              tCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCt                                                                                                      >  4:231623/1‑148 (MQ=255)
                                                cATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCa                                                                                                    <  4:483748/148‑1 (MQ=255)
                                                 aTCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCAc                                                                                                   >  3:436386/1‑148 (MQ=255)
                                                  tCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCAcc                                                                                                  >  4:797046/1‑148 (MQ=255)
                                                    gggCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGGACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaa                                                                                                <  3:873691/148‑1 (MQ=255)
                                                     ggCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaaa                                                                                               >  4:487316/1‑148 (MQ=255)
                                                     ggCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaaa                                                                                               <  3:78976/148‑1 (MQ=255)
                                                     ggCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                                                                              <  1:204164/149‑1 (MQ=255)
                                                          cgcACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                                                                          <  3:925131/148‑1 (MQ=255)
                                                             aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaa                                                                                      >  1:291777/1‑149 (MQ=255)
                                                                 cTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaata                                                                                  >  2:809497/1‑149 (MQ=255)
                                                                  tttAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                                                                                 <  2:54906/149‑1 (MQ=255)
                                                                  tttAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaata                                                                                  >  3:157731/1‑148 (MQ=255)
                                                                      ggCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                             >  2:517031/1‑149 (MQ=255)
                                                                      ggCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                             <  1:870276/149‑1 (MQ=255)
                                                                      ggCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTc                                                                              <  4:67761/148‑1 (MQ=255)
                                                                            taaataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTc                                                                              >  3:775850/1‑142 (MQ=255)
                                                                            taaataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCt                                                                        >  3:430161/1‑148 (MQ=255)
                                                                            taaataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGACATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTc                                                                              <  4:775850/142‑1 (MQ=255)
                                                                              aataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTCGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAATCCCTGCTCCg                                                                     >  2:205625/1‑149 (MQ=255)
                                                                               ataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                                                                                <  2:221809/137‑1 (MQ=255)
                                                                               ataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                                                                                >  1:221809/1‑137 (MQ=255)
                                                                                taaaCGTACCGTCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGc                                                                    >  3:602631/1‑148 (MQ=255)
                                                                                taaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTc                                                                       <  4:102302/145‑1 (MQ=255)
                                                                                taaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTc                                                                       >  3:102302/1‑145 (MQ=255)
                                                                                taaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGc                                                                    >  4:168902/1‑148 (MQ=255)
                                                                                  aaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAg                                                                 >  2:192128/1‑149 (MQ=255)
                                                                                       aCCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtttt                                                             >  4:805302/1‑148 (MQ=255)
                                                                                       aCCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtttt                                                             <  3:168902/148‑1 (MQ=255)
                                                                                            cAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcg                                                          <  3:150146/146‑1 (MQ=255)
                                                                                            cAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcg                                                          >  4:150146/1‑146 (MQ=255)
                                                                                              gagGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATa                                                     >  1:310917/1‑149 (MQ=255)
                                                                                               agGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACTCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    <  2:936784/149‑1 (MQ=255)
                                                                                                     ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagtagt                                               <  3:762813/148‑1 (MQ=255)
                                                                                                      cgccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTgg                                             <  1:616768/149‑1 (MQ=255)
                                                                                                       gccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGc                                            <  2:986091/149‑1 (MQ=255)
                                                                                                        ccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGc                                            >  3:695208/1‑148 (MQ=255)
                                                                                                          gcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAAc                                         <  1:350173/149‑1 (MQ=255)
                                                                                                                aGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGaaa                                    <  4:695208/148‑1 (MQ=255)
                                                                                                                   gcgGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACtaa                                <  2:856421/149‑3 (MQ=255)
                                                                                                                        aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGa                           >  1:129100/1‑149 (MQ=255)
                                                                                                                        aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGa                           >  1:669107/1‑149 (MQ=255)
                                                                                                                          agaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaa                                                                                                >  3:1040875/1‑78 (MQ=255)
                                                                                                                          agaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaa                                                                                                <  4:1040875/78‑1 (MQ=255)
                                                                                                                          agaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaa                          >  4:337813/1‑148 (MQ=255)
                                                                                                                           gaagTAGTCATGCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAACGTGGATCGCCTCACCAAAGGAACAGCAAATAATTCCCTGCTCCGCCAGTTTTTCGCAATAAGTAGAGGCACCCGAAACTGAGTCGaaaa                        >  2:1020344/1‑149 (MQ=255)
                                                                                                                             agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAAt                       <  3:797046/148‑1 (MQ=255)
                                                                                                                                tGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGc                    <  3:940194/148‑1 (MQ=255)
                                                                                                                                  tCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCAt                  >  3:623129/1‑148 (MQ=255)
                                                                                                                                  tCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATc                 <  2:310917/149‑1 (MQ=255)
                                                                                                                                   cATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATcc                <  1:863476/149‑1 (MQ=255)
                                                                                                                                    aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCa               >  2:1049608/1‑149 (MQ=255)
                                                                                                                                     ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCa               <  4:175259/148‑1 (MQ=255)
                                                                                                                                        cccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCg            >  3:972447/1‑148 (MQ=255)
                                                                                                                                         tcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    >  1:858714/2‑107 (MQ=255)
                                                                                                                                         tcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    <  2:858714/106‑1 (MQ=255)
                                                                                                                                         tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTc                                                                                                     <  1:694279/57‑1 (MQ=255)
                                                                                                                                         tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTc                                                                                                     >  2:694279/2‑58 (MQ=255)
                                                                                                                                         tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc          >  2:636626/2‑149 (MQ=255)
                                                                                                                                         ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGt                              >  2:148184/1‑129 (MQ=255)
                                                                                                                                         ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGt                              <  1:148184/129‑1 (MQ=255)
                                                                                                                                          cGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc          >  3:216604/1‑148 (MQ=255)
                                                                                                                                           ggCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCAc        <  2:669107/149‑1 (MQ=255)
                                                                                                                                                 acGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCttt  <  1:54339/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTT  >  NC_000913/1763305‑1763598

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