Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,972 C→G 7.0% intergenic (+21/+10) ybaY → / ← atl PF09619 family lipoprotein YbaY/DNA base‑flipping protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,9720CG 7.0% 182.4 / 6.1 65intergenic (+21/+10)ybaY/atlPF09619 family lipoprotein YbaY/DNA base‑flipping protein
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base G (2/3);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

GCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAC  >  NC_000913/475831‑476113
                                                                                                                                             |                                                                                                                                              
gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTGCACCAACTCAGGTGAATCCGTCGTCGACAGTTACAACTCCTACGCAACAATCAGCCAGACTACCCCTCGCCGccc                                                                                                                                               >  2:620355/1‑143 (MQ=255)
gCCGGGACAGCAAACCGCCGTGCCGGTTCAGGACAGTGGCCGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGccc                                                                                                                                               <  1:462030/143‑1 (MQ=255)
 ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCg                                                                                                                                              >  1:655245/1‑143 (MQ=255)
  cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCCGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCgg                                                                                                                                             <  2:637892/143‑1 (MQ=255)
     tACAGCACTCCGCCCTCCCGGTCCAGGCCAGCGGGGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTGCGGCAGTTCCAGCTTCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGgaga                                                                                                                                          <  1:77805/143‑1 (MQ=255)
       cAGCGAACCGCCGCGCCGGTTCAGGCCAGCGGTGGCGCAACAACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCTGCAGTTCCAGCTCCTACGGAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                        <  2:173942/143‑1 (MQ=255)
        aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                       <  1:366532/143‑1 (MQ=255)
        aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                       <  1:623800/143‑1 (MQ=255)
         gCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGCCGTTAACCCTCTCCGCCCGGCGAGGGt                                                                                                                                      >  1:349351/1‑143 (MQ=255)
          cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCGGGCGAATGCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCGCCCCGCCCGCAGAgggcg                                                                                                                                     >  2:918565/1‑141 (MQ=255)
          cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGtt                                                                                                                                     <  1:615118/143‑1 (MQ=255)
                     gCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAg                                                                                                                          <  2:923918/143‑1 (MQ=255)
                        ggTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCCGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGa                                                                                                                       >  1:1048293/1‑143 (MQ=255)
                           tCAGGGCCGCCGTGGCGCAACGACTACCGGACCTTCGACCTCGCCAACCCAGGTGAATCCGTCCTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATaa                                                                                                                    <  2:911929/143‑1 (MQ=255)
                               gCCAGCGTTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGCTGAACAATGCACCATTATGGTCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCt                                                                                                                <  1:512101/143‑1 (MQ=255)
                               gCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCt                                                                                                                >  1:939114/1‑143 (MQ=255)
                               gCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACGCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTGAACCCTCTCCGCGCGGAGTGGGTTCGTAGTTCCATCGATAACGCt                                                                                                                >  1:10585/1‑143 (MQ=255)
                                ccAGCGGTGGCGCAACGACTACCGTAGCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATGTTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTCCCAGCGAGAACGCTg                                                                                                               >  2:1052869/1‑143 (MQ=255)
                                 cAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                              >  1:530264/1‑143 (MQ=255)
                                   gcgtgggCCCAACGACTACCGTACCTTAGACTGCTCCAACTCAGGTGAATCCGTCTTCAGCAGTTCCAGCTCCTACGCAAGATTAAGCGACGGCAACCCTCGCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCaa                                                                                                            <  2:486120/138‑1 (MQ=255)
                                           gcAGCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAAGATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                    <  1:335684/143‑1 (MQ=255)
                                            cAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCGATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCGAGCGATAACGCTGCACGTCggcttg                                                                                                   >  2:802309/1‑138 (MQ=255)
                                             aaCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCACTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGCTAACGCTGCAAGTCGATTTGc                                                                                                  >  2:145653/1‑143 (MQ=255)
                                                  cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCt                                                                                             >  2:940138/1‑143 (MQ=255)
                                                   tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCtt                                                                                            >  2:489141/1‑143 (MQ=255)
                                                   tACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAACCTTAACCGACGTTAACCCTCTCCGCCCGGAGCGGGTTAGTAGTTCCAGCGATCACGCTGCAAGTCGATTTGCCCGCtt                                                                                            >  1:1080729/1‑143 (MQ=255)
                                                      cGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGAAGGTAACCCTCTCCGCCCGCAGAGGGTTAATAGGTCCAGCGATAACGCTGCAAGTCGATTTGCGCGCGTccc                                                                                         >  2:99776/1‑143 (MQ=255)
                                                            ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                   >  2:981893/1‑143 (MQ=255)
                                                            ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGCTAGTAGCTCCAGCGCTGACGCTGCAAGGCGCTTTGCGCGCTTccacgtacc                                                                                   >  2:904939/1‑136 (MQ=255)
                                                                 cTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGCTTTGCCCGCTTCCCGCTACCATcac                                                                              >  1:220414/1‑143 (MQ=255)
                                                                 cTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGGTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcac                                                                              <  1:923492/143‑1 (MQ=255)
                                                                  ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTc                                                                                           >  2:921394/1‑129 (MQ=255)
                                                                  ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTc                                                                                           <  1:921394/129‑1 (MQ=255)
                                                                     aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACCCCt                                                                          >  1:670585/1‑143 (MQ=255)
                                                                     aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAACTCGATTTGCCCCCTTCCCCATACCATCCCACCt                                                                          >  2:758158/1‑143 (MQ=255)
                                                                       cAACTCAGGTGAATCCGTCTTCGGCAGTTCCCGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTc                                                                        <  1:857774/143‑1 (MQ=255)
                                                                        atgtCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                       <  1:679995/140‑1 (MQ=255)
                                                                                   atccagcTTCGCCAGTTCCAGCGCCTACGCAATATTAGGCCACGTTAACCCTCTCCGCCCTTATAGGTTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑c                                                           <  1:925586/137‑1 (MQ=255)
                                                                                          ttCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGc                                                    <  2:1048293/143‑1 (MQ=255)
                                                                                            cGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCGAGTACTG‑CCTGTCGCTg                                                  >  1:646631/1‑143 (MQ=255)
                                                                                             ggCAGTTCCGGCGCCTACGCGGTATTTGGCGACGTTAAGCCCCGCCGCCGGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCGTGCGTCCCGTTCCCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGa                                                 <  2:1025416/143‑1 (MQ=255)
                                                                                              gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGAGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTCCTGCCAGCAATG‑CCTGTCGCTggc                                                >  2:432417/1‑141 (MQ=255)
                                                                                               cAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGCTAACGCTGCACGTCGATTTGCCCGCTTCCCGATACCATCAGACCTTCTGCCAGTGACG‑CCCGgcgccgcca                                               >  2:594114/1‑138 (MQ=255)
                                                                                                 gTTCCAGCTCCTACGCAAGATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCt                                             <  2:605200/143‑1 (MQ=255)
                                                                                                    ccAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTa                                          >  1:680919/1‑143 (MQ=255)
                                                                                                            tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACGCCTTCTGCCAGGAAGG‑CCTGTCGCTGACGCAGTAAAcacggt                                  >  1:622984/1‑143 (MQ=255)
                                                                                                              cGCAATATTAAGCGACGGTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAGCGCGGCAAGTCGAGTTGCCCGCTTCCCGATACCATCGCGCCCTCT‑CAAGATATGTTGTGTCGGTGACGCTGTAAAtacggtgc                                >  1:298394/1‑141 (MQ=255)
                                                                                                                cAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAAtccggtccgg                              >  1:699082/1‑143 (MQ=255)
                                                                                                                 aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTACCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATccggtccggt                             <  1:940138/143‑1 (MQ=255)
                                                                                                                     ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGc                         >  1:783999/1‑143 (MQ=255)
                                                                                                                         gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTCAGCGaaa                     >  2:993301/1‑143 (MQ=255)
                                                                                                                         gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGCAAATCCGGCCCGGTTAGccaca                     >  2:68713/1‑139 (MQ=255)
                                                                                                                         gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCAAGTAATG‑CCTGTCCCTGACGCTGTAAATCCGGTCCAGATAggcaac                     >  2:68734/1‑138 (MQ=255)
                                                                                                                           gACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                   <  2:655245/143‑1 (MQ=255)
                                                                                                                           gACGTTAACCCTCTCCGCCCGGAGAGGGGTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                   <  1:497350/143‑1 (MQ=255)
                                                                                                                            aCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGCAAGtg                  >  1:823159/1‑143 (MQ=255)
                                                                                                                             cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                 <  1:277588/143‑1 (MQ=255)
                                                                                                                                  aCCCTCTCCGCCCGGAGAGGGTTAGTAGTCCCAGCGCTCCCGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGt            <  1:38088/143‑1 (MQ=255)
                                                                                                                                    ccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTgg          >  1:680715/1‑143 (MQ=255)
                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATCACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTCCCGTGGc         <  2:460556/143‑1 (MQ=255)
                                                                                                                                     ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc         <  1:489141/143‑1 (MQ=255)
                                                                                                                                     ctctCCGCCCCGTGTTGGTTACTAGTTCCAGCTATTACGCTTCAGGTCGATTTGCCCGCTTCCCGATCCCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc         <  2:558393/143‑1 (MQ=255)
                                                                                                                                      tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTCCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATACGGTCCGGTTAGCGAAATTGCGCCGTGGCg        >  1:477845/1‑143 (MQ=255)
                                                                                                                                        gccGCCCGGAGAGGGTTAGCAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTATTG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAt      <  1:66985/142‑1 (MQ=255)
                                                                                                                                          cGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTACATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAATa    >  2:91719/1‑143 (MQ=255)
                                                                                                                                          cGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAACG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTCTGCCGTGCCGATTa    >  2:630386/1‑143 (MQ=255)
                                                                                                                                           gCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa   <  1:850335/143‑1 (MQ=255)
                                                                                                                                            cccGGAGAGGGCTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAGCATCACCCCTTCTGCCCGTAATG‑CCTGTCGCGGCCGCTGTAACGCCGGTCCGGTTAGCGCAACTGTGCCGTGGCGATGAac  >  2:67508/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATG‑CCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAC  >  NC_000913/475831‑476113

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