Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,507 (C)7→8 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.C91.2% 446.1 / 35.3 136intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (9/3);  new base C (57/67);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 7.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

ATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑C‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767401‑2767594
                                                                                                     |                                                                                                
aTGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑ttt                                                                                                   <  1:2062769/100‑1 (MQ=255)
   atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑c‑cc                                                                                              <  2:1767351/100‑1 (MQ=255)
   atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑c‑cc                                                                                              <  2:2077710/100‑1 (MQ=255)
     atTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑c‑cc‑cc                                                                                           <  2:3601752/100‑1 (MQ=255)
     atTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑c‑cc‑cc                                                                                           <  1:350844/100‑1 (MQ=255)
      ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑c‑cc‑ccc                                                                                          <  1:2952205/100‑1 (MQ=255)
            tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑C‑CC‑CCCCcttt                                                                                     <  1:2621577/99‑2 (MQ=255)
              tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTaaa                                                                                   <  1:2000229/100‑1 (MQ=255)
              tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTaaa                                                                                   <  2:2562533/100‑1 (MQ=255)
               ccAATATGATCAAATAGCATTAACCCCCCCTCCCCACGTCATGCATAAGGAACACATTTTCCCCTGTGCACCCCCGACTAAA‑TTc‑c‑cc‑ccccccttaaac                                                                                 >  2:3009626/1‑92 (MQ=255)
               ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAAc                                                                                  <  1:3320881/100‑1 (MQ=255)
               ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAAc                                                                                  <  1:2396994/100‑1 (MQ=255)
                cAATATGATCAAATAGCATTAGCCCCCCCTCACAACGTCCTGCTTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACt                                                                                 <  2:3244936/100‑1 (MQ=255)
                cAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACCTTTTCCCCTGTGCACCCTCGACTAAA‑TTT‑C‑CCTGCCCTTTAAACt                                                                                 >  2:1427295/1‑100 (MQ=255)
                 aaTATGATCAAATAGCATTAAACCCCCCGCACAACATCATGCATAGGGAACACGTTCTCCCCTCTGCCCCCCCGATTACA‑TTT‑C‑AC‑CCCCtcttccac                                                                                 >  2:3572622/1‑91 (MQ=255)
                    atGATCAAATGGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACGCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTtctc                                                                             <  2:2755277/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGTGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTtctc                                                                             <  2:2870816/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTtctc                                                                             <  2:1503391/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTtctc                                                                             <  1:865744/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTtctc                                                                             <  2:260244/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTtctc                                                                             <  1:1588630/100‑1 (MQ=255)
                     tGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCa                                                                            <  1:3373151/100‑1 (MQ=255)
                       aTCAAATAGCATTAACCCCCCCTCACAAAGTCCTGCAAATGAAACACATTTTCCACTGAGCACCCAAGACTCAA‑TTT‑CACA‑CCCCCATAAACCtctc                                                                             >  2:1887991/1‑97 (MQ=255)
                        tCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCaata                                                                         <  1:5346/100‑1 (MQ=255)
                         cAAATAGCATTAACCCCCCCTCACACCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAATTTCTCaataa                                                                        >  2:4074840/1‑100 (MQ=255)
                         cAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCaataa                                                                        <  2:3347494/100‑1 (MQ=255)
                          aaaTAGCATTAACCCCCCCTCACAACGTCCACCATAGGGAACATGTTTTCCCCTGTGCACCCACTACTAAA‑TTA‑CACA‑ACCCTTTAAACAtcttcata                                                                         >  2:3340753/1‑93 (MQ=255)
                           aaTAGGATTACCCCCCCCCCACAACGTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTGCC‑CC‑CCCCTTTAAACTGCTCAATAATg                                                                      <  2:1737486/100‑1 (MQ=255)
                           aaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATg                                                                      <  1:1055464/100‑1 (MQ=255)
                             tAGCATTAACCCCCCCTCACATCGTCCGGCATTGGGAACACGTTCTCCCCTGTGCACCCACGACTAAATTTT‑C‑CC‑CCCCTTTACACTCCTCAGTAATgac                                                                    >  2:3904482/1‑98 (MQ=255)
                             tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTc                                                                    <  1:1888048/100‑1 (MQ=255)
                             tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTc                                                                    <  2:3886773/100‑1 (MQ=255)
                                 aTTAACCCGCCCTCACAACGTCCTGCATAGGGAACACGTTTGCCCCTGTGCACCCACGACGAAA‑TTT‑C‑CCACCCCTCTAAACTTATCTATGATGTCagc                                                                 >  2:1448472/1‑97 (MQ=255)
                                 aTTAACCCCCCCTCACAACGTCCTTCCTAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTCCCCCATACTGTCAc                                                                  >  2:2831318/1‑98 (MQ=255)
                                 aTTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCa                                                                >  2:3189908/1‑100 (MQ=255)
                                 aTTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCCTTAAACTTCTCaataaattaacca                                                                >  1:119563/1‑92 (MQ=255)
                                   tAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGAGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAAt                                                              >  2:2315167/1‑100 (MQ=255)
                                       ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                          <  1:1680357/100‑1 (MQ=255)
                                       ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTCAACCTTCTCAATAAAGTCACCaa                                                               >  1:3275940/1‑95 (MQ=255)
                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         >  2:1560156/1‑100 (MQ=255)
                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         >  1:477157/1‑100 (MQ=255)
                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         >  1:1195576/1‑100 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  1:76941/100‑1 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  1:1301378/100‑1 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  1:563030/100‑1 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCACCTTTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         <  2:1185023/99‑1 (MQ=255)
                                          cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                       <  1:939027/100‑1 (MQ=255)
                                          cccGCACATCGTCCTGCATAAGGAACACGTCTTCCCCTGTGCACCCACTACTAAAGTTC‑C‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         <  2:3709832/98‑1 (MQ=255)
                                           ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                      <  1:100163/100‑1 (MQ=255)
                                           ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                      <  1:3684263/100‑1 (MQ=255)
                                           ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAAATTTCTgg                                                       >  1:3412740/1‑99 (MQ=255)
                                            cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                     >  2:1001556/1‑100 (MQ=255)
                                            cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                     >  2:509482/1‑100 (MQ=255)
                                             tCACAACGTCCTGCATAGGGAACACCTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTGCTCAATAATGTCACCAATTTTCTGGaaa                                                    <  2:410187/100‑1 (MQ=255)
                                                caACGTCCTGCATAGGGAACACGTTTTCCGCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAAGTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 <  2:2511793/100‑1 (MQ=255)
                                                caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 <  1:3566703/100‑1 (MQ=255)
                                                   cggccTGCATAGGGATCCCGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTCAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                               <  2:876820/96‑1 (MQ=255)
                                                   cGTCCTGCATAGGGTACACGTTTTCCTCTGTCCACCCTCGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATTGTGGCACCAATTTTCTGGa                                                      <  2:1854481/92‑1 (MQ=255)
                                                   cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAAGAATGTCACCAATTTTCTGGAAATCCTa                                               <  2:782223/99‑1 (MQ=255)
                                                    gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTGTCTGGAAATCCTatg                                             >  1:4122555/1‑100 (MQ=255)
                                                     tCCTGCATAGGCAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTCTACCCTTCCCACTAATGCCACCCATCTTCTGGAAATCCTat                                              >  2:1843893/1‑98 (MQ=255)
                                                      ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                           <  1:1448483/100‑1 (MQ=255)
                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCCCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                           >  2:691676/1‑99 (MQ=255)
                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        >  2:16230/1‑100 (MQ=255)
                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         >  2:3607330/1‑99 (MQ=255)
                                                         cacaggggAAAACGTGTT‑CCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCGTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                       >  1:616082/6‑100 (MQ=255)
                                                           aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCCCCAATTTTCTGGAAATCCTATGATGCCtt                                      >  2:3229471/1‑100 (MQ=255)
                                                           aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt                                      <  1:2904309/100‑1 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:303102/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:410679/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  1:1876745/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:2608554/1‑100 (MQ=255)
                                                             aGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                    >  2:2753798/1‑100 (MQ=255)
                                                              gggAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                   <  2:861866/100‑1 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  <  2:318890/100‑1 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  >  1:2021181/1‑100 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  >  1:2135633/1‑100 (MQ=255)
                                                                 aaCACGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                                >  1:3539463/1‑100 (MQ=255)
                                                                   cacGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                              >  1:3884970/1‑100 (MQ=255)
                                                                     cGTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                            <  2:1134335/100‑1 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           <  2:1061499/100‑1 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           >  2:3831412/1‑100 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATCTCACCAATTTCCTGGAAATCCTATAATGCCTTTCCCGGTaaa                            >  2:1683054/1‑99 (MQ=255)
                                                                       ttttCCCCTGTGCCCCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                          <  2:3674946/100‑1 (MQ=255)
                                                                          tCCCCTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                       >  2:1902383/1‑100 (MQ=255)
                                                                           ccccTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAACCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGCTGCCTTTCCCGGTAAATGGTGc                      <  2:3963090/100‑1 (MQ=255)
                                                                           ccccTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      >  1:4062571/1‑100 (MQ=255)
                                                                            cccTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCt                     >  1:3267472/1‑100 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                    >  1:3710686/1‑100 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                    >  1:1414540/1‑100 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCt                     <  2:577433/99‑1 (MQ=255)
                                                                              cTGTGCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                   >  2:4154401/1‑100 (MQ=255)
                                                                                gtgCACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGa                 >  2:3588866/1‑100 (MQ=255)
                                                                                gtgCACCCACGACTAAA‑TTT‑C‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                <  1:3473411/100‑1 (MQ=255)
                                                                                   cACCCACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt              >  1:2797412/1‑100 (MQ=255)
                                                                                     cccACGACTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa            >  1:2572856/1‑100 (MQ=255)
                                                                                            cTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     <  1:423091/100‑1 (MQ=255)
                                                                                            cTAAA‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     <  1:1572056/100‑1 (MQ=255)
                                                                                              aaa‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGTTTc   >  2:104736/1‑100 (MQ=255)
                                                                                              aaa‑TTTCC‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc   >  1:632333/1‑100 (MQ=255)
                                                                                                     |                                                                                                
ATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAA‑TTT‑C‑CC‑CCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767401‑2767594

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: ‑