Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,638,584 A→T 13.6% F23L (TTT→TTA bamB ← outer membrane protein assembly factor BamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,638,5840AT13.6% 151.6 / 12.7 63F23L (TTT→TTAbamBouter membrane protein assembly factor BamB
Reads supporting (aligned to +/- strand):  ref base A (32/20);  new base T (4/4);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

TTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATC  >  NC_000913/2638444‑2638722
                                                                                                                                            |                                                                                                                                          
ttGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAACGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTaaa                                                                                                                                          <  1:322644/143‑1 (MQ=255)
   tCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCGACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAg                                                                                                                                       >  2:225196/1‑143 (MQ=255)
   tCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGCTAAACAg                                                                                                                                       >  2:1200775/1‑143 (MQ=255)
      gCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGa                                                                                                                                    <  2:1189947/143‑1 (MQ=255)
      gCCAGTGCCGGATGAAGAGTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGa                                                                                                                                    <  1:970662/143‑1 (MQ=255)
       ccAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAACCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTTAACAGCGaa                                                                                                                                   <  2:1003639/143‑1 (MQ=255)
        cAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACCGAAGTGCTCCACGACGTGGTCGGAGTCAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTAGGACATCTTCTTCGCTGTTGACCAGCGAAc                                                                                                                                  >  1:1213342/1‑143 (MQ=255)
          gttccGGATGCAGATTGGAATAGCAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCTGTGGTCGGCGTAAACTGGTTTTCAACGGTTGACAATGGTGACGTCTTTACCACATCTTCTTCGCTGTTTAACAGCGaagaa                                                                                                                                <  2:390839/140‑4 (MQ=255)
                   gAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaa                                                                                                                       >  1:329065/1‑143 (MQ=255)
                    aaGATTGGAAGAGAAGTTGCCAATGGCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTaaaa                                                                                                                      <  1:340300/143‑1 (MQ=255)
                      gATTGGAATAGAAGTTGCCAATGCCGCGACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGAATTTACACGGTTGGGAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                    <  1:480256/143‑1 (MQ=255)
                      gATTGGAATAGAAGTTGCCAATGCCGATACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCCCATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAgg                                                                                                                    <  1:944809/143‑1 (MQ=255)
                       aTTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCGGCTTAAAAggg                                                                                                                   >  2:4964/1‑143 (MQ=255)
                        ttGGAATTGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAGGGGt                                                                                                                  <  2:750701/143‑1 (MQ=255)
                        ttGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGt                                                                                                                  >  1:662193/1‑143 (MQ=255)
                         tGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTa                                                                                                                 >  2:432090/1‑143 (MQ=255)
                              tAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCGATGGGGACATATTTACGACATCTTCTTCGCTGTTAAAAAGCGAACAGCGGCGTAAAAGGCAAACGGa                                                                                                            >  2:891967/1‑143 (MQ=255)
                                 aatttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGGGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAg                                                                                                         <  1:458494/140‑1 (MQ=255)
                                   gTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCa                                                                                                       >  1:1101658/1‑143 (MQ=255)
                                    ttGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGGTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAAAAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAg                                                                                                      >  1:370130/1‑143 (MQ=255)
                                       ccAATGCCGCTACCAACGGAAGTGCTCCCCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGCACATCTTTACCACAGCTTCTTCTCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTcc                                                                                                   <  1:572547/143‑1 (MQ=255)
                                       ccAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTAACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTcc                                                                                                   >  1:643049/1‑143 (MQ=255)
                                               gCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcagca                                                                                           <  1:675176/143‑1 (MQ=255)
                                               gCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGgcagca                                                                                           >  1:845908/1‑143 (MQ=255)
                                                 tACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTCAACAGGGTAACGGAAAGCCGTCCTGGcagcagc                                                                                         >  2:286590/1‑142 (MQ=255)
                                                  aCCAGCGGACGTGCCCGCCGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAAGGGGGACATCTGCACCGCATCTTCTCCGCGGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTa                                                                                        <  1:438029/143‑1 (MQ=255)
                                                    cAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAt                                                                                      >  1:1085456/1‑143 (MQ=255)
                                                         gAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACg                                                                                 >  1:817248/1‑143 (MQ=255)
                                                            gtccgccACGTCGTCGTCGGCGTAAACTGTGTTTCAACGGTTGGCCATGGGGACATCTTTACCACATCTTCTTCGCTGTCAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCaa                                                                              <  2:1134564/138‑1 (MQ=255)
                                                                tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCGGCAGTAACCTATGCAATTGc                                                                          >  2:256770/1‑143 (MQ=255)
                                                                tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGc                                                                          >  1:947059/1‑143 (MQ=255)
                                                                tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGc                                                                          <  1:432090/143‑1 (MQ=255)
                                                                tCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGc                                                                          >  1:779222/1‑143 (MQ=255)
                                                                  cACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCAt                                                                        >  2:361498/1‑143 (MQ=255)
                                                                     gCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACACCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCgg                                                                     >  1:746280/1‑143 (MQ=255)
                                                                         tGGTCGGCGTAAACTGGTTTTCAGCGGTTGGCAACGGGGACATCTTTACCACGTCTTCGTCGCTGTTAGACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTcc                                                                 <  2:672261/143‑1 (MQ=255)
                                                                         tGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTcc                                                                 <  2:496999/143‑1 (MQ=255)
                                                                           gTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCct                                                               >  2:121789/1‑143 (MQ=255)
                                                                           gTCGGCGTAAACTGGTTTTCAACGGCTGGCGATTCCCACAACTTCACCACATCGTCTTCGCTGTTAAACAGCGAACATCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCct                                                               <  1:866126/143‑1 (MQ=255)
                                                                              ggCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTGAACAGCGAACAGCCGCTTAGAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTGACGCAATTGCATCGGGTCCCtctc                                                            <  2:285805/143‑1 (MQ=255)
                                                                                    aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTACACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCGGGCAGCAGTAATGTACGCAATTGCATCGGCTCCCGCTACGATgg                                                      >  2:193305/1‑143 (MQ=255)
                                                                                    aaCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAACGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTACCCCTCAGATgc                                                      >  2:817739/1‑142 (MQ=255)
                                                                                     aCTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACCGCCGCTTAAAAGGGTAACGGAACGCAGTCCGGGCAGCAGTAAGTTGCGCAATTGCATCGGGTCCCTCGAAGAGGGa                                                     >  2:442372/1‑143 (MQ=255)
                                                                                       tGGTTTTCAACGGTTGGCAATGGGGACCTCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTACCGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATCGCACCGGGTCCCTCTCAGACGGACa                                                   >  2:1116872/1‑143 (MQ=255)
                                                                                       tGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACa                                                   >  2:272963/1‑143 (MQ=255)
                                                                                        ggTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAACGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACaa                                                  <  2:631302/143‑1 (MQ=255)
                                                                                                aCGGTTGGCAATGGGGACATCTTTACCACATCTTCTGCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCTGGTCCCTCTCAGATGGACAAATTATTaa                                          >  1:825832/1‑143 (MQ=255)
                                                                                                           tGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTCCCCAATTGCATCCGGTCCCTCGCAGATGGACAAATTATTAATTTTCATCTGc                               <  1:661284/143‑1 (MQ=255)
                                                                                                              ggACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAGTTATTAATTTTCATCTGcatc                            >  1:863874/1‑143 (MQ=255)
                                                                                                                  aTCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTCAAGAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCAttt                        <  1:587688/143‑1 (MQ=255)
                                                                                                                   tCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAACGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTc                       <  2:632350/143‑1 (MQ=255)
                                                                                                                    cTTTACCACATCGTCTGCTCTGTTGAACAGCTAACACCCGCTTAAAGTCGTAACGTAAAGCAGTCCTTGCAGCAGTAATTTACGCAATTTCATCTGGTCCCTCTCACATGTACAAATTATTAATTTTCATCTGCATCATTTCg                      <  2:747722/143‑1 (MQ=255)
                                                                                                                    cTTTACCACAGCTGCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCg                      <  1:104429/143‑1 (MQ=255)
                                                                                                                     tttACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGATTAAAAGGGTAACGGAGAGCAGTCCTGGCAGCAGTACTCTACCCTGGTGCATCGGCGGCCTCTCAGATGGACAAATTATTAATCTTCATCTGAGTCATTTCGc                     >  1:798693/1‑143 (MQ=255)
                                                                                                                        aCCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCa                  <  1:940306/143‑1 (MQ=255)
                                                                                                                           acaTCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCGCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTg               >  1:44454/1‑143 (MQ=255)
                                                                                                                                 tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGCTGGACAAATTATTAATTTTCATCGGCATCATTTCCCCCAGTGCCGGAg         >  2:196220/3‑143 (MQ=255)
                                                                                                                                 tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAg                 <  1:598010/133‑1 (MQ=255)
                                                                                                                                 tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAg                 >  2:598010/3‑135 (MQ=255)
                                                                                                                                 tgttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAAGTTACGCAATTGCATCGGGGCCGTCTCCTATGCGTGAATTATTAATTGTCATCTGCATCATTTAGCTGAGTGCCagag         >  1:131012/3‑143 (MQ=255)
                                                                                                                                   ttcGCTGTTTAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGAGAGATGGCGAAATGTCATCTGCATCATTTCGCTCAGTGCCGGCGTa       >  2:902342/1‑143 (MQ=255)
                                                                                                                                   ttcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTa       >  2:166479/1‑143 (MQ=255)
                                                                                                                                    tcGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTaa      >  1:237048/1‑143 (MQ=255)
                                                                                                                                        tGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATc  >  2:191117/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TTGTCCGCCAGTGCCGGATGAAGATTGGAATAGAAGTTGCCAATGCCGCTACCAACGGAAGTGCTCCACGCCGTGGTCGGCGTAAACTGGTTTTCAACGGTTGGCAATGGGGACATCTTTACCACATCTTCTTCGCTGTTAAACAGCGAACAGCCGCTTAAAAGGGTAACGGAAAGCAGTCCTGGCAGCAGTAATTTACGCAATTGCATCGGGTCCCTCTCAGATGGACAAATTATTAATTTTCATCTGCATCATTTCGCTCAGTGCCGGAGTAACATC  >  NC_000913/2638444‑2638722

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: