Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,803 T→C 6.8% K201E (AAG→GAG)  ynjF ← putative phosphatidyl transferase, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,840,8030TC 6.8% 265.8 / 7.0 88K201E (AAG→GAG) ynjFputative phosphatidyl transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/45);  new base C (2/4);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCC  >  NC_000913/1840656‑1840950
                                                                                                                                                   |                                                                                                                                                   
acgcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCa                                                                                                                                                   <  2:767935/150‑1 (MQ=255)
acgcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCa                                                                                                                                                   >  2:234353/1‑150 (MQ=255)
 cgcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAg                                                                                                                                                  >  2:71761/1‑150 (MQ=255)
   cacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCg                                                                                                                                                >  1:728889/1‑150 (MQ=255)
     cgcgCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                              <  1:709889/150‑1 (MQ=255)
     cgcgCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                              <  2:200183/150‑1 (MQ=255)
     cgcgCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTc                                                                                                                                              <  1:1227698/150‑1 (MQ=255)
     cgcgCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAACGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                              <  1:1256977/150‑1 (MQ=255)
      gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCa                                                                                                                                             <  1:196301/150‑1 (MQ=255)
      gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCa                                                                                                                                             <  2:1369518/150‑1 (MQ=255)
       cgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAg                                                                                                                                            >  2:959347/1‑150 (MQ=255)
       cgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAg                                                                                                                                            <  1:268367/150‑1 (MQ=255)
         ccATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGAt                                                                                                                                          <  2:1318538/150‑1 (MQ=255)
          cATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATa                                                                                                                                         <  1:789527/150‑1 (MQ=255)
           aTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAc                                                                                                                                        >  2:503142/1‑150 (MQ=255)
           aTGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCACCTTCAGCGTCAGATAc                                                                                                                                        >  1:886905/1‑150 (MQ=255)
            tGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc                                                                                                                                       <  1:1221608/150‑1 (MQ=255)
             gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAccc                                                                                                                                      <  2:884472/150‑1 (MQ=255)
             gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAgtaattg                                                                                                                                      <  2:1094651/150‑7 (MQ=255)
             gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAgtaattg                                                                                                                                      <  1:346549/150‑7 (MQ=255)
                   gaagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAc                                                                                                                                        <  1:1158475/142‑1 (MQ=255)
                   gaagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAc                                                                                                                                        >  2:1158475/1‑142 (MQ=255)
                   gaagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTccc                                                                                                                                >  2:1358381/1‑149 (MQ=255)
                    aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcac                                                                                                                               >  2:1269136/1‑150 (MQ=255)
                     agaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcaca                                                                                                                              >  2:1007045/1‑150 (MQ=255)
                           tttCCGTGTATGACGGCGGCTGGTACGAATCGAGCAGCGATCCAAAAAATCCGGAAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTTACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcc                                                                                                                        >  2:291853/1‑149 (MQ=255)
                            ttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGt                                                                                                                       <  2:1517101/150‑1 (MQ=255)
                             tCCGTGTATGCCGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGCCGCTGGAGCGACTTCAGCGTCAGATACCCACTCCCCCCCCGCGTa                                                                                                                      >  2:950653/1‑150 (MQ=255)
                                  gtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAGCCACTCCACACCCGCGTAAAGgt                                                                                                                 <  2:963689/150‑1 (MQ=255)
                                  gtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgt                                                                                                                 >  2:537107/1‑150 (MQ=255)
                                   tATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                                                <  2:412249/150‑1 (MQ=255)
                                    aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                                               >  2:473664/1‑150 (MQ=255)
                                     tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgt                                                                                                              <  2:923335/150‑1 (MQ=255)
                                       aCGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCa                                                                                                            <  1:537107/150‑1 (MQ=255)
                                          gcggcTGGTACGAATGGAGCAGCTATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATcc                                                                                                         <  1:329546/150‑1 (MQ=255)
                                               tGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGcac                                                                                                    <  2:262724/150‑1 (MQ=255)
                                                 gTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACaa                                                                                                  <  1:234353/150‑1 (MQ=255)
                                                   aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACg                                                                                                >  1:398615/1‑150 (MQ=255)
                                                     gAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCa                                                                                              <  1:1269136/150‑1 (MQ=255)
                                                          gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGaa                                                                                         >  2:1429111/1‑150 (MQ=255)
                                                          gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGaa                                                                                         >  2:1268705/1‑150 (MQ=255)
                                                           agcagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAg                                                                                        >  1:1207014/1‑150 (MQ=255)
                                                            gcagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGa                                                                                       <  2:1579745/150‑1 (MQ=255)
                                                              agcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATc                                                                                     <  2:1499904/150‑1 (MQ=255)
                                                               gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCCCTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGCTcc                                                                                    >  2:1243846/1‑150 (MQ=255)
                                                                 gATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccac                                                                                  >  2:609432/1‑150 (MQ=255)
                                                                 gATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccac                                                                                  >  2:932257/1‑150 (MQ=255)
                                                                   tCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgc                                                                                <  1:959347/150‑1 (MQ=255)
                                                                    ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCAcc                                                                                            >  2:538157/1‑137 (MQ=255)
                                                                    ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCAcc                                                                                            <  1:538157/137‑1 (MQ=255)
                                                                    ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgca                                                                               >  2:872805/1‑150 (MQ=255)
                                                                    ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgca                                                                               <  1:187125/150‑1 (MQ=255)
                                                                     cAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaa                                                                              >  2:1620567/1‑150 (MQ=255)
                                                                          aaTCCGGTAGCAACCGGTGAACGCTGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccac                                                                         <  1:1218640/150‑1 (MQ=255)
                                                                           aTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacg                                                                        <  2:728889/150‑1 (MQ=255)
                                                                               ggTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacgcaaa                                                                    >  2:675333/1‑150 (MQ=255)
                                                                               ggTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacgcaaa                                                                    >  1:762701/1‑150 (MQ=255)
                                                                               ggTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATACAGCACAACGCACCGAAGATccacgcaaaccacgcaaa                                                                    >  2:940054/1‑150 (MQ=255)
                                                                                     aaCCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCAAACACGAGTAAAGGGCGTCATCCAGCACAACGCACCGAATATCCACTCAAACCACGCAAACCATGc                                                              <  2:558043/150‑1 (MQ=255)
                                                                                                gcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAAc                                                   <  1:1620567/150‑1 (MQ=255)
                                                                                                     ccGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATccc                                              <  2:554336/150‑1 (MQ=255)
                                                                                                       ggACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAg                                            >  2:270869/1‑150 (MQ=255)
                                                                                                        gACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGc                                           <  1:1268705/150‑1 (MQ=255)
                                                                                                          cAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGcac                                         >  1:1231011/1‑150 (MQ=255)
                                                                                                          cAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGcac                                         <  2:142331/150‑1 (MQ=255)
                                                                                                              aGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAAc                                     >  1:1450265/1‑150 (MQ=255)
                                                                                                                cAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAg                                   >  1:1600725/1‑150 (MQ=255)
                                                                                                                  aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTa                                 <  2:398615/150‑1 (MQ=255)
                                                                                                                    tAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAg                               >  2:1338389/1‑150 (MQ=255)
                                                                                                                        aTTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCg                           <  1:121531/150‑1 (MQ=255)
                                                                                                                           gCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGttt                        >  1:289531/1‑150 (MQ=255)
                                                                                                                             aaTTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCt                      <  1:513240/150‑1 (MQ=255)
                                                                                                                              gtatcTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtg                     >  2:630421/5‑150 (MQ=255)
                                                                                                                              aTTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtg                     <  1:1007045/150‑1 (MQ=255)
                                                                                                                              aTTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtg                     <  2:386873/150‑1 (MQ=255)
                                                                                                                              aTTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtg                     <  1:473664/150‑1 (MQ=255)
                                                                                                                                tACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtgtg                   <  1:932257/150‑1 (MQ=255)
                                                                                                                                tACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtgtg                   <  1:71761/150‑1 (MQ=255)
                                                                                                                                  cTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGcc                 >  2:541329/1‑150 (MQ=255)
                                                                                                                                  cTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGcc                 >  2:635770/1‑150 (MQ=255)
                                                                                                                                     aCGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGAACAACGCACCGAAGATCCACGCAAACCAAGCAAACCATGCAGGGATTAAACATACCAGCACAAACAGTAAGATCGTTTCTGTGCCTTc              >  2:435174/1‑150 (MQ=255)
                                                                                                                                        cTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGt           <  1:675333/150‑1 (MQ=255)
                                                                                                                                         tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt          <  2:762701/150‑1 (MQ=255)
                                                                                                                                         tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt          <  1:503142/150‑1 (MQ=255)
                                                                                                                                         tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt          <  1:609432/150‑1 (MQ=255)
                                                                                                                                              cGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATcc     <  1:270869/150‑1 (MQ=255)
                                                                                                                                                 cTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATccccc  >  1:1646768/1‑150 (MQ=255)
                                                                                                                                                 cTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATccccc  >  2:1165191/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCC  >  NC_000913/1840656‑1840950

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