Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,146,293 G→A 6.7% A192A (GCC→GCT nicP_9 ← Porin‑like protein NicP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,146,2930GA 6.7% 185.9 / 3.4 60A192A (GCC→GCTnicP_9Porin‑like protein NicP
Reads supporting (aligned to +/- strand):  ref base G (33/23);  new base A (2/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GGCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACCGCGAAAGACCTGAG  >  ACB122/4146159‑4146423
                                                                                                                                      |                                                                                                                                  
ggCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTGAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCggg                                                                                                                                    <  4:21274/135‑1 (MQ=255)
ggCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCgggg                                                                                                                                   >  8:318243/1‑136 (MQ=255)
ggCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCggg                                                                                                                                    <  1:124125/135‑1 (MQ=255)
 gCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGTTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCggg                                                                                                                                    >  6:265816/1‑134 (MQ=255)
 gCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGc                                                                                                                                  <  8:390265/136‑1 (MQ=255)
 gCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTGGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCgggg                                                                                                                                   >  8:168679/1‑135 (MQ=255)
  cTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGc                                                                                                                                  <  2:93177/135‑1 (MQ=255)
  cTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGc                                                                                                                                  <  2:43481/135‑1 (MQ=255)
    tgatgaCTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGc                                                                                                                                  >  1:382058/1‑133 (MQ=255)
     gatgaCTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATa                                                                                                                               >  4:6076/1‑135 (MQ=255)
     gatgaCTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATa                                                                                                                               >  2:372611/1‑135 (MQ=255)
      atgaCTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAg                                                                                                                              <  8:87043/135‑1 (MQ=255)
       tgaCTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACTGGGGCATAGc                                                                                                                             <  4:187033/135‑1 (MQ=255)
       tgaCTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAg                                                                                                                              <  3:60253/134‑1 (MQ=255)
             tAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGc                                                                                                                      >  8:146261/1‑136 (MQ=255)
                        aaGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGAt                                                                                                            >  4:189778/1‑135 (MQ=255)
                            aGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGa                                                                                                        <  2:301228/135‑1 (MQ=255)
                               aTTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAAt                                                                                                     <  6:78512/135‑1 (MQ=255)
                                   ataGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGt                                                                                                 >  5:387459/1‑135 (MQ=255)
                                      gagaTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGaa                                                                                             <  7:375065/136‑1 (MQ=255)
                                        gaTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  8:52301/136‑1 (MQ=255)
                                        gaTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  8:16220/136‑1 (MQ=255)
                                        gaTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGGCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  7:88368/136‑1 (MQ=255)
                                         aTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  3:244864/135‑1 (MQ=255)
                                         aTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  5:339037/135‑1 (MQ=255)
                                         aTCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  5:284751/135‑1 (MQ=255)
                                         aTCTTCAAGTTCGGCGGAGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAAtt                                                                                           <  3:189778/135‑1 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATGGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACg                                                                                        >  3:275440/1‑135 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACgc                                                                                       >  7:182327/1‑136 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACg                                                                                        >  3:303402/1‑135 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGAGCATAGCCAGTGGCtgtg                                                                                                                  >  6:258640/1‑109 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGAGCATAGCCAGTGGCtgtg                                                                                                                  <  5:258640/109‑1 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGAGCATAGCCAGTGGCtgtg                                                                                                                  <  1:282315/109‑1 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGAGCATAGCCAGTGGCtgtg                                                                                                                  >  2:282315/1‑109 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGCCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGAGGTGAGTGATTCGAAATCGGTGGAATTACg                                                                                        >  2:187612/1‑135 (MQ=255)
                                            ttCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGAAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAAAAGGTGGAATTACgc                                                                                       >  7:195031/1‑136 (MQ=255)
                                                 gTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGATATCGGTGGAATTACgctgct                                                                                   >  1:29189/1‑135 (MQ=255)
                                                  ttCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCtt                                                                                  <  3:348870/135‑1 (MQ=255)
                                                                       ggTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGc                                                             <  4:96794/135‑1 (MQ=255)
                                                                        gTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGcc                                                            >  4:308785/1‑135 (MQ=255)
                                                                        gTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGcc                                                            <  2:382058/135‑1 (MQ=255)
                                                                         tCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCg                                                           >  4:237477/1‑135 (MQ=255)
                                                                         tCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCg                                                           >  6:341182/1‑135 (MQ=255)
                                                                           aTGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGtt                                                         >  2:344070/1‑135 (MQ=255)
                                                                              tCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCTGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCg                                                     >  8:308477/1‑136 (MQ=255)
                                                                              tCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTAcc                                                      >  4:345402/1‑135 (MQ=255)
                                                                                cTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCg                                                     >  6:144719/1‑134 (MQ=255)
                                                                                cTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCg                                                     >  4:107622/1‑134 (MQ=255)
                                                                                      aGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAAttt                                              >  6:240311/1‑135 (MQ=255)
                                                                                      aGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTa                                             >  8:180226/1‑136 (MQ=255)
                                                                                                             aGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGc                        >  3:85429/1‑134 (MQ=255)
                                                                                                                  tAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAacca                  >  4:229802/1‑135 (MQ=255)
                                                                                                                      gATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACcg               <  1:372611/134‑1 (MQ=255)
                                                                                                                        tCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACgctgct                                                                                   <  8:66434/64‑1 (MQ=255)
                                                                                                                        tCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACgctgct                                                                                   >  7:66434/1‑64 (MQ=255)
                                                                                                                         ccGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACGACcgcg             >  1:5286/1‑133 (MQ=255)
                                                                                                                          cGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGaa                                                                                             >  6:286456/1‑52 (MQ=255)
                                                                                                                          cGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGaa                                                                                             <  5:286456/52‑1 (MQ=255)
                                                                                                                          cGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACcgcg             >  3:247158/1‑132 (MQ=255)
                                                                                                                          cGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACcgcg             >  4:216031/1‑132 (MQ=255)
                                                                                                                          cGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACcgcg             >  1:397825/1‑132 (MQ=255)
                                                                                                                            ccTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCCCTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGCACCCCCGCGAAAGa        >  5:366290/1‑135 (MQ=255)
                                                                                                                            ccTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCTATATCGTTGACGGTCAGCGAACCACCGCGAAAGa        >  5:77788/1‑135 (MQ=255)
                                                                                                                               cAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACCGCGAAAGACCTg    >  7:93695/1‑136 (MQ=255)
                                                                                                                                aaCCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACCGCGAAAGACCTg    <  1:22372/135‑1 (MQ=255)
                                                                                                                                    ggggCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACCGCGAAAGACCTGAg  >  1:254654/1‑133 (MQ=255)
                                                                                                                                      |                                                                                                                                  
GGCTTGATGACTTTAATGCCGGCGAAGTAGGATTTATAGAGATCTTCAAGTTCGGCGGCGTGCAAACTGAAGGTCATGTCCTTCGCAGGCTTGAAGTCAAATGCCAGGTAGTTGTAGTGATCCGCCTCAACCGGGGCATAGCCAGTGGCTGTGAGTGATTCGAAATCGGTGGAATTACGCTGCTTCACTGCGTTGATTCTTGTCGCCGTTACCGTCAATTTATCGATATCGTTGACGGTCAGCGAACCACCGCGAAAGACCTGAG  >  ACB122/4146159‑4146423

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