Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,616,491 T→C 48.1% D192G (GAT→GGT)  yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,4910TC48.1% ‑5.4 / 269.0 194D192G (GAT→GGT) yjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (51/50);  new base C (38/55);  total (89/105)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 194 total aligned reads displayed.

TTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCG  >  NC_000913/4616393‑4616589
                                                                                                  |                                                                                                  
 tttATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCg                                                                                                  <  2:3205491/100‑1 (MQ=255)
 tttATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCg                                                                                                  <  1:1722069/100‑1 (MQ=255)
  ttATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGt                                                                                                 >  2:580796/1‑100 (MQ=255)
   tATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTa                                                                                                <  2:730954/100‑1 (MQ=255)
    aTCATGCACGGGGCCGTTGGCGAGGTGCGGTTTGCTACATTCTCAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAg                                                                                               <  2:1140165/100‑1 (MQ=255)
    aTCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAg                                                                                               <  1:3443365/100‑1 (MQ=255)
    aTCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAg                                                                                               <  1:156201/100‑1 (MQ=255)
      cATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGAt                                                                                             >  1:2001056/1‑100 (MQ=255)
       aTGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATa                                                                                            >  1:3200523/1‑100 (MQ=255)
       aTGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATa                                                                                            <  1:1147352/100‑1 (MQ=255)
       aTGCACCGGAACGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGa                                                                                              >  2:3445872/1‑98 (MQ=255)
        tGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTTACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATaa                                                                                           <  2:2781788/100‑1 (MQ=255)
           aCCGGACTGTTGGCGATGTGCGGTTTGCTACATCCACAGAGGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCATGATCGTAGATAAAgg                                                                                        <  2:1369469/100‑1 (MQ=255)
           aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAgg                                                                                        <  1:3480852/100‑1 (MQ=255)
           aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAgg                                                                                        <  1:2132677/100‑1 (MQ=255)
           aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAg                                                                                         >  2:3003288/1‑99 (MQ=255)
               gACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGCTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAg                                                                                    <  2:457210/100‑1 (MQ=255)
                aCCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAgg                                                                                   >  1:441643/1‑100 (MQ=255)
                     tGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCgg                                                                              <  1:3053506/100‑1 (MQ=255)
                     tGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCgg                                                                              <  1:1792936/100‑1 (MQ=255)
                     tGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCgg                                                                              <  2:869735/100‑1 (MQ=255)
                      ggCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGt                                                                             <  2:350795/100‑1 (MQ=255)
                      ggCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGt                                                                             <  1:528875/100‑1 (MQ=255)
                      ggCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGt                                                                             >  2:1471763/1‑100 (MQ=255)
                       gCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTg                                                                            <  1:1035725/100‑1 (MQ=255)
                            gtgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAAcc                                                                       <  2:3786362/100‑1 (MQ=255)
                            gtgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAAcc                                                                       <  1:2696918/100‑1 (MQ=255)
                               cGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGc                                                                    <  1:1433195/100‑1 (MQ=255)
                                  tttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTc                                                                 <  1:134526/100‑1 (MQ=255)
                                     gCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAAc                                                              <  1:1583469/100‑1 (MQ=255)
                                     gCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAAc                                                              >  2:485198/1‑100 (MQ=255)
                                       tACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGTTCATCAGGGGTTATTTCAGGATCGTAGATAAAGGACAGGTTCGGTGTACCCTGCTTCAActc                                                            >  2:1974180/1‑100 (MQ=255)
                                       tACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCGGGGGTGATTTCAGGATCGGAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAAct                                                             <  2:3409402/99‑1 (MQ=255)
                                       tACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAActc                                                            >  2:3998757/1‑100 (MQ=255)
                                       tACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCACctc                                                            >  2:3776910/1‑100 (MQ=255)
                                        aCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAActct                                                           >  2:3167566/1‑100 (MQ=255)
                                         cATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctg                                                          <  1:4015437/100‑1 (MQ=255)
                                          aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgc                                                         <  1:2582485/100‑1 (MQ=255)
                                          aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgc                                                         <  2:3081061/100‑1 (MQ=255)
                                           ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca                                                        <  1:1055284/100‑1 (MQ=255)
                                           ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca                                                        >  1:1404536/1‑100 (MQ=255)
                                           ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca                                                        <  1:37333/100‑1 (MQ=255)
                                            tCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcat                                                       >  1:3782111/1‑100 (MQ=255)
                                             cacaGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatc                                                      <  2:3819846/100‑1 (MQ=255)
                                              acaGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatct                                                     <  2:3929943/100‑1 (MQ=255)
                                               caGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctg                                                    <  2:3630346/100‑1 (MQ=255)
                                               caGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctg                                                    <  2:2914849/100‑1 (MQ=255)
                                               caGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctg                                                    <  2:3481908/100‑1 (MQ=255)
                                                aGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctgc                                                   <  1:3764873/100‑1 (MQ=255)
                                                 gATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctgca                                                  <  1:2879270/100‑1 (MQ=255)
                                                 gATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctgca                                                  <  1:1500814/100‑1 (MQ=255)
                                                 gATGTTCTTCGCCACTGCCAGGAGCAGGTGATCAGGGGTGATTTGAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgcatctgca                                                  <  2:2001043/100‑1 (MQ=255)
                                                  aTGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCAc                                                 <  2:4137/100‑1 (MQ=255)
                                                   tGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACg                                                >  1:1920376/1‑100 (MQ=255)
                                                    gTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGt                                               >  1:2386932/1‑100 (MQ=255)
                                                    gTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGt                                               >  1:1836005/1‑100 (MQ=255)
                                                       cttcGCCACTTCCAGCAGTAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCAGGTAAg                                            >  2:676231/1‑100 (MQ=255)
                                                       cttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAg                                            >  2:2609931/1‑100 (MQ=255)
                                                       cttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAg                                            <  1:905677/100‑1 (MQ=255)
                                                         tggCCCCTTCCAGCCGCAGGTCATCGGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTCCGGCGAAACCTGCTTCATCtctgcacctgca                                                  <  2:405079/90‑1 (MQ=255)
                                                          cGCCACTTCCAGCAGTAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATc                                         <  1:4063859/100‑1 (MQ=255)
                                                          cGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATc                                         >  1:66262/1‑100 (MQ=255)
                                                           gCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATcg                                        >  1:489033/1‑100 (MQ=255)
                                                               cTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATcgcgcg                                    <  1:299773/100‑1 (MQ=255)
                                                                 tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg                                  <  2:2468971/100‑1 (MQ=255)
                                                                 tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg                                  <  2:2068197/100‑1 (MQ=255)
                                                                 tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg                                  <  1:624447/100‑1 (MQ=255)
                                                                 tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg                                  <  1:1652416/100‑1 (MQ=255)
                                                                  ccAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGt                                 >  1:2810983/1‑100 (MQ=255)
                                                                   cagcagcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTa                                <  1:986948/100‑1 (MQ=255)
                                                                     gcagcagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAAt                              <  1:2517177/100‑1 (MQ=255)
                                                                     gcagcagGTCATCAGGGGTGATATCAGGATCGTAGATAAAGATCAGGCTCGGTGAAACCTGGTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAAt                              >  2:2254050/1‑100 (MQ=255)
                                                                      cagcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATg                             >  2:402181/1‑100 (MQ=255)
                                                                       agcagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATgg                            <  2:1841025/100‑1 (MQ=255)
                                                                       agcagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTAAGGTTCGGTGAAACCTGCTTCTACTCTTTATCTGCACGTAAGATCGCGCGGGTAATgg                            <  2:143692/100‑1 (MQ=255)
                                                                        gcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATggg                           <  2:3370103/100‑1 (MQ=255)
                                                                         cagGTCATCAGGGGTGATGTCATGACCCTAGCTGCAGGTCAGGCTCGGTGGAACCTGCTTAAACTTTGCATCTGCACGTAAGATCGCGCGGGt                                 <  2:2852124/93‑1 (MQ=255)
                                                                          agGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCg                         <  2:3639328/100‑1 (MQ=255)
                                                                           ggTCCTCAGGGGTGATTTCAGGACCGTAGCTAAAGGTCAGGTTCGGTGCCCCCTGCTTCAACTCTGCATCTGCACGTCACATCGCGCGGGTAATGCGCGa                        >  2:2455650/1‑100 (MQ=255)
                                                                           ggTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGa                        >  2:3791635/1‑100 (MQ=255)
                                                                            gTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGaa                       <  1:1469914/100‑1 (MQ=255)
                                                                             tcatcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAAt                      <  2:2303202/100‑1 (MQ=255)
                                                                              catcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                     >  2:931713/1‑100 (MQ=255)
                                                                              catcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                     >  2:1086104/1‑100 (MQ=255)
                                                                               atcaGGGGTGATTTCAGGACCGTAGATAAAGGTCCGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCa                    <  2:2179041/100‑1 (MQ=255)
                                                                               atcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCa                    <  1:3736931/100‑1 (MQ=255)
                                                                                tcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAg                   >  2:3567121/1‑100 (MQ=255)
                                                                                 caGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGa                  <  2:3988700/100‑1 (MQ=255)
                                                                                   ggggTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACg                <  1:46590/100‑1 (MQ=255)
                                                                                   ggggTGATTTCAGGACCGTAGATAAAGGTCAGCTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACg                >  1:3923314/1‑100 (MQ=255)
                                                                                    gggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACgg               >  1:88763/1‑100 (MQ=255)
                                                                                     ggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGAGGGTAATGGGAGAATCAGACgg               >  2:1165946/1‑99 (MQ=255)
                                                                                       tGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGcc            <  1:2218876/100‑1 (MQ=255)
                                                                                       tGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGcc            <  1:809981/100‑1 (MQ=255)
                                                                                          tgtctggATCGTAGATAAAAGTCAGGTAGGGTGAAACCTGCTTCAACGCTAGACCTGGACGTAAAATCGAGCGGGTAATGGGCAAATCAGACGGGCCGa          <  2:3515651/94‑1 (MQ=255)
                                                                                           ttCAAGATCGTAGATAAAGGTCAGGTTCGGTGAAGCCTGCTTTAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCg           <  2:1150696/97‑1 (MQ=255)
                                                                                             cAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATAtt      >  1:1306372/1‑100 (MQ=255)
                                                                                               ggATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTgg    >  1:3244701/1‑100 (MQ=255)
                                                                                               ggACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTgg    <  2:1567054/100‑1 (MQ=255)
                                                                                                 aTCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCg  <  1:1260438/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCG  >  NC_000913/4616393‑4616589

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