Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 315,244 G→A intergenic (+16/‑47) eaeH → / → insE1 pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913315,2440GA100.0% 266.9 / NA 84intergenic (+16/‑47)eaeH/insE1pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/44);  total (40/44)

TTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAA  >  NC_000913/315111‑315374
                                                                                                                                     |                                                                                                                                  
ttGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGac                                                                                                                               <  2:1183801/139‑1 (MQ=255)
 tGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGaca                                                                                                                              <  1:1360722/139‑1 (MQ=255)
  ggCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             >  2:52170/1‑139 (MQ=255)
     gcAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg                                                                                                                          <  2:1286514/139‑1 (MQ=255)
                      aCAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                                         <  2:871004/139‑1 (MQ=255)
                      aCAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCT‑TTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct                                                                                                        <  1:1361626/139‑1 (MQ=255)
                       cAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct                                                                                                        <  1:215801/139‑1 (MQ=255)
                        aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             <  1:185947/117‑1 (MQ=255)
                        aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             >  2:185947/1‑117 (MQ=255)
                        aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACtctc                                                                                                       >  2:142067/1‑139 (MQ=255)
                          tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                                         <  1:1612538/135‑1 (MQ=255)
                          tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                                         >  2:1612538/1‑135 (MQ=255)
                          tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcag                                                                                                     >  1:17934/1‑139 (MQ=255)
                                ggCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCagag                                                                                               >  1:63625/1‑139 (MQ=255)
                                    aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             <  2:107230/105‑1 (MQ=255)
                                    aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             <  1:1292455/105‑1 (MQ=255)
                                    aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             >  2:1292455/1‑105 (MQ=255)
                                    aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                             >  1:107230/1‑105 (MQ=255)
                                    aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTa                                                                                                              >  2:225288/1‑120 (MQ=255)
                                    aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTa                                                                                                              <  1:225288/120‑1 (MQ=255)
                                         gTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac                                                                                      <  2:63625/139‑1 (MQ=255)
                                              gCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                                 <  1:52170/139‑1 (MQ=255)
                                               cTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACa                                                                                >  2:660146/1‑139 (MQ=255)
                                                    tGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAAc                                                                           >  2:506347/1‑139 (MQ=255)
                                                          aGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATc                                                                     <  1:142067/139‑1 (MQ=255)
                                                                      gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                            >  2:684447/1‑120 (MQ=255)
                                                                      gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                            <  1:684447/120‑1 (MQ=255)
                                                                      gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaaaa                                                         <  2:1626389/139‑1 (MQ=255)
                                                                         aGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGAGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAcc                                                      <  2:844539/139‑1 (MQ=255)
                                                                            cccGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGaca                                                                                                                              <  2:1318914/64‑1 (MQ=255)
                                                                            cccGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGaca                                                                                                                              >  1:1318914/1‑64 (MQ=255)
                                                                              cGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTa                                                 >  2:679744/1‑139 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg                                                                                                                          >  2:1482056/1‑64 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg                                                                                                                          <  1:1482056/64‑1 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGc                                                                                                                        >  1:1568394/1‑66 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGc                                                                                                                        <  2:1568394/66‑1 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                                 <  2:1502548/105‑1 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                                 >  1:1502548/1‑105 (MQ=255)
                                                                                cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa                                               <  2:17934/139‑1 (MQ=255)
                                                                                 aTGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAAc                                              >  1:153488/1‑139 (MQ=255)
                                                                                     tATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCa                                          <  1:660146/139‑1 (MQ=255)
                                                                                       tttCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct                                                                                                        <  1:341216/75‑1 (MQ=255)
                                                                                       tttCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAAatct                                                                                                        >  2:341216/1‑75 (MQ=255)
                                                                                             ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAg                                                                                                                    >  2:1216595/1‑57 (MQ=255)
                                                                                             ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAg                                                                                                                    <  1:1216595/57‑1 (MQ=255)
                                                                                             ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                            <  1:889840/97‑1 (MQ=255)
                                                                                             ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                            >  2:889840/1‑97 (MQ=255)
                                                                                             ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTg                                  >  1:683294/1‑139 (MQ=255)
                                                                                                       cTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAGGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGt                                                  <  2:1087283/113‑1 (MQ=255)
                                                                                                       cTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGt                                                  >  1:1087283/1‑113 (MQ=255)
                                                                                                          taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGc                                                                                                                        <  1:476226/40‑1 (MQ=38)
                                                                                                          taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGc                                                                                                                        >  2:476226/1‑40 (MQ=21)
                                                                                                          taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACt                                                                                         <  2:1377432/71‑1 (MQ=255)
                                                                                                          taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACt                                                                                         >  1:1377432/1‑71 (MQ=255)
                                                                                                          taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGa                     <  1:487868/139‑1 (MQ=255)
                                                                                                          taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGa                     >  2:487868/1‑139 (MQ=255)
                                                                                                            cGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCa                                          <  2:762832/116‑1 (MQ=255)
                                                                                                            cGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCa                                          >  1:762832/1‑116 (MQ=255)
                                                                                                             gCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGc                  >  1:594582/1‑139 (MQ=37)
                                                                                                                aGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATc                                                                     <  1:537622/85‑1 (MQ=255)
                                                                                                                aGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATc                                                                     >  2:537622/1‑85 (MQ=255)
                                                                                                                aGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCt               <  2:153488/139‑1 (MQ=21)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTg                                                                                                                  >  2:305801/1‑39 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTg                                                                                                                  <  1:305801/39‑1 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtca                                                                                                  <  2:865562/55‑1 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtca                                                                                                  >  1:865562/1‑55 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCagag                                                                                               <  1:118202/58‑1 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCagag                                                                                               >  2:118202/1‑58 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac                                                                                      <  2:774741/67‑1 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac                                                                                      >  1:774741/1‑67 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                                    <  2:956646/69‑1 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                                    >  1:956646/1‑69 (MQ=255)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAcc                                                      >  2:106772/1‑99 (MQ=25)
                                                                                                                 gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAcc                                                      <  1:106772/99‑1 (MQ=25)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                            <  2:451597/44‑1 (MQ=255)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                            >  1:451597/1‑44 (MQ=255)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGGGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa             >  2:1654631/1‑139 (MQ=18)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTg                      <  1:1340462/130‑1 (MQ=18)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTg                      >  2:1340462/1‑130 (MQ=18)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa             >  1:1668171/1‑139 (MQ=18)
                                                                                                                  tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa             <  2:1668171/139‑1 (MQ=18)
                                                                                                                   gCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGaa            <  1:506347/139‑1 (MQ=17)
                                                                                                                   gCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGaa            <  2:683294/139‑1 (MQ=17)
                                                                                                                    ccTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAA‑‑‑‑‑‑‑‑‑AAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGaa  >  2:218193/1‑139 (MQ=11)
                                                                                                                                     |                                                                                                                                  
TTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAA  >  NC_000913/315111‑315374

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