Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,614 A→C 15.5% intergenic (+41/‑51) yfjP → / → yfjQ CP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,768,6140AC15.5% 522.8 / 82.1 200intergenic (+41/‑51)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (84/84);  new base C (15/16);  total (100/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 205 total aligned reads displayed.

GGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTG  >  NC_000913/2768516‑2768713
                                                                                                  |                                                                                                   
ggCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAg                                                                                                    >  1:4194230/1‑100 (MQ=255)
ggCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAg                                                                                                    >  2:1986047/1‑100 (MQ=255)
ggCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAg                                                                                                    >  2:633592/1‑100 (MQ=255)
   cGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGc                                                                                                 <  2:3304292/100‑1 (MQ=255)
     tGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCgg                                                                                               >  2:3986663/1‑100 (MQ=255)
      gCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGc                                                                                              <  2:734807/100‑1 (MQ=255)
         ttCCTGAGTTGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                           <  1:2527002/100‑1 (MQ=25)
         ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                           <  1:466606/100‑1 (MQ=25)
         ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                           <  2:283767/100‑1 (MQ=25)
         ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                           <  2:2204878/100‑1 (MQ=25)
          tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                          <  2:1014729/100‑1 (MQ=255)
          tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAggg                                                                                          <  2:3518837/100‑1 (MQ=25)
           ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGc                                                                                         <  2:889474/100‑1 (MQ=255)
            cTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGct                                                                                        <  2:162923/100‑1 (MQ=25)
            cTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGct                                                                                        <  2:1846873/100‑1 (MQ=25)
            cTCGGTGGTGTCAGTGGCCCGCGCGGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGACCCCTGCCGCACTTGCGACAggg                                                                                          <  1:679061/98‑1 (MQ=21)
             tCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGctc                                                                                       <  1:1429107/100‑2 (MQ=25)
              cGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGctcg                                                                                      <  1:3256955/100‑3 (MQ=25)
               ggtggtGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGctcgg                                                                                     <  2:1377428/100‑4 (MQ=25)
                gtggtgTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                    <  2:3749150/100‑1 (MQ=255)
                gtggtgTCAGTGGCCCGCGCTGTCTGGGATTTCTTCGTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                    <  1:3863026/100‑1 (MQ=255)
                   gtgtCAGTGGCCCGCGCGGTCTGGGATTTCTTCTTCTGAGTGTTTACTTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTa                                                                                 >  1:1911855/1‑100 (MQ=255)
                    tgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                >  1:2384811/1‑100 (MQ=255)
                    tgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                >  2:4084902/1‑100 (MQ=255)
                    tgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                >  1:2477719/1‑100 (MQ=255)
                     gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAtt                                                                               <  1:2519917/100‑1 (MQ=255)
                     gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAtt                                                                               >  2:1209114/1‑100 (MQ=255)
                       cAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTg                                                                             <  2:3350677/100‑1 (MQ=255)
                       cAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTg                                                                             <  1:2131708/100‑1 (MQ=255)
                        aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAttt                                                                              <  2:2660902/98‑1 (MQ=255)
                          tGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGttt                                                                          <  2:3065519/100‑1 (MQ=255)
                           ggCCCGCGCTTTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGttt                                                                          <  1:1931181/99‑1 (MQ=255)
                           ggCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGtttt                                                                         <  2:992460/100‑1 (MQ=255)
                            gCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTa                                                                        >  1:1174917/1‑100 (MQ=255)
                              ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                      >  2:3269040/1‑100 (MQ=255)
                              ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                      <  1:1101294/100‑1 (MQ=255)
                                    tGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACt                                                                >  1:801063/1‑100 (MQ=255)
                                     gTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTg                                                               >  1:1223760/1‑100 (MQ=255)
                                       cTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                             >  1:2070042/1‑100 (MQ=255)
                                       cTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                             >  2:2582866/1‑100 (MQ=255)
                                        tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                            <  2:1270284/100‑1 (MQ=255)
                                         gggATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                            <  1:428158/99‑1 (MQ=255)
                                           gATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGt                                                         >  2:27403/1‑100 (MQ=255)
                                           gATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGt                                                         >  1:1061934/1‑100 (MQ=255)
                                           gATTTCTGCTTCCGAGTGTTTAATTGTGCTCCATTATTCTCCGAGTCGTGTCACAAGTGCGGCCGGGCTTTCTTTATTCGTTTTACCCATACTgtggag                                                          <  1:3665548/99‑1 (MQ=255)
                                            aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTc                                                        <  2:1911815/100‑1 (MQ=255)
                                            aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTc                                                        >  2:1639916/1‑100 (MQ=255)
                                               tcttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCGAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGc                                                     >  2:2390822/1‑100 (MQ=255)
                                                 ttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCtt                                                   >  1:415434/1‑100 (MQ=255)
                                                  tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTa                                                  <  1:4265680/100‑1 (MQ=255)
                                                  tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTa                                                  <  2:690265/100‑1 (MQ=255)
                                                   ctactggttgtCTAATTCTGTCCTTTTACTCTGCGAGTACTGTCGAAAGTGCGGCAGGGCTTTCTTTATTTGTTGTACCCTTACTGAGGTGTCTGCCTa                                                  <  1:1331301/91‑1 (MQ=255)
                                                    ttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATg                                                >  2:1597716/1‑100 (MQ=255)
                                                     tctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                               >  1:1958126/1‑100 (MQ=255)
                                                       tgggtgtTTAATTCTGCCCTATGATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAcc                                             <  1:40514/97‑1 (MQ=255)
                                                       tGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                               <  2:3882972/98‑1 (MQ=255)
                                                        gAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAccc                                            >  2:504748/1‑100 (MQ=255)
                                                            gtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCt                                        >  1:1632445/1‑100 (MQ=255)
                                                            gtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCt                                        >  2:2235802/1‑100 (MQ=255)
                                                              ttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                      <  1:3224861/100‑1 (MQ=255)
                                                               tAATTCTGCCCTATTATTCTCCGAGCCCTGTAACAAGTTCGGCAGGGCTTTCTTTAGTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                      >  1:2154831/1‑99 (MQ=38)
                                                                   tCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcg                                 >  1:2367019/1‑100 (MQ=255)
                                                                    cTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgc                                >  2:2928746/1‑100 (MQ=255)
                                                                     tGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgcg                               >  2:3808527/1‑100 (MQ=255)
                                                                     tGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgc                                <  1:3229030/99‑1 (MQ=255)
                                                                        ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCtt                            >  1:1013342/1‑100 (MQ=255)
                                                                         cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttt                           <  1:274214/100‑1 (MQ=255)
                                                                         cTATTATGCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttn                           <  2:3461316/100‑2 (MQ=255)
                                                                          tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttcg                          <  2:2993351/100‑3 (MQ=255)
                                                                          tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                          <  2:421794/100‑1 (MQ=255)
                                                                            ttattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgc                        <  2:2477731/100‑1 (MQ=255)
                                                                            ttattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgc                        <  2:3682056/100‑1 (MQ=255)
                                                                            ttattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgc                        >  2:1320230/1‑100 (MQ=255)
                                                                             tattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcg                       >  1:417193/1‑100 (MQ=255)
                                                                              attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcgc                      >  2:1494220/1‑100 (MQ=255)
                                                                                tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTCTTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTg                    <  1:166623/100‑1 (MQ=255)
                                                                                tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTg                    <  2:2131189/100‑1 (MQ=255)
                                                                                tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTg                    <  1:1320662/100‑1 (MQ=255)
                                                                                 ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                   >  1:1865821/1‑100 (MQ=255)
                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCa                  <  1:3297746/100‑1 (MQ=255)
                                                                                  tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCa                  <  1:1526810/100‑1 (MQ=255)
                                                                                  agaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCa                  >  2:937377/5‑100 (MQ=21)
                                                                                   gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaa                 >  2:117637/4‑100 (MQ=21)
                                                                                    cGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaaa                >  2:3238522/1‑100 (MQ=255)
                                                                                     aaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATAAGGAGGAGTCTGCTTATGACCCATCTGACTTCGCGCTTTGGCGCtcc                   >  1:948884/2‑94 (MQ=18)
                                                                                      aGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATc              >  1:3804417/1‑100 (MQ=255)
                                                                                       gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt             >  1:1858729/1‑100 (MQ=255)
                                                                                       gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt             >  1:3470707/1‑100 (MQ=255)
                                                                                       gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt             >  2:346597/1‑100 (MQ=255)
                                                                                       gCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt             >  2:863261/1‑100 (MQ=21)
                                                                                        cccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATCCTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCacctct             >  1:756549/1‑94 (MQ=18)
                                                                                        cccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCtt            >  2:1455155/1‑100 (MQ=21)
                                                                                          cTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt          >  2:4376682/1‑100 (MQ=255)
                                                                                          cTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt          <  2:2285013/100‑1 (MQ=255)
                                                                                            gTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTc        <  1:3765892/100‑1 (MQ=255)
                                                                                              cGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTc        <  1:2864355/98‑1 (MQ=255)
                                                                                              cGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGc      >  1:1711269/1‑100 (MQ=255)
                                                                                              cGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGc      >  1:2251961/1‑100 (MQ=255)
                                                                                                 aaGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGt   <  2:3588680/100‑1 (MQ=255)
                                                                                                  aGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTg  >  1:4326221/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTG  >  NC_000913/2768516‑2768713

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