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.

ATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767401‑2767594
                                                                                                     |                                                                                                
aTGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAttt                                                                                                   <  1:2062769/100‑1 (MQ=255)
   atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                                              <  2:1767351/100‑1 (MQ=255)
   atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                                              <  2:2077710/100‑1 (MQ=255)
     atTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccccc                                                                                           <  2:3601752/100‑1 (MQ=255)
     atTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccccc                                                                                           <  1:350844/100‑1 (MQ=255)
      ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccccc                                                                                          <  1:2952205/100‑1 (MQ=255)
            tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCcttt                                                                                     <  1:2621577/99‑2 (MQ=255)
              tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                   <  1:2000229/100‑1 (MQ=255)
              tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                   <  2:2562533/100‑1 (MQ=255)
               ccAATATGATCAAATAGCATTAACCCCCCCTCCCCACGTCATGCATAAGGAACACATTTTCCCCTGTGCACCCCCGACTAAATTccccccccccttaaac                                                                                 >  2:3009626/1‑92 (MQ=255)
               ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                  <  1:3320881/100‑1 (MQ=255)
               ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                  <  1:2396994/100‑1 (MQ=255)
                cAATATGATCAAATAGCATTAGCCCCCCCTCACAACGTCCTGCTTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACt                                                                                 <  2:3244936/100‑1 (MQ=255)
                cAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACCTTTTCCCCTGTGCACCCTCGACTAAATTTCCCTGCCCTTTAAACt                                                                                 >  2:1427295/1‑100 (MQ=255)
                 aaTATGATCAAATAGCATTAAACCCCCCGCACAACATCATGCATAGGGAACACGTTCTCCCCTCTGCCCCCCCGATTACATTTCACCCCCtcttccac                                                                                 >  2:3572622/1‑91 (MQ=255)
                    atGATCAAATGGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACGCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  2:2755277/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGTGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  2:2870816/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  2:1503391/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  1:865744/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  2:260244/100‑1 (MQ=255)
                    atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  1:1588630/100‑1 (MQ=255)
                     tGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCa                                                                            <  1:3373151/100‑1 (MQ=255)
                       aTCAAATAGCATTAACCCCCCCTCACAAAGTCCTGCAAATGAAACACATTTTCCACTGAGCACCCAAGACTCAATTTCACACCCCCATAAACCtctc                                                                             >  2:1887991/1‑97 (MQ=255)
                        tCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaata                                                                         <  1:5346/100‑1 (MQ=255)
                         cAAATAGCATTAACCCCCCCTCACACCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAATTTCTCaataa                                                                        >  2:4074840/1‑100 (MQ=255)
                         cAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaataa                                                                        <  2:3347494/100‑1 (MQ=255)
                          aaaTAGCATTAACCCCCCCTCACAACGTCCACCATAGGGAACATGTTTTCCCCTGTGCACCCACTACTAAATTACACAACCCTTTAAACAtcttcata                                                                         >  2:3340753/1‑93 (MQ=255)
                           aaTAGGATTACCCCCCCCCCACAACGTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTGCCCCCCCCTTTAAACTGCTCAATAATg                                                                      <  2:1737486/100‑1 (MQ=255)
                           aaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATg                                                                      <  1:1055464/100‑1 (MQ=255)
                             tAGCATTAACCCCCCCTCACATCGTCCGGCATTGGGAACACGTTCTCCCCTGTGCACCCACGACTAAATTTTCCCCCCCTTTACACTCCTCAGTAATgac                                                                    >  2:3904482/1‑98 (MQ=255)
                             tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                    <  1:1888048/100‑1 (MQ=255)
                             tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                    <  2:3886773/100‑1 (MQ=255)
                                 aTTAACCCGCCCTCACAACGTCCTGCATAGGGAACACGTTTGCCCCTGTGCACCCACGACGAAATTTCCCACCCCTCTAAACTTATCTATGATGTCagc                                                                 >  2:1448472/1‑97 (MQ=255)
                                 aTTAACCCCCCCTCACAACGTCCTTCCTAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTCCCCCATACTGTCAc                                                                  >  2:2831318/1‑98 (MQ=255)
                                 aTTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCa                                                                >  2:3189908/1‑100 (MQ=255)
                                 aTTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTAAACTTCTCaataaattaacca                                                                >  1:119563/1‑92 (MQ=255)
                                   tAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGAGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAAt                                                              >  2:2315167/1‑100 (MQ=255)
                                       ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                          <  1:1680357/100‑1 (MQ=255)
                                       ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTCAACCTTCTCAATAAAGTCACCaa                                                               >  1:3275940/1‑95 (MQ=255)
                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         >  2:1560156/1‑100 (MQ=255)
                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         >  1:477157/1‑100 (MQ=255)
                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         >  1:1195576/1‑100 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  1:76941/100‑1 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  1:1301378/100‑1 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  1:563030/100‑1 (MQ=255)
                                         ccccTCACAACGTCCTGCATAGGGAACACGTTTTCACCTTTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         <  2:1185023/99‑1 (MQ=255)
                                          cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                       <  1:939027/100‑1 (MQ=255)
                                          cccGCACATCGTCCTGCATAAGGAACACGTCTTCCCCTGTGCACCCACTACTAAAGTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                         <  2:3709832/98‑1 (MQ=255)
                                           ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                      <  1:100163/100‑1 (MQ=255)
                                           ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                      <  1:3684263/100‑1 (MQ=255)
                                           ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAAATTTCTgg                                                       >  1:3412740/1‑99 (MQ=255)
                                            cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                     >  2:1001556/1‑100 (MQ=255)
                                            cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                     >  2:509482/1‑100 (MQ=255)
                                             tCACAACGTCCTGCATAGGGAACACCTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTGCTCAATAATGTCACCAATTTTCTGGaaa                                                    <  2:410187/100‑1 (MQ=255)
                                                caACGTCCTGCATAGGGAACACGTTTTCCGCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAGTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 <  2:2511793/100‑1 (MQ=255)
                                                caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 <  1:3566703/100‑1 (MQ=255)
                                                   cggccTGCATAGGGATCCCGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTCAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                               <  2:876820/96‑1 (MQ=255)
                                                   cGTCCTGCATAGGGTACACGTTTTCCTCTGTCCACCCTCGACTAAATTTCCCCCCCCTTTAAACTTCTCAATTGTGGCACCAATTTTCTGGa                                                      <  2:1854481/92‑1 (MQ=255)
                                                   cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAAGAATGTCACCAATTTTCTGGAAATCCTa                                               <  2:782223/99‑1 (MQ=255)
                                                    gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTGTCTGGAAATCCTatg                                             >  1:4122555/1‑100 (MQ=255)
                                                     tCCTGCATAGGCAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTCTACCCTTCCCACTAATGCCACCCATCTTCTGGAAATCCTat                                              >  2:1843893/1‑98 (MQ=255)
                                                      ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                           <  1:1448483/100‑1 (MQ=255)
                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCCCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                           >  2:691676/1‑99 (MQ=255)
                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        >  2:16230/1‑100 (MQ=255)
                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         >  2:3607330/1‑99 (MQ=255)
                                                         cacaggggAAAACGTGTTCCCTGTGCACCCACGACTAAATTTCCCCCCCCGTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                       >  1:616082/6‑100 (MQ=255)
                                                           aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCCCCAATTTTCTGGAAATCCTATGATGCCtt                                      >  2:3229471/1‑100 (MQ=255)
                                                           aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt                                      <  1:2904309/100‑1 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:303102/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:410679/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  1:1876745/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:2608554/1‑100 (MQ=255)
                                                             aGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                    >  2:2753798/1‑100 (MQ=255)
                                                              gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                   <  2:861866/100‑1 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  <  2:318890/100‑1 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  >  1:2021181/1‑100 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  >  1:2135633/1‑100 (MQ=255)
                                                                 aaCACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                                >  1:3539463/1‑100 (MQ=255)
                                                                   cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                              >  1:3884970/1‑100 (MQ=255)
                                                                     cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                            <  2:1134335/100‑1 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           <  2:1061499/100‑1 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           >  2:3831412/1‑100 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATCTCACCAATTTCCTGGAAATCCTATAATGCCTTTCCCGGTaaa                            >  2:1683054/1‑99 (MQ=255)
                                                                       ttttCCCCTGTGCCCCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                          <  2:3674946/100‑1 (MQ=255)
                                                                          tCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                       >  2:1902383/1‑100 (MQ=255)
                                                                           ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAACCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGCTGCCTTTCCCGGTAAATGGTGc                      <  2:3963090/100‑1 (MQ=255)
                                                                           ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      >  1:4062571/1‑100 (MQ=255)
                                                                            cccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCt                     >  1:3267472/1‑100 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                    >  1:3710686/1‑100 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                    >  1:1414540/1‑100 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCt                     <  2:577433/99‑1 (MQ=255)
                                                                              cTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                   >  2:4154401/1‑100 (MQ=255)
                                                                                gtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGa                 >  2:3588866/1‑100 (MQ=255)
                                                                                gtgCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                <  1:3473411/100‑1 (MQ=255)
                                                                                   cACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt              >  1:2797412/1‑100 (MQ=255)
                                                                                     cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa            >  1:2572856/1‑100 (MQ=255)
                                                                                            cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     <  1:423091/100‑1 (MQ=255)
                                                                                            cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     <  1:1572056/100‑1 (MQ=255)
                                                                                              aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGTTTc   >  2:104736/1‑100 (MQ=255)
                                                                                              aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc   >  1:632333/1‑100 (MQ=255)
                                                                                                     |                                                                                                
ATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  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: