Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,137,180 C→G 16.7% A160P (GCA→CCA)  yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,1800CG16.7% 504.2 / 86.4 192A160P (GCA→CCA) yijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base C (84/75);  new base G (17/15);  total (102/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 197 total aligned reads displayed.

TAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGT  >  NC_000913/4137081‑4137279
                                                                                                   |                                                                                                   
 aaaCAGCACATCTTCTTCCTGCGCTCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGcc                                                                                                    <  1:1881970/100‑1 (MQ=255)
  aaCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCa                                                                                                   >  2:679587/1‑100 (MQ=255)
    cAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGc                                                                                                 <  2:2770749/100‑1 (MQ=255)
    cAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGc                                                                                                 <  1:679591/100‑1 (MQ=255)
      gCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGccc                                                                                               >  1:3823868/1‑100 (MQ=255)
        acaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGt                                                                                             <  2:2443029/100‑1 (MQ=255)
         caTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCg                                                                                              <  2:2522902/98‑1 (MQ=255)
         caTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                                                            >  2:1852881/1‑100 (MQ=255)
           tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTAt                                                                                          >  2:3809540/1‑100 (MQ=255)
           tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTAt                                                                                          <  1:72939/100‑1 (MQ=255)
           tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTAt                                                                                          >  1:330801/1‑100 (MQ=255)
           tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTAt                                                                                          <  2:3910649/100‑1 (MQ=255)
            cttcttCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATc                                                                                         <  1:1497554/100‑1 (MQ=255)
             ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATcc                                                                                        >  1:4117530/1‑100 (MQ=255)
              tcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCa                                                                                       <  1:3896401/100‑1 (MQ=255)
              tcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCa                                                                                       >  2:26690/1‑100 (MQ=255)
                  ccTGCGCGCCTGCGCCAATGTTCTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTc                                                                                    >  2:275455/1‑99 (MQ=255)
                  ccTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                                                                   >  2:1762878/1‑100 (MQ=255)
                   cTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGc                                                                                  <  1:2721159/100‑1 (MQ=255)
                   cTGCGCGCCTGCGCCAATGTTATGTATCACCAGCAGCGTACCGGCGCGGGCGAAGCGATCTGAAACCACCCCATTATGTGCCAGCCCGTTATCCAGTCGc                                                                                  >  2:885429/1‑100 (MQ=255)
                    tGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATGTGTGGCAGCCCGTTATCCAGTCGcc                                                                                 <  2:3521703/100‑1 (MQ=255)
                    tGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGcc                                                                                 >  2:190935/1‑100 (MQ=255)
                      cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAg                                                                               <  1:1269360/100‑1 (MQ=255)
                         gcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                                                            <  1:547720/100‑1 (MQ=255)
                         gcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                                                            <  1:1177058/100‑1 (MQ=255)
                          ccTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                                                           <  1:882929/100‑1 (MQ=255)
                             gcgcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAACTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACaa                                                                        >  2:595113/1‑100 (MQ=255)
                              cgcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAat                                                                       >  1:2031729/1‑100 (MQ=255)
                                ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAatat                                                                     >  2:2224116/1‑100 (MQ=255)
                                    tGTTATGTTTCACCAGCGGCGTACTGTCGCGGGCGAAGCGATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                                                                 <  1:2209372/100‑1 (MQ=255)
                                    tGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGCAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                                                                 >  2:3372929/1‑100 (MQ=255)
                                     gTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                                                                <  2:3074759/100‑1 (MQ=255)
                                      ttATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                                                               >  2:3179014/1‑100 (MQ=255)
                                         tGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCtt                                                            >  2:1479685/1‑100 (MQ=255)
                                         tGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATATGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATAGCGCCCGCtt                                                            >  2:425911/1‑100 (MQ=255)
                                          gTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTg                                                           <  1:4042523/100‑1 (MQ=255)
                                          gTACCACCAGCGGCGTACAATCGCGGGCGAAACCATCTGAAACCACCCCAATATGTGCCAACCCGTTATACAGTCGCCAAGAGACAATATc                                                                    >  2:493654/1‑91 (MQ=255)
                                           tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                                                          >  2:3303585/1‑100 (MQ=255)
                                             tCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa                                                        >  1:217625/1‑100 (MQ=255)
                                              cACCAGCTGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                       <  1:2243507/100‑1 (MQ=255)
                                              cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                       <  2:2980777/100‑1 (MQ=255)
                                               aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGt                                                      >  1:1013387/1‑100 (MQ=255)
                                                   gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACt                                                  <  2:405832/100‑1 (MQ=255)
                                                    cggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTg                                                 <  1:892636/100‑1 (MQ=255)
                                                       cgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg                                              <  2:4144165/100‑1 (MQ=255)
                                                        gTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGt                                             <  1:3800499/100‑1 (MQ=255)
                                                        gTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGt                                             <  2:3064354/100‑1 (MQ=255)
                                                         tACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg                                              <  2:3868790/98‑1 (MQ=255)
                                                         tACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGtt                                            <  2:796881/100‑1 (MQ=255)
                                                          acaggcgcgGGCGAAGCCAGCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg                                              <  2:3753568/92‑1 (MQ=255)
                                                           ccGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt                                          >  1:1828739/1‑100 (MQ=255)
                                                           ccGTCGCGGGCGAAGCCATATGAAACCACCCCAATATGTGGTAGCCCGGTATCCAGTCGCCAGCAGGCAATATCGGCCGCTTGATAGGCACTGGGCTTCt                                          >  2:118290/1‑100 (MQ=255)
                                                             gTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTg                                        <  2:576193/100‑1 (MQ=255)
                                                                gcgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtg                                     <  1:3858403/100‑1 (MQ=255)
                                                                gcgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtg                                     >  1:2842277/1‑100 (MQ=255)
                                                                 cgGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtgg                                    <  1:4174387/100‑1 (MQ=255)
                                                                 cgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtgg                                    <  1:3482362/100‑1 (MQ=255)
                                                                 cgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtgg                                    >  1:1248213/1‑100 (MQ=255)
                                                                  gggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggt                                   <  1:74948/100‑1 (MQ=255)
                                                                  gggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggt                                   >  2:1340856/1‑100 (MQ=255)
                                                                  gggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAAGATCGGCCTCTTGATAGTCACTGGTGTTCTTGCtggt                                   <  1:396444/100‑1 (MQ=255)
                                                                   ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggtg                                  >  1:806969/1‑100 (MQ=255)
                                                                   ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggtg                                  >  1:502693/1‑100 (MQ=255)
                                                                        aGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCg                             <  1:108855/100‑1 (MQ=255)
                                                                        aGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCg                             <  2:1013383/100‑1 (MQ=255)
                                                                         gCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgt                            >  1:214788/1‑100 (MQ=255)
                                                                         gCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgt                            >  1:2620496/1‑100 (MQ=255)
                                                                         gCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAAACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgt                            >  1:1581908/1‑100 (MQ=255)
                                                                          ccATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgtg                           <  1:3514468/100‑1 (MQ=255)
                                                                          ccATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgtg                           <  2:2151532/100‑1 (MQ=255)
                                                                          ccATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgt                            <  1:3071704/99‑1 (MQ=255)
                                                                               tGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                      >  2:567524/1‑100 (MQ=255)
                                                                                gAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTAt                     >  1:2553130/1‑100 (MQ=255)
                                                                                gAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTAt                     >  1:1310830/1‑100 (MQ=255)
                                                                                gAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTAt                     >  1:2876295/1‑100 (MQ=255)
                                                                                 aaaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATc                    <  2:4008718/100‑1 (MQ=255)
                                                                                  aaCCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCa                   >  2:1541928/1‑100 (MQ=255)
                                                                                   accaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCAt                  <  1:3737390/100‑1 (MQ=255)
                                                                                    ccaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCAt                  <  2:3160328/99‑1 (MQ=255)
                                                                                     caccCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGCCAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGGTTTATCATg                 >  2:1809466/1‑99 (MQ=255)
                                                                                     caccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATgg                >  1:3464326/1‑100 (MQ=255)
                                                                                     caccCCAATATGGGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATc                    <  2:3630285/96‑1 (MQ=255)
                                                                                      accCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGTCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggc               <  1:1674214/100‑1 (MQ=255)
                                                                                      accCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggc               <  1:1535649/100‑1 (MQ=255)
                                                                                       ccccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              >  1:3211727/1‑100 (MQ=255)
                                                                                       ccccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              >  1:2263626/1‑100 (MQ=255)
                                                                                       ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              >  1:2451592/1‑100 (MQ=255)
                                                                                       ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              >  1:222712/1‑100 (MQ=255)
                                                                                         ccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcggc            <  1:2955149/100‑1 (MQ=255)
                                                                                              atGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc       >  2:980766/1‑100 (MQ=255)
                                                                                              atGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc       >  2:3555219/1‑100 (MQ=255)
                                                                                               tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc      >  2:2240120/1‑100 (MQ=255)
                                                                                               tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc      <  1:1109742/100‑1 (MQ=255)
                                                                                               tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc      >  2:1295832/1‑100 (MQ=255)
                                                                                               tgggCCAGCCCGTTATCCAGTCGCCAGGAGACAATAGCGCCCCCTTCATCGGCACCGGGGTTCTTGCTGGTTGGGAGTGATTTATCCTtgcgggtgaggc      >  2:2942089/4‑93 (MQ=255)
                                                                                                gtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa     >  1:3219955/1‑100 (MQ=255)
                                                                                                gtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa     >  1:48420/1‑100 (MQ=255)
                                                                                                gtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa     >  1:10112/1‑100 (MQ=255)
                                                                                                  gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAgg   <  1:1750704/100‑1 (MQ=255)
                                                                                                  gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAgg   >  2:1689222/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGT  >  NC_000913/4137081‑4137279

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