Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,758,656 A→T 5.1% G402G (GGT→GGA 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,6560AT 5.1% 234.1 / 2.8 78G402G (GGT→GGAselBselenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base A (43/31);  new base T (2/2);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGCA  >  NC_000913/3758512‑3758804
                                                                                                                                                |                                                                                                                                                    
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                                                                                                        >  2:1130468/1‑149 (MQ=255)
                                          aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                        >  2:208500/1‑149 (MQ=255)
                                          aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                        >  1:1566969/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:1548502/1‑149 (MQ=255)
                                                                gTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAg                                                                                  >  2:239166/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                                                                    >  1:1346120/1‑149 (MQ=255)
                                                                              gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                                    >  2:666539/1‑149 (MQ=255)
                                                                              gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                                    >  1:1340039/1‑149 (MQ=255)
                                                                              gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                                    >  1:1325678/1‑149 (MQ=255)
                                                                              gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                                    >  2:393303/1‑149 (MQ=255)
                                                                              gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                                    >  2:124680/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                                                            >  1:414278/1‑131 (MQ=255)
                                                                                                        cAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                            <  2:414278/131‑1 (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                            >  1:529196/1‑149 (MQ=255)
                                                                                                                      aacaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgcgc                            <  2:1566969/149‑1 (MQ=255)
                                                                                                                        caaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                          <  1:1130468/149‑1 (MQ=255)
                                                                                                                        caaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                          >  2:230258/1‑149 (MQ=255)
                                                                                                                        caaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                          <  1:208500/149‑1 (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)
                                                                                                                                            tataACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCAt      >  2:1024427/1‑149 (MQ=255)
                                                                                                                                            tataACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCAt      >  2:63834/1‑149 (MQ=255)
                                                                                                                                               aaCCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGc   <  2:317950/149‑1 (MQ=255)
                                                                                                                                                aCCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGCa  >  2:1170553/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGCA  >  NC_000913/3758512‑3758804

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: