Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,256,441 T→G 7.5% T71P (ACG→CCG)  cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,4410TG 7.5% 263.5 / 13.3 93T71P (ACG→CCG) cyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base T (39/47);  new base G (2/5);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTT  >  NC_000913/3256292‑3256589
                                                                                                                                                     |                                                                                                                                                    
gATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGt                                                                                                                                                      <  1:424245/150‑1 (MQ=255)
  tCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGc                                                                                                                                                    <  2:1551453/150‑1 (MQ=255)
   cTTAACGGAGACTTTCCCCGCCGCCAGCAGTTCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGcc                                                                                                                                                   <  2:163287/150‑1 (MQ=255)
   cTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGcc                                                                                                                                                   >  2:602298/1‑150 (MQ=255)
    ttAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCTGCAGCCCCACCATTCCCGTGCCg                                                                                                                                                  <  2:1233976/150‑1 (MQ=255)
    ttAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGGCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCg                                                                                                                                                  <  2:1307280/150‑1 (MQ=255)
     tAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCgg                                                                                                                                                 <  2:1589297/150‑1 (MQ=255)
      aaCGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCggg                                                                                                                                                <  2:1360962/150‑1 (MQ=255)
           agaCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACgg                                                                                                                                           <  1:467279/150‑1 (MQ=255)
           agaCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACATCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACgg                                                                                                                                           <  1:193362/150‑1 (MQ=255)
            gaCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGt                                                                                                                                          >  2:1263773/1‑150 (MQ=255)
             aCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTg                                                                                                                                         <  1:1617986/150‑1 (MQ=255)
              cTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTgg                                                                                                                                        <  1:1280909/150‑2 (MQ=255)
              cTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGa                                                                                                                                        <  2:1499573/150‑1 (MQ=255)
               tttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTggg                                                                                                                                       <  1:236818/150‑3 (MQ=255)
               tttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTggg                                                                                                                                       <  2:1398591/150‑3 (MQ=255)
               tttCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGGTGCGTCTTGCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCGTTCCCGGCACGGGAATGGTggg                                                                                                                                       <  1:1402971/150‑3 (MQ=255)
                 tCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGc                                                                                                                                     >  2:207972/1‑150 (MQ=255)
                   cccGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGccc                                                                                                                                   >  1:682852/1‑150 (MQ=255)
                    ccgccgCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCa                                                                                                                                  >  2:454055/1‑150 (MQ=255)
                      gccgccAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGa                                                                                                                                <  1:603522/150‑1 (MQ=255)
                       ccgccAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGGCCAGag                                                                                                                               <  2:639549/150‑2 (MQ=255)
                        cgccAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGAcc                                                                                                                              <  1:207972/150‑1 (MQ=255)
                         gccAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCg                                                                                                                             >  1:186408/1‑150 (MQ=255)
                                gTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtca                                                                                                                      <  1:965126/150‑1 (MQ=255)
                                  gCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                                                                    <  2:664033/150‑1 (MQ=255)
                                  gCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGGCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                                                                    <  2:944920/150‑1 (MQ=255)
                                    tttGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAg                                                                                                                  <  2:1260932/150‑1 (MQ=255)
                                        gCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAttt                                                                                                              >  2:1328413/1‑150 (MQ=255)
                                        gCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCAGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAttt                                                                                                              <  2:923289/150‑1 (MQ=255)
                                           tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGc                                                                                                           <  1:55832/150‑1 (MQ=255)
                                           tCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGc                                                                                                           <  1:1622511/150‑1 (MQ=255)
                                            cGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCCCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACACCCAGACCGTTCTTCATCAGATTTGGCg                                                                                                          <  2:1290398/150‑1 (MQ=255)
                                             ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                         >  2:892407/1‑150 (MQ=255)
                                             ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                         <  2:1175195/150‑1 (MQ=255)
                                             ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                         >  2:1631953/1‑150 (MQ=255)
                                             ggCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                                                                         >  2:1018568/1‑150 (MQ=255)
                                                      cTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGc                                                                                                >  1:1053948/1‑150 (MQ=255)
                                                       tGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCt                                                                                               >  1:317349/1‑150 (MQ=255)
                                                          gcTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCt                                                                                            <  2:1341070/150‑1 (MQ=255)
                                                          gcTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCt                                                                                            <  1:1289951/150‑1 (MQ=255)
                                                          gcTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCt                                                                                            <  1:220859/150‑1 (MQ=255)
                                                             gTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTaca                                                                                         >  2:62388/1‑150 (MQ=255)
                                                                gCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGt                                                                                      >  2:660308/1‑150 (MQ=255)
                                                                gCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACAGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGt                                                                                      >  1:291294/1‑150 (MQ=255)
                                                                 cGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGtt                                                                                     <  2:72589/150‑1 (MQ=255)
                                                                  gTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTc                                                                                    <  1:454055/150‑1 (MQ=255)
                                                                     tttCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCACTCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGCCGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCacc                                                                                 >  1:1129833/1‑148 (MQ=255)
                                                                      ttCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAcc                                                                                >  1:1490723/1‑150 (MQ=255)
                                                                      ttCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAcc                                                                                <  1:1335468/150‑1 (MQ=255)
                                                                          gCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAc                                                                                 >  1:486071/1‑145 (MQ=255)
                                                                          gCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAc                                                                                 <  2:486071/145‑1 (MQ=255)
                                                                          gCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAc                                                                            >  2:153031/1‑150 (MQ=255)
                                                                          gCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAc                                                                            <  1:846499/150‑1 (MQ=255)
                                                                               tCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGCCCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTcc                                                                       >  2:171304/1‑150 (MQ=255)
                                                                               tCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTcc                                                                       <  2:1264606/150‑1 (MQ=255)
                                                                               tCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTcc                                                                       <  1:243021/150‑1 (MQ=255)
                                                                                 cAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAg                                                                     <  2:1358841/150‑1 (MQ=255)
                                                                                  aGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGt                                                                    >  1:1432071/1‑150 (MQ=255)
                                                                                   gCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCATACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGtt                                                                   >  1:422409/1‑150 (MQ=255)
                                                                                   gCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGtt                                                                   >  2:1425631/1‑150 (MQ=255)
                                                                                     ccGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGc                                                                                                           >  1:1031864/1‑108 (MQ=255)
                                                                                     ccGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGc                                                                                                           <  2:1031864/108‑1 (MQ=255)
                                                                                             ggCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCgcagc                                                         <  1:1237210/150‑1 (MQ=255)
                                                                                               cATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCaa                                                       >  2:989918/1‑150 (MQ=255)
                                                                                                  ttCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccg                                                    <  1:602298/150‑1 (MQ=255)
                                                                                                     cACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccgccc                                                 >  2:1632788/1‑149 (MQ=255)
                                                                                                      aCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGCACCTTCCAGTTCTGCCGCAGCAAccgcccc                                                >  2:607820/1‑148 (MQ=255)
                                                                                                          aaCGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCgccgccgcc                                            <  2:1490723/150‑1 (MQ=255)
                                                                                                              cccccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAATTGGCGAAACCCAGGCTTCTACAGGTTAAACCGGACCTTCCAGTTCAGCCGCAGCAACCCCCGCcgccagcg                                        >  2:514119/1‑150 (MQ=255)
                                                                                                               ccccAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCcgccagcgc                                       >  1:1531887/1‑150 (MQ=255)
                                                                                                                   aGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGt                                   <  1:892407/150‑1 (MQ=255)
                                                                                                                    gcgcCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTg                                  <  2:161020/150‑1 (MQ=255)
                                                                                                                         gccgcAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATc                             >  1:815689/1‑150 (MQ=255)
                                                                                                                          acaccaTCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCg                            <  2:1144240/145‑1 (MQ=255)
                                                                                                                            gcAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGt                          >  1:893442/1‑150 (MQ=255)
                                                                                                                              aaTCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTc                        <  1:1263773/150‑1 (MQ=255)
                                                                                                                                tCGGCAGCCCCACCATTCCCGTGCCGGGAACTGTGACGCCCATACCTTTCTTCATCAGATTTGGCGAAACCGAGTCTTCTTGTCGTTCAACCGGACCTTCCAGCTCTGCCGCTTCTACCGCCGCCGCCAGCTCCAGTGAAATCGGTTCAg                      <  1:1445207/150‑1 (MQ=255)
                                                                                                                                 cGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCCTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGt                     <  1:317265/150‑1 (MQ=255)
                                                                                                                                  ggCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTa                    <  1:1425631/150‑1 (MQ=255)
                                                                                                                                  ggCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCACTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTa                    >  2:644963/1‑150 (MQ=255)
                                                                                                                                      gCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATc                >  2:1344103/1‑150 (MQ=255)
                                                                                                                                       cgtcaccgTTCCCGGCACGGGAATGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATcc               >  1:96647/4‑150 (MQ=255)
                                                                                                                                       cgtcaccgTTCCCGGCACGGGAATGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATcc               >  2:130793/4‑150 (MQ=255)
                                                                                                                                        cccACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATccc              >  2:1051171/1‑150 (MQ=255)
                                                                                                                                            ccaTTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAgcg          >  1:676084/1‑150 (MQ=255)
                                                                                                                                              aTTCCCGTGCCGGGAACGGTGACGCCCATACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGcc        <  2:422409/150‑1 (MQ=255)
                                                                                                                                                 cccGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt     >  1:601474/1‑150 (MQ=255)
                                                                                                                                                 cccGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt     >  1:465217/1‑150 (MQ=255)
                                                                                                                                                 cccGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGt     >  2:686871/1‑150 (MQ=255)
                                                                                                                                                  ccGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGtt    >  2:594370/1‑150 (MQ=255)
                                                                                                                                                  ccGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGtt    <  2:186408/150‑1 (MQ=255)
                                                                                                                                                   cGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGttt   <  1:686871/150‑1 (MQ=255)
                                                                                                                                                    gTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGtttt  <  1:1344103/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGCCAGCGCCAGTGAAATCGGTTCAGTACATCCCAGCGCCGGTTTT  >  NC_000913/3256292‑3256589

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