Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,758,651 A→T 6.4% I404N (ATT→AAT) ‡ selB ← selenocysteinyl‑tRNA‑specific translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,758,6510AT 6.4% 227.8 / 3.2 76I404N (ATT→AAT) ‡selBselenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base A (41/30);  new base T (3/2);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCC  >  NC_000913/3758503‑3758798
                                                                                                                                                    |                                                                                                                                                   
ttCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGaa                                                                                                                                                     >  2:865952/1‑149 (MQ=255)
ttCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGaa                                                                                                                                                     >  2:105740/1‑149 (MQ=255)
ttCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGaa                                                                                                                                                     >  2:484967/1‑149 (MQ=255)
ttCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGaa                                                                                                                                                     <  2:1142697/149‑1 (MQ=255)
         cAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTtatattcag                                                                                                                                            >  1:1397493/1‑144 (MQ=255)
         cAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAg                                                                                                                                            <  1:403489/149‑1 (MQ=255)
          aaCGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAgg                                                                                                                                           <  2:1448779/149‑1 (MQ=255)
            cGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCt                                                                                                                                         >  2:1495014/1‑149 (MQ=255)
            cGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGGCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCt                                                                                                                                         <  2:390222/149‑1 (MQ=255)
             gCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTg                                                                                                                                        <  2:1397493/149‑1 (MQ=255)
               cATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGtt                                                                                                                                      >  2:1087822/1‑149 (MQ=255)
                aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTg                                                                                                                                     >  2:203865/1‑149 (MQ=255)
                   cGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgcag                                                                                                                                  <  1:169181/149‑1 (MQ=255)
                        cAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAAtt                                                                                                                             >  1:178584/1‑149 (MQ=255)
                          gACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTcg                                                                                                                           <  1:1205437/149‑1 (MQ=255)
                             gTTCGCGCCCAGGGCCAGGTTAATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCAATAACCAGCCTGAATATAACAAGGCTGTTGCAGCAATTCtcgc                                                                                                                        >  1:423693/1‑149 (MQ=255)
                                     ccAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTc                                                                                                                <  1:1484130/149‑1 (MQ=255)
                                                   aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                  >  1:1566969/1‑149 (MQ=255)
                                                   aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                  >  2:208500/1‑149 (MQ=255)
                                                   aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                  >  2:1130468/1‑149 (MQ=255)
                                                    tCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGc                                                                                                              >  1:631703/1‑136 (MQ=255)
                                                    tCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGc                                                                                                              <  2:631703/136‑1 (MQ=255)
                                                      gcgATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgcgcg                                                                                               >  2:407866/1‑149 (MQ=255)
                                                           gTTGCTCATTATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCAAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGTCGCGCCCAg                                                                                          <  2:842311/149‑1 (MQ=255)
                                                            ttGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAgg                                                                                         <  2:1135533/149‑1 (MQ=255)
                                                                   tGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATc                                                                                  >  1:1449229/1‑149 (MQ=255)
                                                                     aTAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCg                                                                                >  1:697476/1‑149 (MQ=255)
                                                                        aGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCaa                                                                             <  1:1242645/149‑1 (MQ=255)
                                                                         gTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAg                                                                            >  2:239166/1‑149 (MQ=255)
                                                                         gTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAg                                                                            >  2:1548502/1‑149 (MQ=255)
                                                                          tCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAgg                                                                           >  1:317950/1‑149 (MQ=255)
                                                                           cGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGt                                                                          >  2:997458/1‑149 (MQ=255)
                                                                           cGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGt                                                                          >  1:445616/1‑149 (MQ=255)
                                                                            gCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGtt                                                                         >  1:456595/1‑149 (MQ=255)
                                                                              tAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATACCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTaa                                                                       >  2:679330/1‑149 (MQ=255)
                                                                               aaTGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAc                                                                      <  1:1423823/149‑1 (MQ=255)
                                                                                     tCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGTAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgcc                                                                >  2:709118/1‑149 (MQ=255)
                                                                                       gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                              >  2:393303/1‑149 (MQ=255)
                                                                                       gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                              >  1:1346120/1‑149 (MQ=255)
                                                                                       gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                              >  1:1340039/1‑149 (MQ=255)
                                                                                       gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                              >  1:1325678/1‑149 (MQ=255)
                                                                                       gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                              >  2:124680/1‑149 (MQ=255)
                                                                                       gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                              >  2:666539/1‑149 (MQ=255)
                                                                                         gaATTTTCCGCTGCCAGCGGGCTGCAACCGGCGCATTCAACAAGCTTTAACCAGCGTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAACTCCGCCTGGTTAACCGCGCCGCGt                                                            >  2:258977/1‑149 (MQ=255)
                                                                                          aaTTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGtt                                                           >  2:878274/1‑149 (MQ=255)
                                                                                           aTTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTc                                                          <  2:216381/149‑1 (MQ=255)
                                                                                           aTTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTc                                                          <  2:697476/149‑1 (MQ=255)
                                                                                            ttttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAgg                                                                                         >  1:309824/1‑117 (MQ=255)
                                                                                            ttttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAgg                                                                                         <  2:309824/117‑1 (MQ=255)
                                                                                               tCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                      <  2:178584/149‑1 (MQ=255)
                                                                                                ccGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGAt                                                     <  1:1163714/149‑1 (MQ=255)
                                                                                                   cTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACTAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGaa                                                  <  2:487730/149‑1 (MQ=255)
                                                                                                    tGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAAc                                                 <  2:1526685/149‑1 (MQ=255)
                                                                                                     gCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACa                                                <  2:592358/149‑1 (MQ=255)
                                                                                                        aGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGAt                                             <  1:1024898/149‑1 (MQ=255)
                                                                                                        aGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGAt                                             >  1:210052/1‑149 (MQ=255)
                                                                                                        aGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGAt                                             >  2:1216968/1‑149 (MQ=255)
                                                                                                          ggggcggCAACCGTCGCATTCAACAATATATAACCAGCCTTAATATAACCAGGCTGTTGCAGCAATTCTCGCATCCCTTCGCCGTTGAGCTGGCGCGCACAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATaa                                           <  2:423693/148‑1 (MQ=255)
                                                                                                               ggcAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcat                                      >  1:1245268/1‑149 (MQ=255)
                                                                                                                 cAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                      <  2:414278/131‑1 (MQ=255)
                                                                                                                 cAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                      >  1:414278/1‑131 (MQ=255)
                                                                                                                      ggCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGTAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggcatc                               <  1:709118/149‑1 (MQ=255)
                                                                                                                       gcgcATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggcatcg                              <  2:210052/149‑1 (MQ=255)
                                                                                                                            ttCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTg                         >  2:834521/1‑149 (MQ=255)
                                                                                                                            ttCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTg                         >  2:1098200/1‑149 (MQ=255)
                                                                                                                             tCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgc                        <  1:105740/149‑1 (MQ=255)
                                                                                                                             tCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgc                        <  1:791909/149‑1 (MQ=255)
                                                                                                                              caacaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgcg                       >  1:345785/1‑149 (MQ=255)
                                                                                                                               aacaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgcgc                      <  2:1566969/149‑1 (MQ=255)
                                                                                                                               aacaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgcgc                      >  1:529196/1‑149 (MQ=255)
                                                                                                                                 caaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                    <  1:1130468/149‑1 (MQ=255)
                                                                                                                                 caaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                    <  1:208500/149‑1 (MQ=255)
                                                                                                                                 caaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                    >  2:230258/1‑149 (MQ=255)
                                                                                                                                      tataACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCa               >  1:583358/1‑149 (MQ=255)
                                                                                                                                        taACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAag             >  1:767533/1‑149 (MQ=255)
                                                                                                                                         caacAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAaga            >  2:505704/4‑149 (MQ=255)
                                                                                                                                         caacAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCCAAATCCGCAAGGTTAACCGCGCCGCGGTCCAGATGAACAGATAACGCATCGCCATCGGTCTGCGCCCGTGCAaga            >  2:1480186/4‑149 (MQ=255)
                                                                                                                                         aaCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAaga            <  2:473839/149‑1 (MQ=255)
                                                                                                                                             aGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACg        <  1:673603/149‑1 (MQ=255)
                                                                                                                                                  gAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgccagc   >  2:1052339/1‑149 (MQ=255)
                                                                                                                                                   aaTATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgccagcc  <  1:1087822/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCC  >  NC_000913/3758503‑3758798

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