Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,507:1 (C)7→8 100% intergenic (+152/‑203) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,5001.C83.7% 348.6 / 63.7 123intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (5/15);  new base C (54/49);  total (59/64)
Fisher's exact test for biased strand distribution p-value = 2.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

TGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTG  >  NC_000913/2767402‑2767595
                                                                                                    |                                                                                                    
tGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑c                                                                                                    <  1:2937303/100‑1 (MQ=255)
   tatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCTACGACTAAATTG‑‑cccc                                                                                              <  2:591115/100‑1 (MQ=255)
   tatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccc                                                                                              <  1:2010068/100‑1 (MQ=255)
    atTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑ccccc                                                                                             <  2:1442412/100‑1 (MQ=255)
     ttAACTTGTCCAATATGATCAAATAGCATTAAACAACCCTCCCAACCTCCTTCATAGAAAACACATTTTACCCTGTGCACCCCCGACTAAAGTT‑‑ccaccc                                                                                           >  2:3497505/1‑100 (MQ=255)
       aaCTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccccccc                                                                                         <  1:1191686/100‑2 (MQ=255)
       aaCTTGTCCAATATGATCAAATAGCATTAAAACCCCCCCACAAAATACTACATAAGGAACAAGTTATCCCCTGTACACACACAACCCAATCT‑‑cccccaca                                                                                          >  2:3112266/1‑97 (MQ=255)
          ttGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttt                                                                                      <  2:279195/100‑2 (MQ=255)
          ttGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAGTT‑‑CCCCCCCctt                                                                                       <  2:1180581/99‑1 (MQ=255)
           tGTCCAATATGATCAAATAGCATTATCCCCCCCTCACAACGTCCTTCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                     <  1:1164382/100‑3 (MQ=255)
           tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                     <  1:3818494/100‑3 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                   <  1:2147919/100‑1 (MQ=255)
               cAATATGATCAAATAGCATTAACCCAAACTCACAACGTACTGCATAGGGAACACGTTTTCCCTTGTCCACAAACAACTAAATAT‑‑CCCACCCCTTTaga                                                                                    >  2:182500/1‑96 (MQ=255)
                aaTATGATCAAATAGCATTAACCCCCCATAAAAACATCCTTCAAAAGGAACACATTTTCCCCTGTGCACCCACCACTAAATTT‑‑CCTCCCCCTATAAAc                                                                                   >  2:194339/1‑98 (MQ=255)
                 atatGATCAAATATCATTAACCCCCCCTCACAACGTCCTTCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtc                                                                                <  1:1347598/100‑1 (MQ=255)
                 atatGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtc                                                                                <  1:743122/100‑1 (MQ=255)
                   atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                              <  1:1389322/100‑1 (MQ=255)
                    tcaTCAAATAGCATTAGCGCCCCCTCAGAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTACGTTT‑‑CCCCCTCCTTTAAACTTCTCa                                                                             <  2:895389/98‑1 (MQ=255)
                    tGATCAAATAGCATTAACACCCCCTAACACAGTACTGCATAGGGAACAATTATACCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTATACTcctga                                                                             >  2:1614346/1‑98 (MQ=255)
                     gATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaa                                                                            <  1:350681/100‑1 (MQ=255)
                      aTCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGCGAACACGTTTTACCCTGTGCACCCACAACTAAATTTCCCCCCCCTTTCAACTTCTCaa                                                                            >  2:876383/1‑99 (MQ=255)
                       tCAAATAGCATTAGCTCCCCCTTACATTTTTTTTTTTAGTGAATACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaat                                                                           <  1:1745116/99‑1 (MQ=255)
                           atatCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATg                                                                       <  2:2510927/99‑1 (MQ=255)
                            tAGCATTAACCCCCCCTCACAACGTCCTGCATAGTGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                     <  1:1173137/100‑1 (MQ=255)
                            tAGCAATATCCCCCCCTCACCACGTCGTGAATAGGGAACACGTTTTCCCATTTGCACCCACGACTTAATTT‑‑CTCCCCACTTTAAACTTCAAaataa                                                                         <  2:1472128/96‑1 (MQ=255)
                              gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCAc                                                                   <  1:2873260/100‑1 (MQ=255)
                               cATTAACCCCCCCTCAAAACGCCCTGCATAGTGAACACATTTTTCACTATGCACCCACGCCTAAATTTCCCCCCCCTTTAAACTACTCtataat                                                                        >  2:1093319/1‑94 (MQ=255)
                               cATTAACACCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCAc                                                                   <  2:224105/99‑1 (MQ=255)
                                aTTAACCCCCACTCACAACGTCCTGCATAGGGAACACGTTGTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATTTCACCa                                                                 >  2:2024092/1‑100 (MQ=255)
                                       cccccTCACAACGTCCTGCATAGCGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                          <  2:3701773/100‑1 (MQ=255)
                                        ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                         <  1:3008930/100‑1 (MQ=255)
                                        ccccTCACAACGTCCTGCATAGGCAACACGTTTTCCCCTATGACCACACGACTAAATTT‑‑CCCACCCCTTTAAACTTATACACATTGTCACAAATTTTCt                                                          >  2:3464449/1‑99 (MQ=255)
                                         atcTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                         <  2:2307459/98‑1 (MQ=255)
                                             cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTTAAACTTCTCAATAATATCACCAATTTTCTgagat                                                     >  1:2307481/1‑96 (MQ=255)
                                              acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                   >  1:2356235/1‑100 (MQ=255)
                                              acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                   <  1:1279120/100‑1 (MQ=255)
                                              acaACGTCATGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACGAAATTT‑‑CACCCCCCCTTAAACTTCTCAATAATGTCACCACTTTTCTGGAAAt                                                    >  2:1247616/1‑99 (MQ=255)
                                               caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                  >  2:3174395/1‑100 (MQ=255)
                                               caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                  <  1:1214180/100‑1 (MQ=255)
                                                 acttCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTGTGAACTTCTCAATAATGTCATCGATTTTCTGGAGATCCTa                                                <  2:148490/97‑1 (MQ=255)
                                                 aCGTCCTGCATATGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                                <  2:1685803/100‑1 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                                >  1:3170272/1‑100 (MQ=255)
                                                  cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                               <  2:2631657/100‑1 (MQ=255)
                                                   gggtgcacaggggAAAACGTGTTCCCTATGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                             >  1:1178995/11‑100 (MQ=21)
                                                     ccTGCATAGTGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                               <  2:2220210/97‑1 (MQ=255)
                                                     ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                            >  2:3104611/1‑100 (MQ=255)
                                                     ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                            >  2:2864930/1‑100 (MQ=255)
                                                      cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                           <  1:3164764/100‑1 (MQ=255)
                                                      cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                           <  1:2120129/100‑1 (MQ=255)
                                                      cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                           >  2:104441/1‑100 (MQ=255)
                                                       tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                          >  1:1978525/1‑100 (MQ=255)
                                                         cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAATTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                        >  1:339290/1‑100 (MQ=255)
                                                           tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                      >  2:1182584/1‑100 (MQ=255)
                                                           tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt                                       >  1:3301569/1‑99 (MQ=255)
                                                             gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                    >  2:2401064/1‑100 (MQ=255)
                                                             gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                    >  2:2112087/1‑100 (MQ=255)
                                                             gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTAAGATGCCTTTcc                                    >  1:828028/1‑100 (MQ=255)
                                                              ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                   <  1:928262/100‑1 (MQ=255)
                                                               aaaCACGTTTTCCCCTATTCCCCCACCTCTAATTTT‑‑CCCCCCTTTTAAACTTTTCAATAATTTCACCAATTTTCTGTAAATCCTATTATGCCTTTCCCgg                                 <  2:2928711/99‑1 (MQ=255)
                                                                aaCACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                                 <  1:3514326/100‑1 (MQ=255)
                                                                aaCACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                <  2:716715/100‑1 (MQ=255)
                                                                  cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATATTGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                               <  2:3472953/100‑1 (MQ=255)
                                                                  cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATTTCACCAATTTTCAGGAAAACCTATAATGCCTTTCCCGGTa                               >  1:3558044/1‑100 (MQ=255)
                                                                    cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                             >  1:51856/1‑100 (MQ=255)
                                                                    cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                             >  2:3745392/1‑100 (MQ=255)
                                                                      ttttCCCCTGTGCACCCACGACTAAATTCC‑‑CCCCCCCTTTAAACTTCTCAATAATGTCAACAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                           <  2:2686486/100‑1 (MQ=255)
                                                                       tttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATCATGCCTTTCCCGGTAAATgg                          <  2:2392263/100‑1 (MQ=255)
                                                                       tttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTCCTGGAAATCCTATGATGCCTTTCCCGAAAAAt                            >  1:1636857/1‑98 (MQ=255)
                                                                          cttccgatctcccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                       <  2:197110/90‑1 (MQ=255)
                                                                          ccccTGTTCACCCTCGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                       <  1:1707954/100‑1 (MQ=255)
                                                                          ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTCTCGGTAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                       <  2:2417899/100‑1 (MQ=255)
                                                                            ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                     <  1:235217/100‑1 (MQ=255)
                                                                             cTGTACACCCCCGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAAGAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCTGTAAAGGGTGCtgtg                   <  2:523535/100‑1 (MQ=255)
                                                                              tgtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                   <  1:104444/100‑1 (MQ=255)
                                                                              tgtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                   <  1:1182597/100‑1 (MQ=255)
                                                                               gtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGa                  <  2:3394537/100‑1 (MQ=255)
                                                                                tgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                 <  1:3258594/100‑1 (MQ=255)
                                                                                tgCACCCACGACTAAATTTCCCCCCCCTCTAAACTTCTCAATAATGTCACCAATTTTATGGAAATCCTATAATGCCTTTACCGGTAAAAGGTGATGTGa                  >  1:1737345/1‑99 (MQ=255)
                                                                                 ncgcccACGACTAAATTTCCCCCCCCTTTAAACTTCCTCACAATGTCACCAATATTCTGGAAATCTTACAATGACTTAACCGGTAAAAAGTGCt                      >  1:1283048/4‑94 (MQ=255)
                                                                                 gCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                 >  1:222037/1‑99 (MQ=255)
                                                                                    cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATTCCTTTCCCGGTAAATGGTGCTGTGAGTTCa             >  2:3117764/1‑100 (MQ=255)
                                                                                    cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa             >  1:2818087/1‑100 (MQ=255)
                                                                                    cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAACCTATGATACCTTTCCCGGTAAATGGTGcagatcggaa              >  1:197115/1‑90 (MQ=255)
                                                                                      cACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa           >  2:3693514/1‑100 (MQ=255)
                                                                                      cACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa           >  2:815436/1‑100 (MQ=255)
                                                                                        cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa         <  1:32014/100‑1 (MQ=255)
                                                                                         gACTAAATTTCCCCCCCCTTTATACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCGTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg        <  1:1049125/100‑1 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTTTGAGTTCAGAAAGa       >  1:90187/1‑100 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa       >  2:2002849/1‑100 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa       >  1:1076943/1‑100 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGATAAATGGTGCTGTGAGTTCAGAAAGa       >  1:2517908/1‑100 (MQ=255)
                                                                                           cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATACTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt      >  1:1192271/1‑100 (MQ=255)
                                                                                           cTAAATTTCCCCCCCCTTTAAACTTATCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTAAGTTCAGAAAGAt      >  2:2193907/1‑100 (MQ=255)
                                                                                            tAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCCCCAATTTTCTGGAAATCCTATAATGCCTTTCCCGGTAAATGGTGATGTGAGTTCAGAAAGAtt     >  1:152221/1‑100 (MQ=255)
                                                                                             aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAtt     >  1:696621/1‑99 (MQ=255)
                                                                                             aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc    >  1:2676348/1‑100 (MQ=255)
                                                                                             aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCATTTCACGGTAAATGGTGCTGTTAGTTAAGAAAGATTc    >  2:1061370/1‑100 (MQ=255)
                                                                                              aaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATTCCTTTCCCGGTAAATGGTGCTATGAGTTAAGAAAGAt      >  1:3573092/1‑97 (MQ=255)
                                                                                              aaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCt   >  1:2777114/1‑100 (MQ=255)
                                                                                               aTTT‑‑CCCCCGTCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATTATGCCTTTCCCGTTAAATGGTTCTGTGAGTTCAGAAAGATTCTg  <  2:3309049/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTG  >  NC_000913/2767402‑2767595

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