Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,945 A→G 6.1% M185T (ATG→ACG)  PP_3571 ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,051,9450AG 6.1% 192.0 / 2.4 66M185T (ATG→ACG) PP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base G (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

GAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAAC  >  NC_002947/4051800‑4052091
                                                                                                                                                 |                                                                                                                                                  
gAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGc                                                                                                                                                <  1:729437/150‑1 (MQ=255)
 aaGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGATGCGCCCTGCACCTGGCAACGCATTGcc                                                                                                                                               >  2:260739/1‑150 (MQ=255)
 aaGCGCTCGTAGTTATGATCCGGCACTGCAACGGCCCCGCCGCTTTCACCCAACCTCCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCTTCCGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGcc                                                                                                                                               <  2:117680/150‑1 (MQ=255)
            gtGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCtt                                                                                                                                    <  1:594560/150‑1 (MQ=255)
               atgaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTtgtg                                                                                                                                 >  1:917745/1‑150 (MQ=255)
                 gaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAgg                                                                                                                                            >  1:759939/1‑136 (MQ=255)
                 gaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAgg                                                                                                                                            <  2:759939/137‑2 (MQ=255)
                  aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAg                                                                                                                                             >  2:456405/1‑135 (MQ=255)
                  aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAg                                                                                                                                             <  1:456405/135‑1 (MQ=255)
                   aCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGcccc                                                                                                                             <  2:744779/150‑1 (MQ=255)
                      ggCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGt                                                                                                                          >  1:678719/1‑150 (MQ=255)
                              cAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGttt                                                                                                                  <  1:462685/150‑1 (MQ=255)
                                gggCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGc                                                                                                                                <  2:256474/134‑1 (MQ=255)
                                gggCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGc                                                                                                                                >  1:256474/1‑134 (MQ=255)
                                gggCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGa                                                                                                                >  2:1070553/1‑150 (MQ=255)
                                  gCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGt                                                                                                              >  1:738811/1‑150 (MQ=255)
                                       gCCGATTTCACCCAACTTGCCCATCAGTTTCGCCGTATACGCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGt                                                                                                         <  2:917745/150‑1 (MQ=255)
                                       gCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                              >  1:1195942/1‑129 (MQ=255)
                                       gCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                              <  2:1195942/129‑1 (MQ=255)
                                       gCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGt                                                                                                         <  1:497484/150‑1 (MQ=255)
                                          gATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCACATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTggg                                                                                                      <  2:738811/150‑1 (MQ=255)
                                              tCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCa                                                                                                  <  2:862946/150‑1 (MQ=255)
                                               cACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                                 <  1:692333/150‑1 (MQ=255)
                                               cACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCACATACGGAGTGGACGCCTGAGGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                                 <  1:445245/150‑1 (MQ=255)
                                                  ccAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGa                                                                                              >  2:408680/1‑150 (MQ=255)
                                                    aaCTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATc                                                                                            <  1:408680/150‑1 (MQ=255)
                                                     aCTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                           <  1:675606/150‑1 (MQ=255)
                                                          ccAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGcac                                                                                      >  2:289570/1‑150 (MQ=255)
                                                           cAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGcaca                                                                                     >  1:1105457/1‑150 (MQ=255)
                                                              tCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGcc                                                                                  <  1:289570/150‑1 (MQ=255)
                                                              gcAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGcc                                                                                  >  2:924429/2‑150 (MQ=255)
                                                               cAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCg                                                                                 <  2:1100421/150‑1 (MQ=255)
                                                                  tttCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACATGGCAACGCATTGCCAGCTTGCCCTTGTACCCCTGTAACGGGTTGAGTTGGGTGTGGTGTCGGGTCAGGGACAGCCGTgg                                                                              >  2:62581/1‑150 (MQ=255)
                                                                   ttCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGAAACGCATTGCCAGCTTGGCCTTGTGCCGCAGTACCGGGTTGAGTTCGGTGGGTGCACCGATCAGGCACAGCCGTGGt                                                                             >  1:498966/1‑150 (MQ=255)
                                                                      tcttgttACGCGATGGCGTTCGCCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTa                                                                          <  2:498966/144‑1 (MQ=255)
                                                                       ccgcATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                           <  1:857070/132‑1 (MQ=255)
                                                                       ccgcATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                           >  2:857070/1‑132 (MQ=255)
                                                                            tACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACggg                                                                    <  1:260739/150‑1 (MQ=255)
                                                                             aCTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGgc                                                                   <  1:1070553/150‑1 (MQ=255)
                                                                                  aTGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGt                                                              <  2:941939/150‑1 (MQ=255)
                                                                                       gTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAg                                                         >  1:778043/1‑150 (MQ=255)
                                                                                             cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagc                                                   >  1:810377/1‑150 (MQ=255)
                                                                                             cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagc                                                   >  1:1106864/1‑150 (MQ=255)
                                                                                             cccAGACGGCGCGGACGCCGAATGGCTCGCTGCGCCCTGCACCTGGCAACGTATTCGCCGCTGGCCCTTGGGCGCCAGGACCGGTTTTAGTTCGTTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagc                                                   <  2:507810/150‑1 (MQ=255)
                                                                                              ccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCa                                                                                                  >  2:1030364/1‑102 (MQ=255)
                                                                                              ccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCa                                                                                                  <  1:1030364/102‑1 (MQ=255)
                                                                                               cATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagcca                                                 >  1:1108267/1‑150 (MQ=255)
                                                                                               cATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagcca                                                 >  2:880354/1‑150 (MQ=255)
                                                                                                aTACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagccag                                                >  1:127891/1‑150 (MQ=255)
                                                                                                               cGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGcc                                 >  1:48912/1‑150 (MQ=255)
                                                                                                                   tGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCaaa                             <  2:1106864/150‑1 (MQ=255)
                                                                                                                     gCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTg                                                                                                        <  2:1149276/73‑1 (MQ=255)
                                                                                                                     gCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTg                                                                                                        >  1:1149276/1‑73 (MQ=255)
                                                                                                                      cACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAaacaa                          <  2:1076722/150‑1 (MQ=255)
                                                                                                                          cTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGa                      <  2:778043/150‑1 (MQ=255)
                                                                                                                              gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGtt                                                                                                    >  2:1199312/1‑68 (MQ=255)
                                                                                                                              gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGtt                                                                                                    <  1:1199312/68‑1 (MQ=255)
                                                                                                                              gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCa                                     >  2:945624/1‑131 (MQ=255)
                                                                                                                              gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCa                                     <  1:945624/131‑1 (MQ=255)
                                                                                                                                ccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGa                                                          >  1:890208/1‑108 (MQ=255)
                                                                                                                                ccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGa                                                          <  2:890208/108‑1 (MQ=255)
                                                                                                                                  tGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCg              <  1:43286/150‑1 (MQ=255)
                                                                                                                                   gCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGt             <  2:127891/150‑1 (MQ=255)
                                                                                                                                            aaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGaaa    >  2:236109/1‑150 (MQ=255)
                                                                                                                                             aCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGaaaa   >  2:533128/1‑150 (MQ=255)
                                                                                                                                              cGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAAc  <  2:1105457/150‑1 (MQ=255)
                                                                                                                                                 aTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa           >  2:1177388/1‑138 (MQ=255)
                                                                                                                                                 aTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCAGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa           <  1:1177388/138‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAAC  >  NC_002947/4051800‑4052091

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