Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,404,075 A→G 5.5% N350S (AAT→AGT)  CJ019_03222 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,404,0750AG 5.5% 240.9 / 3.5 73N350S (AAT→AGT) CJ019_03222hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (2/2);  total (34/39)
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.25e-01

GACATGAGAATGATCCGACGCGCCCTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCAGGTACAGGAAGATGA  >  ACB122/3403942‑3404207
                                                                                                                                     |                                                                                                                                    
gACATGAGAATGATCCGACGCGCCCTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATg                                                                                                                                    >  7:274398/1‑136 (MQ=255)
   aTGAGAATGATCCGACGCGCCCTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCt                                                                                                                                  <  5:156995/135‑1 (MQ=255)
               cGACGCGCCCTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAggg                                                                                                                      >  1:171053/1‑135 (MQ=255)
               cGACGCGCCCTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAggg                                                                                                                      >  1:89691/1‑135 (MQ=255)
                      cccTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACtt                                                                                                              <  3:236574/136‑1 (MQ=255)
                      cccTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACtt                                                                                                              <  3:241331/136‑1 (MQ=255)
                         ttGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTg                                                                                                            >  1:364811/1‑135 (MQ=255)
                            aGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGAGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCgg                                                                                                        <  8:94565/136‑1 (MQ=255)
                                 ggCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGaa                                                                                                   <  3:303472/136‑1 (MQ=255)
                                  gCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGaa                                                                                                   <  5:211487/135‑1 (MQ=255)
                                  gCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGaa                                                                                                   <  5:109039/135‑1 (MQ=255)
                                  gCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAg                                                                                                  <  4:33367/136‑1 (MQ=255)
                                   ccGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAg                                                                                                  <  1:149088/135‑1 (MQ=255)
                                   ccGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAg                                                                                                  <  6:359482/135‑1 (MQ=255)
                                    cGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGCAGc                                                                                                 >  6:320343/1‑135 (MQ=255)
                                    cGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGc                                                                                                 >  6:122677/1‑135 (MQ=255)
                                    cGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGc                                                                                                 >  1:152215/1‑135 (MQ=255)
                                    cGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGc                                                                                                 >  1:80021/1‑135 (MQ=255)
                                     gttcccGCACCGATCAATGCGATCGAGGGTAGGGGTTGGCTGTGATCGTGATGAAGGCTGCGGTAGAACTGAGTGGGAGACTGCATGTCTACAGCTATGCTGGTGTGGTTGGGGTTGACTTTGGCGGTGGaattc                                                                                                <  1:243095/132‑4 (MQ=255)
                                      tACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCcct                                                                                              <  7:381737/136‑1 (MQ=255)
                                       aCCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCcct                                                                                              <  6:346593/135‑1 (MQ=255)
                                         ccGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCcctacc                                                                                           >  3:359479/1‑136 (MQ=255)
                                         ccGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCcctac                                                                                            >  5:299629/1‑135 (MQ=255)
                                               cGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATAcc                                                                                     >  4:41444/1‑136 (MQ=255)
                                                   cAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCTGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGc                                                                                                                                   <  8:79423/86‑1 (MQ=255)
                                                   cAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGc                                                                                                                                   >  7:79423/1‑86 (MQ=255)
                                                    aaGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAGTGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATAcc                                                                                     >  6:70697/1‑131 (MQ=255)
                                                    aaGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAGTGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATAcc                                                                                     >  6:76545/1‑131 (MQ=255)
                                                    aaGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATAccc                                                                                    >  8:146468/1‑132 (MQ=255)
                                                          aTCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGcca                                                                          <  7:322882/136‑1 (MQ=255)
                                                          aTCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGcca                                                                          <  4:322386/136‑1 (MQ=255)
                                                           tCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGcca                                                                          <  6:265400/135‑1 (MQ=255)
                                                            cGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTATGGGTTGACTTTGCGGGTGGAAGCCCTACCTATACCCAAGATGGccac                                                                         <  2:301457/135‑1 (MQ=255)
                                                              aGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGccacca                                                                       >  1:191416/1‑135 (MQ=255)
                                                              aGGAAAGTGCTTTGCTCAGAACGTGATGAAGGATGCGGAAGAACTGAGTGGGATACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCATGTGGAAGCCCTACATATACCCAAGATGGACACCaa                                                                      >  8:48468/1‑136 (MQ=255)
                                                               ggCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAAc                                                                     <  8:274398/136‑1 (MQ=255)
                                                                gCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGACCGGAGTGGGAGACTGCATGTCGACAGCAATGCTGGGGTGGGAGGGGTTGACTTTGCCGGGGGAAGCCCTACCTATACCCAAGATGGCCACCAAc                                                                     <  1:393435/135‑1 (MQ=255)
                                                                gCAAGGGCTTGGCTCAGAACGGGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAAc                                                                     <  6:191867/135‑1 (MQ=255)
                                                                    gggCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCt                                                                <  4:359479/136‑1 (MQ=255)
                                                                     ggCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCt                                                                <  5:192949/135‑1 (MQ=255)
                                                                           gCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGGTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTgg                                                         <  4:38139/136‑1 (MQ=255)
                                                                            cTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGTGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGTTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATTGCCACCAACGATCTATACTgg                                                         <  6:91604/135‑1 (MQ=255)
                                                                                   cGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAg                                                   >  5:208930/1‑134 (MQ=255)
                                                                                    gTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACtt                                                                                                              >  5:196772/1‑74 (MQ=255)
                                                                                    gTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACtt                                                                                                              <  6:196772/74‑1 (MQ=255)
                                                                                     tgatgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCAt                                               <  3:182179/136‑1 (MQ=255)
                                                                                     tgatgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCAt                                               <  7:314603/136‑1 (MQ=255)
                                                                                      gatgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCAt                                               <  1:7988/135‑1 (MQ=255)
                                                                                        tgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCAtt                                              >  3:68876/1‑134 (MQ=255)
                                                                                        tgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCAtt                                              >  3:324649/1‑134 (MQ=255)
                                                                                        tgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTCGACTTTGCCGGTGGaa                                                                                                   <  3:310765/81‑1 (MQ=255)
                                                                                        tgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTCGACTTTGCCGGTGGaa                                                                                                   >  4:310765/1‑81 (MQ=255)
                                                                                        tgaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATCCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACACAAGATGGCCACCAACGATCTATACTGGGGACAGTCAt                                               >  1:218449/1‑133 (MQ=255)
                                                                                         gaAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAg                                           >  4:226302/1‑136 (MQ=255)
                                                                                           aGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGt                                          <  2:89691/135‑1 (MQ=255)
                                                                                              cTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTc                                      >  4:149/1‑136 (MQ=255)
                                                                                                  ggAAGAACTGAGTGGGAGACTGCATGTCTACAGCAGTGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCg                                     <  5:70697/133‑1 (MQ=255)
                                                                                                  ggAAGAACTGAGTGGGAGACTGCATGTCTACAGCAGTGCTGGTGTGGTAGGGGTTGACTNTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCg                                     <  5:76545/133‑1 (MQ=255)
                                                                                                   gaagaaCTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTAc                                 <  8:218736/136‑1 (MQ=255)
                                                                                                    aagaaCTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGTTCGGTAc                                 >  1:256633/1‑135 (MQ=255)
                                                                                                     agaaCTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTAcc                                <  5:150802/135‑1 (MQ=255)
                                                                                                              tGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACTATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGAc                      >  4:131423/1‑136 (MQ=255)
                                                                                                               cgcagaCTGCATGTCTACAGCAATTATTTTGTTGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGAc                      >  8:209417/4‑135 (MQ=255)
                                                                                                                   gaCTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATAACCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGa                    >  5:136904/1‑133 (MQ=255)
                                                                                                                    aCTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCa                 >  2:236251/1‑135 (MQ=255)
                                                                                                                      tGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCAgg               <  2:171053/135‑1 (MQ=255)
                                                                                                                         aTGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCccta                                                                                             <  6:11083/54‑1 (MQ=255)
                                                                                                                         aTGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCccta                                                                                             >  5:11083/1‑54 (MQ=255)
                                                                                                                                 aGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCcctac                                                                                            >  7:80319/1‑47 (MQ=255)
                                                                                                                                 aGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCcctac                                                                                            <  8:80319/47‑1 (MQ=255)
                                                                                                                                 aGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCAGGTACAGGAAgatg   <  7:73046/136‑1 (MQ=255)
                                                                                                                                   cAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCAGGTACAGGAAgatga  >  7:59876/1‑135 (MQ=255)
                                                                                                                                   cAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCAGGTACAGGAAgatga  >  3:145791/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GACATGAGAATGATCCGACGCGCCCTTGAGCTAGGCCGTACCCGCACCGATCAAGGCCATCGAGGCAAGGGCTTGGCTCAGAACGTGATGAAGGCTGCGGAAGAACTGAGTGGGAGACTGCATGTCTACAGCAATGCTGGTGTGGTAGGGGTTGACTTTGCCGGTGGAAGCCCTACCTATACCCAAGATGGCCACCAACGATCTATACTGGGGACAGTCATTCAGTGGTCGGTACCGATCACAGACGAGCAGGTACAGGAAGATGA  >  ACB122/3403942‑3404207

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: ‑