Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 294,990 T→C H444R (CAC→CGC)  intF ← putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913294,9900TC99.6% 907.1 / ‑5.3 278H444R (CAC→CGC) intFputative phage integrase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (160/117);  total (161/117)
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.91e-01

Only 100 of 282 total aligned reads displayed.

TATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATA  >  NC_000913/294900‑295079
                                                                                          |                                                                                         
catGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         >  2:702728/2‑93 (MQ=255)
 atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGcc                                                                                        >  2:4019832/1‑93 (MQ=255)
 atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGcc                                                                                        <  1:6006290/93‑1 (MQ=255)
   tttGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAt                                                                                      <  2:198543/92‑1 (MQ=255)
   gTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAt                                                                                      <  2:1451766/93‑1 (MQ=255)
   gTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAt                                                                                      <  1:3386515/93‑1 (MQ=255)
   gTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTAGCGCCAt                                                                                      >  1:3241409/1‑93 (MQ=255)
    ttGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAtt                                                                                     >  1:5325816/1‑93 (MQ=255)
      gATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAtttt                                                                                   >  1:1068856/1‑93 (MQ=255)
      gATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAtttt                                                                                   >  1:96986/1‑93 (MQ=255)
      gATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAtttt                                                                                   <  2:2608666/93‑1 (MQ=255)
        tttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCg                                                                                 <  2:4729540/93‑1 (MQ=255)
        tttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCg                                                                                 <  2:2365285/93‑1 (MQ=255)
         ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGt                                                                                >  2:3685143/1‑93 (MQ=255)
           ttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTaa                                                                              <  2:6682791/93‑1 (MQ=255)
           ttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTaa                                                                              <  1:2875011/93‑1 (MQ=255)
            tAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAAc                                                                             >  2:4045084/1‑93 (MQ=255)
             aaTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACa                                                                            >  1:849078/1‑93 (MQ=255)
             aaTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACa                                                                            >  1:1478585/1‑93 (MQ=255)
             aaTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACa                                                                            >  1:3759476/1‑93 (MQ=255)
              aTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAg                                                                           >  1:1854119/1‑93 (MQ=255)
              aTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAg                                                                           >  1:1466044/1‑93 (MQ=255)
                ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGAt                                                                         >  2:3293805/1‑93 (MQ=255)
                ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGAt                                                                         <  1:5539766/93‑1 (MQ=255)
                 gAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATc                                                                        >  2:447827/1‑93 (MQ=255)
                  aaTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCt                                                                       >  2:3134199/1‑93 (MQ=255)
                  aaTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCt                                                                       <  1:2116377/93‑1 (MQ=255)
                  aaTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCt                                                                       >  1:1611594/1‑93 (MQ=255)
                   aTCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTa                                                                      >  2:3204171/1‑93 (MQ=255)
                    tCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTaa                                                                     >  2:2935635/1‑93 (MQ=255)
                    tCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTaa                                                                     <  2:6185788/93‑1 (MQ=255)
                     cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAAc                                                                    <  1:3535527/93‑1 (MQ=255)
                     cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAAc                                                                    <  1:4913755/93‑1 (MQ=255)
                      aaCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACg                                                                   >  2:5892438/1‑93 (MQ=255)
                        cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGga                                                                 <  2:6380283/93‑1 (MQ=255)
                         gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgac                                                                <  2:2314733/93‑1 (MQ=255)
                           atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacga                                                              >  2:831230/1‑93 (MQ=255)
                           atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacga                                                              >  1:432215/1‑93 (MQ=255)
                            tatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacgac                                                             <  2:3576076/93‑1 (MQ=255)
                               cAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGc                                                          <  2:6282849/93‑1 (MQ=255)
                                aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCg                                                         >  1:5751141/1‑93 (MQ=255)
                                 aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGa                                                        <  1:4689698/93‑1 (MQ=255)
                                   gttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                      >  1:294029/1‑93 (MQ=255)
                                   gttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                      >  1:1900048/1‑93 (MQ=255)
                                    ttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAg                                                     >  1:1339985/1‑93 (MQ=255)
                                     tgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgt                                                    <  2:426802/93‑1 (MQ=255)
                                      gttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtg                                                   <  2:336010/93‑1 (MQ=255)
                                      gttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtg                                                   <  2:5126149/93‑1 (MQ=255)
                                          tttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                               <  1:5039749/93‑1 (MQ=255)
                                          tttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                               <  1:3626937/93‑1 (MQ=255)
                                           ttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGttttt                                              >  2:1528486/1‑93 (MQ=255)
                                            tATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                             <  2:4312659/93‑1 (MQ=255)
                                            tATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                             <  2:209446/93‑1 (MQ=255)
                                             aTGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCg                                            >  1:4794500/1‑93 (MQ=255)
                                                gTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCa                                         >  1:2494282/1‑93 (MQ=255)
                                                  aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAag                                       <  2:4672965/93‑1 (MQ=255)
                                                   aTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                      <  2:1205093/93‑1 (MQ=255)
                                                   aTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                      <  2:1336823/93‑1 (MQ=255)
                                                   aTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                      >  2:1775538/1‑93 (MQ=255)
                                                   aTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGCTTTTCGGTAaga                                      >  2:5715239/1‑93 (MQ=255)
                                                    tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                     >  1:2625463/1‑93 (MQ=255)
                                                    tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                     <  2:2888607/93‑1 (MQ=255)
                                                    tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                     >  1:3469801/1‑93 (MQ=255)
                                                     tccTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgc                                    <  1:1113042/92‑1 (MQ=255)
                                                     aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgc                                    >  1:794065/1‑93 (MQ=255)
                                                     aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgc                                    >  1:258892/1‑93 (MQ=255)
                                                       cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                  >  1:5635436/1‑93 (MQ=255)
                                                       cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                  >  2:2995007/1‑93 (MQ=255)
                                                       cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                  >  2:2255587/1‑93 (MQ=255)
                                                       cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                  <  2:1095706/93‑1 (MQ=255)
                                                        tGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                 >  1:2083018/1‑93 (MQ=255)
                                                         gCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGTCAAGAGCGCTg                                >  2:3776339/1‑93 (MQ=255)
                                                         gCTTCATTAAAGATCCATGCCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTg                                >  1:6506869/1‑93 (MQ=255)
                                                          cTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTgg                               <  1:1788245/93‑1 (MQ=255)
                                                           ttCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTggg                              >  1:2340128/1‑93 (MQ=255)
                                                           ttCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTggg                              >  2:5333709/1‑93 (MQ=255)
                                                            tCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGt                             <  1:397173/93‑1 (MQ=255)
                                                              aTTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGttt                           >  1:3324903/1‑93 (MQ=255)
                                                              aTTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGttt                           >  1:5872359/1‑93 (MQ=255)
                                                              aTTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGttt                           >  2:2828259/1‑93 (MQ=255)
                                                                ttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                         >  1:6718570/1‑93 (MQ=255)
                                                                   aGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTc                      <  2:4428455/93‑1 (MQ=255)
                                                                   aGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTc                      <  2:1475653/93‑1 (MQ=255)
                                                                    gATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTcc                     >  2:5282258/1‑93 (MQ=255)
                                                                      tCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCa                   <  2:557535/93‑1 (MQ=255)
                                                                        cATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                 <  1:2789373/93‑1 (MQ=255)
                                                                        cATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                 >  1:5953179/1‑93 (MQ=255)
                                                                         aTGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                >  2:4457626/1‑93 (MQ=255)
                                                                          tGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAt               >  1:4884190/1‑93 (MQ=255)
                                                                           gTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAtt              >  1:4366338/1‑93 (MQ=255)
                                                                            tctcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATATTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAttt             >  2:3954286/1‑93 (MQ=255)
                                                                            tctcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAttt             <  1:2754408/93‑1 (MQ=255)
                                                                            tctcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAttt             <  1:1146317/93‑1 (MQ=255)
                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCCGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           >  2:1609520/1‑93 (MQ=255)
                                                                                aaTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGc         >  1:5901016/1‑93 (MQ=255)
                                                                                    ttCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca     <  1:6308251/93‑1 (MQ=255)
                                                                                    ttCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaaca     <  1:2163063/93‑1 (MQ=255)
                                                                                      ctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAAt   <  2:4904333/93‑1 (MQ=255)
                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  >  1:4087215/1‑93 (MQ=255)
                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  <  1:1249147/93‑1 (MQ=255)
                                                                                          |                                                                                         
TATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATA  >  NC_000913/294900‑295079

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