Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,672,481 A→T 7.6% A174A (GCA→GCT ydgD → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,4810AT 7.6% 313.6 / 8.4 114A174A (GCA→GCTydgDputative peptidase
Reads supporting (aligned to +/- strand):  ref base A (52/53);  new base T (5/4);  total (57/57)
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

Only 100 of 114 total aligned reads displayed.

GAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAG  >  NC_000913/1672333‑1672627
                                                                                                                                                    |                                                                                                                                                  
gAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCAt                                                                                                                                                   >  1:1699872/1‑150 (MQ=255)
gAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCAt                                                                                                                                                   >  2:891409/1‑150 (MQ=255)
 aaaGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCAtt                                                                                                                                                  <  1:1609513/150‑1 (MQ=255)
 aaaGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCAtt                                                                                                                                                  <  1:1213018/150‑1 (MQ=255)
   aGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTaa                                                                                                                                                >  2:1488273/1‑150 (MQ=255)
    gCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTaaa                                                                                                                                               <  2:1484043/150‑1 (MQ=255)
      ggTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAg                                                                                                                                             <  2:773430/150‑1 (MQ=255)
      ggTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAg                                                                                                                                             >  1:801339/1‑150 (MQ=255)
       gTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGc                                                                                                                                            >  1:1212027/1‑150 (MQ=255)
          aaaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGTGGc                                                                                                                                         <  1:1298079/150‑1 (MQ=255)
          aaaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGc                                                                                                                                         <  2:631170/150‑1 (MQ=255)
           aaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCa                                                                                                                                        <  1:1378704/150‑1 (MQ=255)
           aaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCa                                                                                                                                        <  1:788108/150‑1 (MQ=255)
           aaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCa                                                                                                                                        >  2:653688/1‑150 (MQ=255)
           aaaGCAGATGGGGATGGTTGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAg                                                                                                                                       >  2:1320437/1‑150 (MQ=255)
            aaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCTGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAg                                                                                                                                       <  1:1379651/150‑1 (MQ=255)
            aaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGGGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcgg                                                                                                                                       <  1:1610047/150‑3 (MQ=255)
             cgCAGATGGTGATGGTTGGATTGTACCTCGCGCAGCCGCGACGTTGGACTTCGGATTGATTGTGCTAGGTATTCCCCCTTCAGGCATTACGCCGGTGCCGTTATTTGAGGGATATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcggt                                                                                                                                      <  2:1077490/149‑4 (MQ=255)
             aGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcggt                                                                                                                                      <  2:182835/150‑4 (MQ=255)
             aGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAACGCGGCAgg                                                                                                                                      <  1:745361/150‑1 (MQ=255)
             aGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAGCGGCAGGt                                                                                                                                     >  1:1286222/1‑150 (MQ=255)
             aGCAGATGGGGATGGTGGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCAGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcggt                                                                                                                                      <  1:986770/150‑4 (MQ=255)
                aGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCg                                                                                                                                   <  2:1476279/150‑1 (MQ=255)
                     gggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGGGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAg                                                                                                                              >  2:1475038/1‑150 (MQ=255)
                       gATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTg                                                                                                                            <  2:766843/150‑1 (MQ=255)
                       gATGGTTGGATTGTACCTACCGCAGCGGCTCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGTCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTg                                                                                                                            <  2:420574/150‑1 (MQ=255)
                         tGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGAc                                                                                                                          <  1:98576/150‑1 (MQ=255)
                           gTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTc                                                                                                                        >  2:1383945/1‑150 (MQ=255)
                              ggATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcagg                                                                                                                     <  2:960932/150‑1 (MQ=255)
                               gATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcaggc                                                                                                                    >  2:700295/1‑150 (MQ=255)
                                  tGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcaggcagg                                                                                                                 >  1:1172708/1‑150 (MQ=255)
                                     aCCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTa                                                                                                              >  1:1476915/1‑150 (MQ=255)
                                      ccTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTAc                                                                                                             <  1:1685647/150‑1 (MQ=255)
                                         cccGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCt                                                                                                          >  2:240435/1‑150 (MQ=255)
                                           cgcgGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGa                                                                                                        >  1:943902/1‑150 (MQ=255)
                                           cGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTCATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGa                                                                                                        <  1:96433/150‑1 (MQ=255)
                                             cAGCCGCGCCGTGTGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAg                                                                                                      >  2:49039/1‑150 (MQ=255)
                                              aGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGa                                                                                                     <  1:1320437/150‑1 (MQ=255)
                                              aGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCATCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGa                                                                                                     >  1:979254/1‑150 (MQ=255)
                                                 cgcgCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatca                                                                                                  >  1:13187/1‑150 (MQ=255)
                                                     ccGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAtctc                                                                                              >  1:1642155/1‑150 (MQ=255)
                                                     ccGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAtctc                                                                                              <  1:653688/150‑1 (MQ=255)
                                                          ggACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATAc                                                                                         <  1:1136085/150‑1 (MQ=255)
                                                          ggACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATAc                                                                                         <  1:1488273/150‑1 (MQ=255)
                                                             cTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGtt                                                                                      <  2:119956/150‑1 (MQ=255)
                                                              ttCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTg                                                                                     <  2:282629/150‑1 (MQ=255)
                                                              ttCGGATTGATTGTGCTACATAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTg                                                                                     <  1:381305/150‑1 (MQ=255)
                                                                  gattgattGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                 >  1:1535816/1‑150 (MQ=255)
                                                                      gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                             >  1:521181/1‑150 (MQ=255)
                                                                       attgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat                                                                            >  2:1557895/1‑150 (MQ=255)
                                                                              tACGTAATCCCCCTTCTGGAATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACt                                                                     >  1:265452/1‑150 (MQ=255)
                                                                              tACGTAATCACCCTTCTGGCATTACGCCGTTGCGGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGATGATCATCTCGATACGTTGTACAGTCATCAAAACt                                                                     >  1:825091/1‑150 (MQ=255)
                                                                               aCGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtg                                                                    >  1:307941/1‑150 (MQ=255)
                                                                                 ttAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCTGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgtg                                                                  <  1:177485/149‑1 (MQ=255)
                                                                                  tAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGa                                                                 >  2:208623/1‑150 (MQ=255)
                                                                                  tAATCCCCATTATGGCATTACGCCCGTGCCGTTTTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCAGTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGGCATCAAAACTGTGa                                                                 <  1:284291/150‑1 (MQ=255)
                                                                                    aTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAg                                                               <  1:1294990/150‑1 (MQ=255)
                                                                                         ccTTCTGGCATTACGACGTTGCAGTTATTTGAGGGAGATAAAGCCGCGATTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGAct                                                          <  2:1274538/150‑1 (MQ=255)
                                                                                            tCTGGCAGTACGCCGTGGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActggc                                                       <  2:1402287/150‑1 (MQ=255)
                                                                                             cTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActggct                                                      >  2:1115683/1‑150 (MQ=255)
                                                                                                gCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTggg                                                   <  2:1172708/150‑1 (MQ=255)
                                                                                                gCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGAGCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTggg                                                   <  1:168270/150‑1 (MQ=255)
                                                                                                 cATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGgc                                                  <  2:1468947/150‑1 (MQ=255)
                                                                                                  aTTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGgcg                                                 <  1:891409/150‑1 (MQ=255)
                                                                                                        ccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACg                                           >  2:891652/1‑150 (MQ=255)
                                                                                                         cgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTc                                         <  2:1286222/150‑1 (MQ=255)
                                                                                                           ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                        >  1:750829/1‑150 (MQ=255)
                                                                                                           ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                        >  2:1104314/1‑150 (MQ=255)
                                                                                                            tgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCgg                                       <  2:1224342/150‑1 (MQ=255)
                                                                                                            tgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCgg                                       <  2:1409461/150‑1 (MQ=255)
                                                                                                               cgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAGCTCGATATGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGa                                    <  1:1070658/150‑1 (MQ=255)
                                                                                                                  tATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGTCTGGGCGCAAACGTCGGTGATGt                                 <  2:979254/150‑1 (MQ=255)
                                                                                                                   aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTc                                <  2:801339/150‑1 (MQ=255)
                                                                                                                   aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCAGCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTc                                >  1:1135132/1‑150 (MQ=255)
                                                                                                                      tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTcaca                             <  2:307941/150‑1 (MQ=255)
                                                                                                                        tgggAGATAAAGCCGCGCTTACCGCCGCATTATAAGCGGCAGGTCGTAAAGTGAATCAGGCAGGCTACCCGGAAGGTCAGCTCGATACGATGTAAAGTCAACAAAACTGTGAAGTGACTGGGTGGGAGCAAACGTCGGGGCTGTCAAATc                           >  2:71506/2‑150 (MQ=255)
                                                                                                                        aGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATc                           <  1:208623/150‑1 (MQ=255)
                                                                                                                        aGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATc                           >  1:470766/1‑150 (MQ=255)
                                                                                                                          ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAg                         >  2:397416/1‑150 (MQ=255)
                                                                                                                          ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAg                         >  2:566199/1‑150 (MQ=255)
                                                                                                                          ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCG‑‑GATGTCACATCAGTg                       >  2:1094513/1‑150 (MQ=255)
                                                                                                                           gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGt                        <  2:659864/150‑1 (MQ=255)
                                                                                                                             gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGc                      <  1:1557895/150‑1 (MQ=255)
                                                                                                                              aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCg                     <  2:459309/150‑1 (MQ=255)
                                                                                                                                        gcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGcc           >  1:1571155/1‑150 (MQ=255)
                                                                                                                                         cTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCg          >  2:376848/1‑150 (MQ=255)
                                                                                                                                          ttACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCgg         <  2:943902/150‑1 (MQ=255)
                                                                                                                                          gacctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCgg         >  1:1407449/6‑150 (MQ=255)
                                                                                                                                           tACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCggg        >  1:752819/1‑150 (MQ=255)
                                                                                                                                           tACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCggg        >  2:9253/1‑150 (MQ=255)
                                                                                                                                           acctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGAACAGTCATCAAAACAGTGAAGTGACTGGCAGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCATGCCggg        >  2:503561/5‑150 (MQ=255)
                                                                                                                                            aCCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGt       <  2:1476915/150‑1 (MQ=255)
                                                                                                                                            aCCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGt       >  2:440162/1‑150 (MQ=255)
                                                                                                                                             ccgccgCATTAAAAGCGGCAGGTCGTAAAGTGACTGAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGCTGATGTCACATCAGTGCGATACCTTGCCGGGTg      <  2:1022458/150‑1 (MQ=255)
                                                                                                                                             ccgccgCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTg      <  1:534850/150‑1 (MQ=255)
                                                                                                                                              cgccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGa     >  1:527754/1‑150 (MQ=255)
                                                                                                                                               gccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGAc    <  2:1290489/150‑1 (MQ=255)
                                                                                                                                                 cgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGCGATGTCACATCAGTGCGATACCTTGCCGGGTGACAg  >  2:228578/1‑150 (MQ=255)
                                                                                                                                                  gcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTcaca                             <  2:121139/122‑1 (MQ=255)
                                                                                                                                                  gcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTcaca                             >  1:121139/1‑122 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAG  >  NC_000913/1672333‑1672627

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