Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,763,450 A→G 7.7% V187A (GTG→GCG) ‡ 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,4500AG 7.7% 282.0 / 11.2 91V187A (GTG→GCG) ‡sufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base G (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 = 7.68e-01

TCAAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTT  >  NC_000913/1763303‑1763598
                                                                                                                                                   |                                                                                                                                                    
tCAAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTGGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGcacca                                                                                                                                                      >  4:866040/1‑148 (MQ=255)
  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                                                                                                                              >  1:986091/1‑149 (MQ=255)
                       gTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCa                                                                                                                              >  2:54339/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:866040/148‑1 (MQ=255)
                                          gAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGa                                                                                                            <  3:369260/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                                                                             <  1:870276/149‑1 (MQ=255)
                                                                        ggCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                             >  2:517031/1‑149 (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                                                                                >  1:221809/1‑137 (MQ=255)
                                                                                 ataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                                                                                <  2:221809/137‑1 (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                                                             <  3:168902/148‑1 (MQ=255)
                                                                                         aCCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtttt                                                             >  4:805302/1‑148 (MQ=255)
                                                                                              cAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcg                                                          >  4:150146/1‑146 (MQ=255)
                                                                                              cAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcg                                                          <  3:150146/146‑1 (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                                                    <  2:858714/106‑1 (MQ=255)
                                                                                                                                           tcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    >  1:858714/2‑107 (MQ=255)
                                                                                                                                           tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTc                                                                                                     >  2:694279/2‑58 (MQ=255)
                                                                                                                                           tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTc                                                                                                     <  1:694279/57‑1 (MQ=255)
                                                                                                                                           tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc          >  2:636626/2‑149 (MQ=255)
                                                                                                                                           ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGt                              <  1:148184/129‑1 (MQ=255)
                                                                                                                                           ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGt                              >  2:148184/1‑129 (MQ=255)
                                                                                                                                            cGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc          >  3:216604/1‑148 (MQ=255)
                                                                                                                                             ggCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCAc        <  2:669107/149‑1 (MQ=255)
                                                                                                                                                   acGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCttt  <  1:54339/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTT  >  NC_000913/1763303‑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: ‑