Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,580,442 C→T 6.4% W452* (TGG→TGA)  parE ← DNA topoisomerase 4 subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,580,4420CT 6.4% 191.1 / 2.5 62W452* (TGG→TGA) parEDNA topoisomerase 4 subunit B
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base T (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCAC  >  ACB122/5580308‑5580575
                                                                                                                                      |                                                                                                                                     
gAGTCAGCGTCGGCTAGGATGCAGATCTTGCCATAGCGCATCTGCGCAAGGTCGGTGGCACCCGGATCGACGACAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTcc                                                                                                                                       <  6:55202/135‑1 (MQ=255)
      gCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGgt                                                                                                                                  <  6:44423/134‑1 (MQ=255)
       cGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGgtgt                                                                                                                                <  1:243364/135‑1 (MQ=255)
        gTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTccc                                                                                                                                      <  8:39786/128‑1 (MQ=255)
        gTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTccc                                                                                                                                      >  7:39786/1‑128 (MQ=255)
        gTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGgtgt                                                                                                                                >  5:132232/1‑134 (MQ=255)
        gTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGgtgt                                                                                                                                >  2:243364/1‑134 (MQ=255)
          cGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGa                                                                                                                             >  4:343248/1‑135 (MQ=255)
                    gCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCTCAGGTGTTGAGGATCTTGcc                                                                                                                   >  8:234755/1‑135 (MQ=255)
                    gCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCTCAGGTGTTGAGGATCTTGcc                                                                                                                   >  8:240650/1‑135 (MQ=255)
                       gATCTTGCAATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCcgcg                                                                                                                <  5:145200/135‑1 (MQ=255)
                                aTAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAg                                                                                                       >  8:326648/1‑135 (MQ=255)
                                     gcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATgg                                                                                                  >  7:169129/1‑135 (MQ=255)
                                     gcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATgg                                                                                                  >  2:293263/1‑135 (MQ=255)
                                       aGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGcc                                                                                                >  6:379979/1‑135 (MQ=255)
                                       aGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGcc                                                                                                >  8:364901/1‑135 (MQ=255)
                                        gCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCTCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                               <  7:240650/135‑1 (MQ=255)
                                        gCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCTCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                               <  7:234755/135‑1 (MQ=255)
                                        gCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                               <  6:339707/135‑1 (MQ=255)
                                            cgcAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGaa                                                                                           >  4:294283/1‑135 (MQ=255)
                                                 ggTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAGTTCCt                                                                                      >  5:36482/1‑135 (MQ=255)
                                                        tGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCcg                                                                               <  3:133386/135‑1 (MQ=255)
                                                         gCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCcgc                                                                              >  3:225198/1‑135 (MQ=255)
                                                         gCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGAATTCGCCCCCGTCCACTTCCCAGGAGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCcgc                                                                              >  1:222138/1‑135 (MQ=255)
                                                            cccGGCTCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGcc                                                                           <  1:331266/135‑1 (MQ=255)
                                                            cccGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGGGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGcc                                                                           <  3:114006/135‑1 (MQ=255)
                                                              cGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTg                                                                         <  7:77701/135‑1 (MQ=255)
                                                              cGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTg                                                                         >  6:159690/1‑135 (MQ=255)
                                                                  tCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTg                                                                     >  8:206254/1‑135 (MQ=255)
                                                                    gACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                                                   >  8:92458/1‑135 (MQ=255)
                                                                        ccAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGaa                                                               >  6:130827/1‑135 (MQ=255)
                                                                             ggCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccacc                                                          <  6:327560/135‑1 (MQ=255)
                                                                                   ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCgc                                                    >  2:194291/1‑134 (MQ=255)
                                                                                      aaTGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTc                                                 <  7:206254/135‑1 (MQ=255)
                                                                                      aaTGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTc                                                 <  3:267370/135‑1 (MQ=255)
                                                                                      aaTGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGACGCGCAGCCGCAGGATGGCCTGGAATTCCATGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTc                                                 >  8:223592/1‑135 (MQ=255)
                                                                                      aaTGTTGTGCACTTCCTGGCTGGACAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTc                                                 >  7:11403/1‑135 (MQ=255)
                                                                                           tgtgCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCt                                            <  5:379979/135‑1 (MQ=255)
                                                                                           tgtgCACTGCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCt                                            <  7:104299/135‑1 (MQ=255)
                                                                                              gCACTTCCTGGCTGGCCAGGACTTCGCCCCCGCCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcg                                          <  7:14388/134‑1 (MQ=255)
                                                                                                aCTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGCCCCCCCGGCCGAGTCACCCTcga                                         >  5:374446/1‑133 (MQ=255)
                                                                                                aCTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGCACCACCGGCCGAGTCACCCTcga                                         >  1:216707/1‑133 (MQ=255)
                                                                                                aCTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcga                                         >  3:285046/1‑133 (MQ=255)
                                                                                                aCTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcga                                         >  4:153832/1‑133 (MQ=255)
                                                                                                       ggctggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaacc                                >  7:356040/1‑134 (MQ=255)
                                                                                                       ggctggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaacc                                >  6:245785/1‑134 (MQ=255)
                                                                                                       ggctggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaaca                                >  7:88121/1‑135 (MQ=255)
                                                                                                          tggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGCACAACt                             >  3:33943/1‑135 (MQ=255)
                                                                                                            gcCCGGACTTCGCCCCCGTCCACTTCCCAGGTGTGGAGGATCGTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTTCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCg                           <  7:276505/135‑1 (MQ=255)
                                                                                                               aGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCa                        >  3:60951/1‑135 (MQ=255)
                                                                                                                ggACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAc                       <  3:252859/135‑1 (MQ=255)
                                                                                                                  cctgcGCCCCCGTCCAGTTCCCCGGTGTTGTGGATCTCGCCGCGCAGCGGCTGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     <  2:358473/131‑1 (MQ=255)
                                                                                                                  aCTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     <  2:219249/135‑1 (MQ=255)
                                                                                                                  aCTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     <  5:196652/135‑1 (MQ=255)
                                                                                                                  aCTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     <  7:92458/135‑1 (MQ=255)
                                                                                                                  aCTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCACGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     <  4:60951/135‑1 (MQ=255)
                                                                                                                       gCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCgg                <  4:225198/135‑1 (MQ=255)
                                                                                                                        cccccGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCgg                <  1:153484/134‑1 (MQ=255)
                                                                                                                         ccccGTCCACATCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTCTCCCGCGCCTGCTTGGCCGAACCACCGCCCGAGTCACCCTCGCCGAGGAACACCTCGGCACGCATCGGGt              >  5:96176/1‑135 (MQ=255)
                                                                                                                            cGTCCACTTCCCAGGTGTTGTGGATCTTGCCGCGCCGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCt           <  2:193509/135‑1 (MQ=255)
                                                                                                                              tCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTgg         >  7:117853/1‑135 (MQ=255)
                                                                                                                                 aCTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     >  3:315129/1‑120 (MQ=255)
                                                                                                                                 aCTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGc                     <  4:315129/120‑1 (MQ=255)
                                                                                                                                  cTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTg     <  3:65154/135‑1 (MQ=255)
                                                                                                                                  cTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTg     <  7:200212/135‑1 (MQ=255)
                                                                                                                                      ccAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGcac  >  4:333856/1‑134 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCAC  >  ACB122/5580308‑5580575

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