Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 318,458 G→A 6.6% S37L (TCA→TTA)  rclB ← reactive chlorine species (RCS) stress resistance periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913318,4580GA 6.6% 178.8 / 2.1 65S37L (TCA→TTA) rclBreactive chlorine species (RCS) stress resistance periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TTCCATTGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGAAGTGATTAAGAA  >  NC_000913/318314‑318597
                                                                                                                                                |                                                                                                                                           
ttCCATTGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATCACAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATAc                                                                                                                                         <  1:588981/149‑1 (MQ=255)
ttCCATTGTGAACATCCTTATTTATACAAGACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATAc                                                                                                                                         >  1:487952/1‑149 (MQ=255)
 tCCATTGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTGCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACt                                                                                                                                        <  2:51424/149‑1 (MQ=255)
  ccATTGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTg                                                                                                                                       <  1:552573/149‑1 (MQ=255)
    aTTGTGAACATCCTTTTTTTTACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAAGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGAc                                                                                                                                     <  2:257058/149‑1 (MQ=255)
    aTTGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGAc                                                                                                                                     >  2:454419/1‑149 (MQ=255)
     ttGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACg                                                                                                                                    >  1:211853/1‑149 (MQ=255)
         gaagagccTTATTTATGTAATACCGCTGTTTCCGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTTTATTTTCTGGCTTTTTCAGCCAGGGCATCTTCCAGAACTTTCGGTGATTCTGATACTGACGAGAc                                                                                                                                <  2:712341/143‑1 (MQ=255)
         gAACATCCTTATTTATACAATACCGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGAc                                                                                                                                <  1:218620/149‑1 (MQ=255)
             aTCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGATTTTTAAGTATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTg                                                                                                                            >  1:385935/1‑149 (MQ=255)
             aTCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTg                                                                                                                            <  1:462550/149‑1 (MQ=255)
                 ttatttatACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGatat                                                                                                                        >  1:46626/1‑149 (MQ=255)
                  tatttatACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATg                                                                                                                       <  1:209620/149‑1 (MQ=255)
                   atttatACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATAAGCAACTGAGGCCTCCAGAACTGTCGGTGATGCTCATACTGACGAGACGCTGaaaggg                                                                                                                      >  1:350187/1‑144 (MQ=255)
                     ttatACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGttt                                                                                                                    >  1:26395/1‑149 (MQ=255)
                       ataCAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGCGACGCTGATATGTTTCa                                                                                                                  >  2:51741/1‑149 (MQ=255)
                       ataCAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGATGATATGTTTCa                                                                                                                  >  2:111127/1‑149 (MQ=255)
                         accaTACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGa                                                                                                                <  1:151298/146‑1 (MQ=255)
                         aCAATACTGCGGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGa                                                                                                                <  2:558312/149‑1 (MQ=255)
                           aaTACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGCGACGCTGATATGTTTCAGAat                                                                                                              >  1:16529/1‑149 (MQ=255)
                                 gctgTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTAtttt                                                                                                        <  1:132760/149‑1 (MQ=255)
                                   tgTTGCGGTTGAATTATTATCAGTTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTCCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTat                                                                                                      <  1:68811/149‑1 (MQ=255)
                                   tgTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTat                                                                                                      <  1:454419/149‑1 (MQ=255)
                                         ggTTGAATTATTATCAATTCGCATCGAAGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTCATGCTGATACTGACGCGACGCAGAAATGTTTCAGAATTATTATATAATCTg                                                                                                >  2:15708/1‑149 (MQ=255)
                                                    tatCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTCTATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGaa                                                                                     >  2:437908/1‑149 (MQ=255)
                                                    tatCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGaa                                                                                     >  2:496324/1‑149 (MQ=255)
                                                    tatCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTAATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGaa                                                                                     >  2:119902/1‑149 (MQ=255)
                                                    tatCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTAATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGaa                                                                                     >  2:513107/1‑149 (MQ=255)
                                                     atCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGaaa                                                                                    >  2:505641/1‑149 (MQ=255)
                                                        aaTTCGCATCGATGTGACTTTCCAGGATGAAGCATTATTTATGCGGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATg                                                                                 >  1:193046/1‑149 (MQ=255)
                                                        aaTTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATg                                                                                 >  1:539904/1‑149 (MQ=255)
                                                        aaTTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAAGAGAATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCATAATTATTTTATCATCTGCATTGGTTGAAAATg                                                                                 >  2:441194/1‑149 (MQ=255)
                                                                 cGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGc                                                                        <  2:26395/149‑1 (MQ=255)
                                                                                 gATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGc                                                        >  2:101057/1‑149 (MQ=255)
                                                                                      aGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATCGTGTTGCAAATa                                                   >  2:151483/1‑149 (MQ=255)
                                                                                      aGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGACATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATa                                                   >  2:418948/1‑149 (MQ=255)
                                                                                      aGCATTATATTTTCGGGCTATATCAGCAATGGTATCGTCCAGTACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATCTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAAGAGTGTTGCAAATa                                                   <  1:564707/149‑1 (MQ=255)
                                                                                          ttATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTTTTGCAAATAAAAc                                               <  2:555272/149‑1 (MQ=255)
                                                                                                tttCTGGCTATATCAGAAATGGTAGCCTGCAGACATGGCGGTGATGCTGATACTGACGAGACGCTGATATGGTTCAGAACTATATTATGCTCCTCATACGGGGTCAAGGGCGTGACGCCAGATAATAGTGTTGCAAATAAAACAGAttt                                         <  2:769393/149‑1 (MQ=255)
                                                                                                 ttCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGAtttt                                        >  2:59269/1‑149 (MQ=255)
                                                                                                  tCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGAttttt                                       <  1:497237/149‑1 (MQ=255)
                                                                                                      gtttgagcAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTaaa                                   <  2:211853/142‑1 (MQ=255)
                                                                                                         atatCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACat                                <  1:59269/149‑1 (MQ=255)
                                                                                                           atCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCGGACATATTTCAGAATTATTTAAACATCTGAAATGATGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACagat                              <  2:385935/149‑1 (MQ=255)
                                                                                                            tCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACAtata                             <  2:567090/149‑1 (MQ=255)
                                                                                                            tCAGCAATGGTATCCTCCAGAACTGTCGGGGATGCTGATACTGACGAGACGTTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACAtata                             <  2:388572/149‑1 (MQ=255)
                                                                                                            tCAGCAATCGTATCCCCCAGAACTGTAGGAGATACACATACCCACGATACGCTGATAAGTTACATAATTAGTTTAGCATCTACATTCGATGAAAATTCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACAtata                             <  2:466020/149‑1 (MQ=255)
                                                                                                             cAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTTATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACATATTTTTTAAACATATaa                            <  2:273060/149‑1 (MQ=255)
                                                                                                             cAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATaa                            >  1:8942/1‑149 (MQ=255)
                                                                                                             cAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAACAAAACAGATTTTTTAAACATATaa                            >  1:138005/1‑149 (MQ=255)
                                                                                                               tcAATTGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATaaaa                          <  1:230216/148‑1 (MQ=255)
                                                                                                               gCAATGGTATCCTCCAGAACTGTCGTTGATGCTTATACTGACGAGTTGTTGAAATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTCTTTAAACATATaaaa                          <  1:306123/149‑1 (MQ=255)
                                                                                                                   tGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATAcc                      >  1:127603/1‑149 (MQ=255)
                                                                                                                         cctccAGAACTGTCGGTGATGCTAATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAAt                <  1:513107/149‑1 (MQ=255)
                                                                                                                         cctccAGAACTGTCGGTGATGCTAATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAAt                <  1:119902/149‑1 (MQ=255)
                                                                                                                          ctccAGAACTGTCGGTGATGCTGATACTGACGAGACGGTGATATGTTTCAGAATTATTTTATCGTCTGCATTGGTGGAAAATGCCATATCGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATg               >  1:517566/1‑149 (MQ=255)
                                                                                                                          ctccAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATg               >  2:375385/1‑149 (MQ=255)
                                                                                                                           tccAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGa              <  1:749934/149‑1 (MQ=255)
                                                                                                                           tccAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGa              <  1:667491/149‑1 (MQ=255)
                                                                                                                             cAGAACTGTCGGTGATGCTGATAATAACGAGACGCTGATATGTTTAAGACTAATTTTATCATCTGCATTGGTGGAAAATGCAATAACGCCAGATAATAGTGTTGCACACAAAACAGATTTTTAAAACATATAAAATACAATCAATGAAg            >  2:10206/1‑149 (MQ=255)
                                                                                                                              aGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGAAGt           >  1:142990/1‑149 (MQ=255)
                                                                                                                                 aCTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGAAGTGAt        <  2:777597/149‑1 (MQ=255)
                                                                                                                                  cTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGAAGTGAtt       <  1:101057/149‑1 (MQ=255)
                                                                                                                                     tCGGTGATGCTGATACTGACGAGACGCTAAAAAGCTTAATTAGTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAAGGAAGTGATTaag    >  1:193883/1‑149 (MQ=255)
                                                                                                                                       ggTGATGCTTATACTGACGAGGCGCTGAGATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGAAGTGATTaagaa  <  2:90855/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
TTCCATTGTGAACATCCTTATTTATACAATACTGCTGTTGCGGTTGAATTATTATCAATTCGCATCGATGTGACTTTCCAGGATGAAGCATTATATTTTCTGGCTATATCAGCAATGGTATCCTCCAGAACTGTCGGTGATGCTGATACTGACGAGACGCTGATATGTTTCAGAATTATTTTATCATCTGCATTGGTGGAAAATGCCATAACGCCAGATAATAGTGTTGCAAATAAAACAGATTTTTTAAACATATAAAATACCATCAATGAAGTGATTAAGAA  >  NC_000913/318314‑318597

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: