Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,291,109 A→T 5.1% K766* (AAA→TAA) ‡ yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,1090AT 5.1% 466.3 / 11.2 131K766* (AAA→TAA) ‡yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (57/67);  new base T (5/2);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

CTCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGT  >  NC_000913/3290961‑3291258
                                                                                                                                                    |                                                                                                                                                     
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      <  2:1307283/150‑1 (MQ=255)
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      <  2:15299/150‑1 (MQ=255)
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      <  1:1531088/150‑1 (MQ=255)
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      <  2:1047144/150‑1 (MQ=255)
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      <  1:244725/150‑1 (MQ=255)
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      <  2:794016/150‑1 (MQ=255)
  cGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAAc                                                                                                                                                    <  1:453403/150‑1 (MQ=255)
  cGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAAc                                                                                                                                                    <  1:1579565/150‑1 (MQ=255)
  cGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAAc                                                                                                                                                    >  2:2174667/1‑150 (MQ=255)
      aGATACTACGCAGCTGCCCGATAACTTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGtt                                                                                                                                                <  2:1527655/150‑1 (MQ=255)
       gATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTg                                                                                                                                               >  1:2331812/1‑150 (MQ=255)
       gATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTg                                                                                                                                               <  2:1091362/150‑1 (MQ=255)
       gATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTg                                                                                                                                               >  2:249016/1‑150 (MQ=255)
       gATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTg                                                                                                                                               <  2:1091432/150‑1 (MQ=255)
        aTACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCTACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGc                                                                                                                                              <  1:2238868/150‑1 (MQ=255)
        aTACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGc                                                                                                                                              <  2:211220/150‑1 (MQ=255)
        aTACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGc                                                                                                                                              >  2:527420/1‑150 (MQ=255)
         tactacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGcc                                                                                                                                             <  2:1616110/150‑1 (MQ=255)
          actacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGccc                                                                                                                                            >  2:714719/1‑150 (MQ=255)
          actacGAAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGAGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAAGATCGGTAATCGCGTACTAGTTACAGTCAGCGATCGCAACGGTAAACCGATGccc                                                                                                                                            >  2:1914656/1‑150 (MQ=255)
            tacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCtt                                                                                                                                          <  1:1748440/150‑1 (MQ=255)
            tacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCtt                                                                                                                                          >  2:1753244/1‑150 (MQ=255)
            tacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCtt                                                                                                                                          >  2:2368194/1‑150 (MQ=255)
             acGCAGCTGCCCGCGAACGTCGATCTTGAACCCACAACACAGTTTGTGGTGCCCAACAGCGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCttt                                                                                                                                         <  2:2158664/150‑1 (MQ=255)
              cGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTg                                                                                                                                        <  2:819630/150‑1 (MQ=255)
              cGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTg                                                                                                                                        <  1:254230/150‑1 (MQ=255)
                     gCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGctct                                                                                                                                 <  1:1148521/150‑1 (MQ=255)
                      cccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGTGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                                >  2:1184167/1‑150 (MQ=255)
                      cccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                                >  1:1870413/1‑150 (MQ=255)
                      ccagATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                                <  1:1235781/147‑1 (MQ=255)
                       ccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTg                                                                                                                               >  2:1943364/1‑150 (MQ=255)
                       ccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTg                                                                                                                               <  2:796484/150‑1 (MQ=255)
                        cGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                              <  2:747878/150‑1 (MQ=255)
                        cGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                              <  2:345660/150‑1 (MQ=255)
                              cGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                        <  2:862220/150‑1 (MQ=255)
                              cGTCGATCCTGAACAAACAACACAGTGTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCCACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGCTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                        <  2:1164260/150‑1 (MQ=255)
                               gTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAAc                                                                                                                       >  2:1030899/1‑150 (MQ=255)
                                    tCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatga                                                                                                                  >  2:815503/1‑150 (MQ=255)
                                        gAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACg                                                                                                              >  1:2297291/1‑150 (MQ=255)
                                         aacaaacaACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                             <  2:1241961/150‑1 (MQ=255)
                                         aacaaacaACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                             >  2:2132914/1‑150 (MQ=255)
                                         aacaaacaACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                             <  2:1241900/150‑1 (MQ=255)
                                         aacaaacaACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                             >  1:1911230/1‑150 (MQ=255)
                                          acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACggg                                                                                                            >  1:389774/1‑150 (MQ=255)
                                          acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACggg                                                                                                            <  2:389211/150‑1 (MQ=255)
                                           caaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                           <  1:1730554/150‑1 (MQ=255)
                                                  cacaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                    <  2:1353375/150‑1 (MQ=255)
                                                       tttGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtc                                                                                               >  2:1650763/1‑150 (MQ=255)
                                                       tttGTGGTGCCCAACAGAGGTGAAATGGTAGCGGCGCGTTTCAACGCAAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtc                                                                                               >  2:1298767/1‑150 (MQ=255)
                                                            ggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTTTCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                          <  2:1856873/150‑1 (MQ=255)
                                                            ggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                          <  2:2205267/150‑1 (MQ=255)
                                                            ggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                          <  1:778612/150‑1 (MQ=255)
                                                            ggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                          <  2:1856877/150‑1 (MQ=255)
                                                            ggtgCCCAACAGAGGTGCAATGGGAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                          <  2:2205306/150‑1 (MQ=255)
                                                                cccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAg                                                                                      >  2:1692757/1‑150 (MQ=255)
                                                                   aaCAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACggg                                                                                                            <  2:1212541/125‑1 (MQ=255)
                                                                   aaCAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggc                                                                                   >  1:1773222/1‑150 (MQ=255)
                                                                      agagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggcggc                                                                                <  1:2212199/150‑1 (MQ=255)
                                                                       gagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCa                                                                               <  1:31237/150‑1 (MQ=255)
                                                                       gagGTGCAATGGTAGCGGCGCGTTACAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAATTATCGTCGATGAGGGCGGCa                                                                               <  2:1830021/150‑1 (MQ=255)
                                                                           tGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACt                                                                           <  2:536992/150‑1 (MQ=255)
                                                                             cAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtat                                                                         >  1:17214/1‑150 (MQ=255)
                                                                                        gcgcGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTggg                                                              <  1:1091918/150‑1 (MQ=255)
                                                                                        gcgcGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTggg                                                              >  2:1390182/1‑150 (MQ=255)
                                                                                        gcgcGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTggg                                                              <  1:1054278/150‑1 (MQ=255)
                                                                                         cgcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGa                                                             <  2:1462144/150‑1 (MQ=255)
                                                                                          gcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGat                                                            >  1:1134123/1‑150 (MQ=255)
                                                                                          gcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGAGGCATACTATATCTCTCTGGGat                                                            >  1:1695811/1‑150 (MQ=255)
                                                                                              ttCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCg                                                        >  1:788466/1‑150 (MQ=255)
                                                                                               tCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGa                                                       >  2:362356/1‑150 (MQ=255)
                                                                                                 aaCGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGt                                                     >  2:1247308/1‑150 (MQ=255)
                                                                                                     ccAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAAt                                                 <  1:1753244/150‑1 (MQ=255)
                                                                                                     ccAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAAt                                                 <  1:2368194/150‑1 (MQ=255)
                                                                                                          atCGGTTATCGCGTACTTGTTACAGTCAGCGTTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAGCGTCGATGAGGGCGGCATACTATATCTCTCTGGGAGCTCGAGTAAATCACaa                                            >  2:386936/1‑150 (MQ=255)
                                                                                                                  tCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGAc                                    <  1:740986/150‑1 (MQ=255)
                                                                                                                  tCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGAc                                    <  1:740980/150‑1 (MQ=255)
                                                                                                                  tCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGAc                                    <  1:282321/150‑1 (MQ=255)
                                                                                                                    gcgTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTCTATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTg                                  >  1:865174/1‑150 (MQ=255)
                                                                                                                             ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                         >  2:768095/1‑150 (MQ=255)
                                                                                                                             ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                         <  2:1233510/150‑1 (MQ=255)
                                                                                                                             ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                         <  1:2243595/150‑1 (MQ=255)
                                                                                                                             ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                         >  1:2234150/1‑150 (MQ=255)
                                                                                                                                 aGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGaaa                     >  2:1373016/1‑150 (MQ=255)
                                                                                                                                  gTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAAt                    >  2:2081680/1‑150 (MQ=255)
                                                                                                                                   tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATc                   >  2:1615638/1‑150 (MQ=255)
                                                                                                                                        gATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCa              >  1:1293772/1‑150 (MQ=255)
                                                                                                                                        gATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCa              >  2:1502989/1‑150 (MQ=255)
                                                                                                                                        gATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCa              <  1:1885658/150‑1 (MQ=255)
                                                                                                                                         aTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAg             >  1:1137213/1‑150 (MQ=255)
                                                                                                                                           cGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGAt           <  1:1390182/150‑1 (MQ=255)
                                                                                                                                           cGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGAGGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGTAATCAGGCAGAt           <  1:2212323/150‑1 (MQ=255)
                                                                                                                                            gCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATc          >  2:2265980/1‑150 (MQ=255)
                                                                                                                                            gCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATAGCTCTCTGGGCTCTCGAGTCAATCACCCAGCTGGACTGTACGCGGGGGACCTCAGGCAgcgc          >  2:2336898/1‑147 (MQ=255)
                                                                                                                                             cAACGGTAAACCGTTGCCCTTTGGCGCTCTTTCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATca         >  2:304247/1‑150 (MQ=255)
                                                                                                                                               aCGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaac       >  1:1969622/1‑150 (MQ=255)
                                                                                                                                                  gTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAAt    <  2:94508/150‑1 (MQ=255)
                                                                                                                                                   tAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATg   >  2:568779/1‑150 (MQ=255)
                                                                                                                                                   tAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATg   >  1:425393/1‑150 (MQ=255)
                                                                                                                                                   tAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATg   >  2:2065112/1‑150 (MQ=255)
                                                                                                                                                    aaaCCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGt  <  2:1132865/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGT  >  NC_000913/3290961‑3291258

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