Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,054,977 T→G 5.6% intergenic (+18/‑84) shiA → / → amn shikimate transporter/AMP nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,054,9770TG 5.6% 249.1 / 3.7 83intergenic (+18/‑84)shiA/amnshikimate transporter/AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base G (2/3);  total (40/43)
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.33e-01

CACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGA  >  NC_000913/2054829‑2055125
                                                                                                                                                    |                                                                                                                                                    
cacCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatcatc                                                                                                                                                     <  1:429892/150‑1 (MQ=255)
cacCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatcatc                                                                                                                                                     <  2:738870/150‑1 (MQ=255)
 acCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATcc                                                                                                                                                    <  1:938819/150‑1 (MQ=255)
     tttATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGGCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATc                                                                                                                                                <  1:380865/150‑1 (MQ=255)
          tGCCGCTGCTCTCATGACTTACTTTTCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGc                                                                                                                                           <  1:1004250/150‑1 (MQ=255)
          tGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTAGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGc                                                                                                                                           <  1:720159/150‑1 (MQ=255)
            ccGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCt                                                                                                                                         >  1:171549/1‑150 (MQ=255)
            ccGCTGAACTCATCACTTACTTTGACGGGAACTGGCATAGCGTCGCCAATTATTTGCTTGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACACTCAACGCGCTTGATAGCCTCGCTAAGAACAACCGATCTTCGCCt                                                                                                                                         >  2:671518/1‑150 (MQ=255)
              gCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTa                                                                                                                                       >  1:1155722/1‑150 (MQ=255)
               cTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTac                                                                                                                                      <  2:590222/150‑1 (MQ=255)
                 gCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCTCCTTacac                                                                                                                                    >  2:366656/1‑150 (MQ=255)
                  cACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGccag                                                                                                                                        >  2:150580/1‑143 (MQ=255)
                  cACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCGGATGATCTTCGccag                                                                                                                                        <  1:150580/145‑3 (MQ=255)
                  cACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACt                                                                                                                                   >  1:191610/1‑150 (MQ=255)
                       atcaCTTACTTTGCCGGGAACTGGCATAGCGTCGTCATTTATTTGCTGGCGGGAAGCCTGACTTCCGCAATGACCGCTTTGTTGATGACAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGTTCTTCGCCTTACACTTTTGt                                                                                                                              <  1:205510/150‑1 (MQ=255)
                        tcaCTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGtt                                                                                                                             >  1:1312709/1‑150 (MQ=255)
                          acttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCTATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTAGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTc                                                                                                                           <  1:669358/150‑1 (MQ=255)
                          acttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGttta                                                                                                                           >  2:379508/1‑149 (MQ=255)
                          acttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTc                                                                                                                           <  2:1128635/150‑1 (MQ=255)
                           cttacttTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTca                                                                                                                          >  1:270163/1‑150 (MQ=255)
                            ttacttTGCCGTGATCTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTTGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTcac                                                                                                                         <  1:1117674/150‑1 (MQ=255)
                                    ccGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgt                                                                                                                 >  2:1307893/1‑150 (MQ=255)
                                     cGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCtgtg                                                                                                                <  1:494309/150‑1 (MQ=255)
                                       ggAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGAc                                                                                                              <  1:1295599/150‑1 (MQ=255)
                                        gAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACa                                                                                                             <  1:220096/150‑1 (MQ=255)
                                          aCTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATa                                                                                                           <  2:624632/150‑1 (MQ=255)
                                           cTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATAc                                                                                                          >  1:1084347/1‑150 (MQ=255)
                                            tGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACt                                                                                                         <  1:475498/150‑1 (MQ=255)
                                               cATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATc                                                                                                      <  2:1053561/150‑1 (MQ=255)
                                                aTAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCg                                                                                                     <  2:850222/150‑1 (MQ=255)
                                                    cGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAtg                                                                                                 <  1:613392/150‑1 (MQ=255)
                                                          cATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGt                                                                                           >  1:277181/1‑150 (MQ=255)
                                                            tttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTaa                                                                                         >  1:402882/1‑150 (MQ=255)
                                                             ttatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAAt                                                                                        >  1:1352762/1‑150 (MQ=255)
                                                              tatttGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCATCGCGCTTGATAGCCTGGCGAAGATCATCCGCTCTTCGCCTTACACTTTTGTTTAACATTTCTGTGACATACTATCGGATGTGCGGTAAtt                                                                                       >  1:191684/1‑150 (MQ=255)
                                                                   gctggctgGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTAtt                                                                                  >  2:741789/1‑149 (MQ=255)
                                                                    ctggctggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTTTATgg                                                                                 >  1:568080/1‑150 (MQ=255)
                                                                    ctggctggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATgg                                                                                 <  2:862224/150‑1 (MQ=255)
                                                                          ggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATgcgccagg                                                                           <  1:1315075/150‑8 (MQ=255)
                                                                          ggATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAgg                                                                           <  2:1312709/150‑1 (MQ=255)
                                                                           gATGCGTGATGTCCGCAATGACCGCTTTGTTGATGAATGCCAGTCATCGAGATTGATAGCGTGTCGAAGATCATCCGTTCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGATGTAATTGTATGGAACAGga                                                                          <  1:68569/150‑1 (MQ=255)
                                                                           gATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGga                                                                          <  1:981840/150‑1 (MQ=255)
                                                                           gATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGga                                                                          >  2:1266419/1‑150 (MQ=255)
                                                                            aTGCCTGATTTCCGCAATGACCGCTTTGTTGACGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACTTTTCTGTGACATACTATCGGTTGTGCGGTAATTGTATGGAACAGgag                                                                         <  2:1162703/150‑1 (MQ=255)
                                                                               ccTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGaca                                                                      >  2:988821/1‑150 (MQ=255)
                                                                                cTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGacac                                                                     >  1:758432/1‑150 (MQ=255)
                                                                                  gATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCACCGCGCTTGATAGCCTGGCGAAGATCAGCCGATCTTCGCCTTACACTGTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTTTATGTAACAGGAGacacac                                                                   >  1:638901/1‑150 (MQ=255)
                                                                                     ttCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAAAAGGAGACACACATg                                                                >  2:131638/1‑150 (MQ=255)
                                                                                         gCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaata                                                            >  1:423597/1‑150 (MQ=255)
                                                                                         gCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaata                                                            <  1:475630/150‑1 (MQ=255)
                                                                                         gCAATGACCGCTTTGTTGATGAAAGACAGTCAAAGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGaata                                                            <  1:220139/150‑1 (MQ=255)
                                                                                             tGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAAtaataa                                                        <  1:321722/150‑1 (MQ=255)
                                                                                             tGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAAtaataa                                                        >  1:170542/1‑150 (MQ=255)
                                                                                              gACCGCTTTGTTGATGAAAGACAGTCACCGCGCTTGATAGCCTGGCGAAGATCATCCGATCATCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTTCGGAAATTGTAAGGATCAGGAGACCAACATGAATAATAAg                                                       >  1:202602/1‑150 (MQ=255)
                                                                                              gACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAg                                                       <  2:1106136/150‑1 (MQ=255)
                                                                                              gACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAg                                                       >  2:253830/1‑150 (MQ=255)
                                                                                                acGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCTGTAATTGTATGGAACAGGAGACACACATGAATAATAAggt                                                     <  1:401683/149‑2 (MQ=255)
                                                                                                  gCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCt                                                   >  1:1219897/1‑150 (MQ=255)
                                                                                                    tttGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTcc                                                 <  2:478371/150‑1 (MQ=255)
                                                                                                         tgatgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATGGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCt                                            <  1:402301/150‑1 (MQ=255)
                                                                                                           atgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTGTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGa                                          <  2:191610/150‑1 (MQ=255)
                                                                                                            tgaAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAc                                         <  2:431569/150‑1 (MQ=255)
                                                                                                             taaaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcc                                        <  2:518710/149‑1 (MQ=255)
                                                                                                              aaaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAccc                                       >  1:1273708/1‑150 (MQ=255)
                                                                                                               aaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcccc                                      >  2:829884/1‑150 (MQ=255)
                                                                                                               aaGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTGTTGTTTCACATTTCTGGGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGAcccc                                      >  2:218863/1‑150 (MQ=255)
                                                                                                                  aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATGATAAGGGCTCCGGTCTGCCCCCAGc                                   >  2:1263923/1‑150 (MQ=255)
                                                                                                                  aCAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGc                                   <  1:988821/150‑1 (MQ=255)
                                                                                                                   cAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCt                                  <  2:1155722/150‑1 (MQ=255)
                                                                                                                       cAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAgg                              <  2:756986/150‑1 (MQ=255)
                                                                                                                          cgcgCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCAc                           >  2:1143833/1‑150 (MQ=255)
                                                                                                                             gcTTGATAGCCTGGGTAAGATCATCCTATCTTCGACTTACACATTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTgg                        <  2:817301/150‑1 (MQ=255)
                                                                                                                             gcTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTgg                        <  2:397459/150‑1 (MQ=255)
                                                                                                                             gcTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTgg                        >  1:486584/1‑150 (MQ=255)
                                                                                                                              cTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACAGGa                       >  2:361084/1‑150 (MQ=255)
                                                                                                                                tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                     >  2:800575/1‑150 (MQ=255)
                                                                                                                                tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATa                     >  2:1271627/1‑150 (MQ=255)
                                                                                                                                tGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACAGGATa                     >  2:1078211/1‑150 (MQ=255)
                                                                                                                                    aGCCTGGCGAAGATCATCCGATCTTAGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTTATTGTATGGAACAGGAGACACACATTAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACt                 <  1:671518/150‑1 (MQ=255)
                                                                                                                                     gCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGCCCCCAGCGCAGGCACTGGATTAACTc                >  2:665133/1‑150 (MQ=255)
                                                                                                                                     gCCTGGCGAAGAACATCCGATCTTCGCCTTACACTTTTGTTTAAAATTTCTGTTACATAATATCGGTTGTGCGGTAATTGTATGGAACAGGAGACACCCATGAAAAATAATGGCTCCGCTCTGACCCCAGCTCAGGCAATTGATTAACta                >  2:732861/1‑149 (MQ=255)
                                                                                                                                        tGGCTAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCTGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAc             <  1:1270440/150‑1 (MQ=255)
                                                                                                                                         ggCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCTAcg            >  1:546709/1‑150 (MQ=255)
                                                                                                                                           cGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGAcgcg          <  2:270163/150‑1 (MQ=255)
                                                                                                                                              aGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTg       <  1:253830/150‑1 (MQ=255)
                                                                                                                                                   atcCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGa  >  2:1322059/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACCTTTTATTGCCGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACACACATGAATAATAAGGGCTCCGGTCTGACCCCAGCTCAGGCACTGGATAAACTCGACGCGCTGTATGA  >  NC_000913/2054829‑2055125

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: