Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,638,584 A→T 8.3% F23L (TTT→TTA)  bamB ← BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,638,5840AT 8.3% 198.0 / 9.1 60F23L (TTT→TTA) bamBBamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base T (3/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

AGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCAC  >  NC_000913/2638436‑2638731
                                                                                                                                                    |                                                                                                                                                   
aGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTa                                                                                                                                                     >  1:236772/1‑149 (MQ=255)
  aCAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCACTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCATATCTTCTTCGCTGTTcaa                                                                                                                                                   <  4:186163/149‑1 (MQ=255)
  aCAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                                   >  2:580350/1‑149 (MQ=255)
         tGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGaa                                                                                                                                            <  1:148886/149‑1 (MQ=255)
              gCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGcc                                                                                                                                       <  4:432375/149‑1 (MQ=255)
                cAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGc                                                                                                                                     <  3:192668/149‑1 (MQ=255)
                        gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                             <  3:84356/149‑1 (MQ=255)
                        gATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                             <  2:224122/149‑1 (MQ=255)
                             aGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAAc                                                                                                                        >  3:539066/1‑149 (MQ=255)
                                ttGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGa                                                                                                                     >  2:257818/1‑149 (MQ=255)
                                  ggAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGaaa                                                                                                                   <  1:252358/149‑1 (MQ=255)
                                    aaTAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGc                                                                                                                 >  2:178206/1‑149 (MQ=255)
                                     aTAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                            >  1:266271/1‑121 (MQ=255)
                                     aTAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                            >  2:518818/1‑121 (MQ=255)
                                     aTAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                            <  2:266271/121‑1 (MQ=255)
                                     aTAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                            <  1:518818/121‑1 (MQ=255)
                                      tAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAg                                                                                                               >  1:410435/1‑149 (MQ=255)
                                       aGAAGTTGCCAATACCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGACGTAAACTGGTTTTCAACGGTCGGCAATGGGGAAATCTTTACCACATCTTCTTCGCTGTTAAACATCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGt                                                                                                              <  3:405889/149‑1 (MQ=255)
                                          aGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCt                                                                                                           <  4:179909/149‑1 (MQ=255)
                                           gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTg                                                                                                          >  3:342435/1‑149 (MQ=255)
                                           gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTg                                                                                                          >  4:101918/1‑149 (MQ=255)
                                            ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTgg                                                                                                         <  3:443556/149‑1 (MQ=255)
                                              gCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgca                                                                                                       >  4:331189/1‑149 (MQ=255)
                                                 aaTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCa                                                                                                                >  2:195491/1‑137 (MQ=255)
                                                 aaTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCa                                                                                                                <  1:195491/137‑1 (MQ=255)
                                                 aaTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcagca                                                                                                    >  3:445071/1‑149 (MQ=255)
                                                  aTGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGcagcag                                                                                                   <  2:236772/149‑1 (MQ=255)
                                                      cGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAt                                                                                               >  4:100791/1‑149 (MQ=255)
                                                      cGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAt                                                                                               <  1:141792/149‑1 (MQ=255)
                                                       gCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAtt                                                                                              >  3:213042/1‑149 (MQ=255)
                                                          aCCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTAc                                                                                           <  4:479885/149‑1 (MQ=255)
                                                                  aaGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGc                                                                                   <  1:369574/149‑1 (MQ=255)
                                                                   aGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCATATCTTCTTCGCTGTTAAACAGCGAACAGCGGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCa                                                                                  <  3:499016/149‑1 (MQ=255)
                                                                        tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCggg                                                                             >  3:203500/1‑149 (MQ=255)
                                                                                         gTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACa                                                            >  1:373825/1‑149 (MQ=255)
                                                                                            aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAat                                                         <  4:354562/149‑1 (MQ=255)
                                                                                            aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAat                                                         <  4:441460/149‑1 (MQ=255)
                                                                                            aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAat                                                         <  2:410171/149‑1 (MQ=255)
                                                                                               tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattat                                                      <  4:108530/149‑1 (MQ=255)
                                                                                               tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattat                                                      >  3:382909/1‑149 (MQ=255)
                                                                                               tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattat                                                      >  1:227387/1‑149 (MQ=255)
                                                                                               tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAattat                                                      >  2:311660/1‑149 (MQ=255)
                                                                                                 gTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAAttatta                                                    <  4:233145/149‑1 (MQ=255)
                                                                                                    ttCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAAtt                                                 <  1:178206/149‑1 (MQ=255)
                                                                                                       aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                              >  1:116794/1‑149 (MQ=255)
                                                                                                       aaCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTc                                              >  1:481413/1‑149 (MQ=255)
                                                                                                        aCGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCa                                             <  1:580350/149‑1 (MQ=255)
                                                                                                         cGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCAt                                            <  1:220939/149‑1 (MQ=255)
                                                                                                         cGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCAt                                            <  2:410435/149‑1 (MQ=255)
                                                                                                              ggCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGca                                       <  4:213042/149‑1 (MQ=255)
                                                                                                                  aTGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTCCATCTGcatcat                                   >  3:320263/1‑149 (MQ=255)
                                                                                                                                aCCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCg                     <  1:257818/149‑1 (MQ=255)
                                                                                                                                  cacaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                   >  2:291760/1‑149 (MQ=255)
                                                                                                                                  cacaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                   <  1:80349/149‑1 (MQ=255)
                                                                                                                                  cacaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGa                   <  4:297634/149‑1 (MQ=255)
                                                                                                                                      tcttctTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTaa               >  1:100401/1‑149 (MQ=255)
                                                                                                                                        ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAAc              >  1:515483/1‑148 (MQ=255)
                                                                                                                                        ttcttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAAc              <  2:515483/148‑1 (MQ=255)
                                                                                                                                          cttcGCTGTTAAACAGCAAACAGTCGCTTAGAAGTGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGATAAATTATTAATTTACATCTGCATCATTTCGCTCAGTGCCGGAGTAACATc           <  2:304175/149‑1 (MQ=255)
                                                                                                                                           ttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCg          >  1:268595/1‑149 (MQ=255)
                                                                                                                                             cGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCt        >  1:279260/1‑149 (MQ=255)
                                                                                                                                               cTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCttt      <  3:331189/149‑1 (MQ=255)
                                                                                                                                                   tAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCAc  <  3:224347/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGACAACGTTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATCGCTTTTCAC  >  NC_000913/2638436‑2638731

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