Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,613,051 A→G 5.2% Y214C (TAC→TGC)  yjjU → putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,613,0510AG 5.2% 237.8 / 3.0 77Y214C (TAC→TGC) yjjUputative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (2/2);  total (39/38)
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.33e-01

GCGATACCTGGCTTTTATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGCG  >  NC_000913/4612903‑4613192
                                                                                                                                                    |                                                                                                                                             
gCGATACCTGGCTTTTATCGTAGCGGAGGGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGta                                                                                                                                               <  6:30459/149‑1 (MQ=255)
 cGATACCTGGCTTTTATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGta                                                                                                                                               <  7:74354/148‑1 (MQ=255)
       cTGGCTTTTATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACAcg                                                                                                                                        >  4:170544/1‑149 (MQ=255)
        tGGCTTTTATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACAcgc                                                                                                                                       >  3:169635/1‑149 (MQ=255)
          gCTGTTATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACAcgccg                                                                                                                                     <  4:210095/149‑1 (MQ=255)
            ttttATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACAcgccg                                                                                                                                     >  3:234951/1‑147 (MQ=255)
                aTCGTTGCGGAGTGTCGCTGGAAGGCATTTACTACCTGGATGGCGGGATCAGTGATTCGATTCCGGTTAAAGAGGCTGAAAGGCTGGTCGCTAAAACGTTGGTCTTCATTCTAAATTTGCCGTCACAAATGTACTACACGCCGCAGTg                                                                                                                                <  2:33101/148‑1 (MQ=255)
                aTCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGTTCAGTGATTCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGAACTACACGCAGCAGTgg                                                                                                                               >  3:262431/1‑149 (MQ=255)
                aTCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTgg                                                                                                                               >  7:21263/1‑149 (MQ=255)
                  cGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGt                                                                                                                              <  6:287256/148‑1 (MQ=255)
                   gTAGCGGAGTGTCGCTGGAAGGCATTATCTACCTGGATGGCGGGATCATTGATTCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGATGGTCGGCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTc                                                                                                                            <  4:262431/149‑1 (MQ=255)
                   gTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGTATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGGCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGtt                                                                                                                             <  6:255661/148‑1 (MQ=255)
                   gTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGtt                                                                                                                             <  3:170544/148‑1 (MQ=255)
                      gCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCaaa                                                                                                                         <  3:25825/149‑1 (MQ=255)
                             gtCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATg                                                                                                                   >  1:220569/1‑148 (MQ=255)
                              tCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATgg                                                                                                                  >  3:307329/1‑148 (MQ=255)
                                 cTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTGCTACACGCCGCAGTGGTTCAAACGCATGGAACg                                                                                                              >  4:47325/1‑149 (MQ=255)
                                 cTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTGCTACACGCCGCAGTGGTTCAAACGCATGGAACg                                                                                                              >  4:52187/1‑149 (MQ=255)
                                         cATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGctggctgg                                                                                                      <  3:137876/149‑1 (MQ=255)
                                          aTTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTggg                                                                                                     >  6:124895/1‑149 (MQ=255)
                                          aTTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTggg                                                                                                     >  7:391757/1‑149 (MQ=255)
                                           ttAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGt                                                                                                    <  1:366252/149‑1 (MQ=255)
                                             aaCTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCAGAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGa                                                                                                  >  1:266416/1‑149 (MQ=255)
                                                  ccTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGt                                                                                              <  8:145073/148‑1 (MQ=255)
                                                  ccTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGt                                                                                              <  6:97879/148‑1 (MQ=255)
                                                  ccTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTa                                                                                             <  4:307329/149‑1 (MQ=255)
                                                      gATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACTAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTCGCCt                                                                                         >  2:336015/1‑149 (MQ=255)
                                                      gATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCt                                                                                         <  4:169635/149‑1 (MQ=255)
                                                      gATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCt                                                                                         >  7:118702/1‑149 (MQ=255)
                                                       aTGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTg                                                                                        <  7:378061/149‑1 (MQ=255)
                                                         ggCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCCAACGCATGGAACGCTGGCTGGGTGACCGTAGCCTg                                                                                        >  5:276467/1‑147 (MQ=255)
                                                             ggAGCAGTGATGCGATTCCGGGTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTGGGGCGTCATTCGCACTGTGCCGTCACAAATGGACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCg                                                                                  <  6:244763/149‑1 (MQ=255)
                                                                tCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTg                                                                               <  8:391757/149‑1 (MQ=255)
                                                                tCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGACGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTg                                                                               >  1:324653/1‑149 (MQ=255)
                                                                      aTGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAAt                                                                         >  5:169310/1‑149 (MQ=255)
                                                                       tGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATc                                                                        >  6:234330/1‑149 (MQ=255)
                                                                         cGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTg                                                                      <  7:213232/149‑1 (MQ=255)
                                                                               cGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTgca                                                                 <  4:220335/148‑1 (MQ=255)
                                                                               cGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTgca                                                                 <  4:274024/148‑1 (MQ=255)
                                                                                    aaaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTGCTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatc                                                            <  3:47325/148‑1 (MQ=255)
                                                                                    aaaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTGCTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatc                                                            <  3:52187/148‑1 (MQ=255)
                                                                                    aaaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatca                                                           >  6:378010/1‑149 (MQ=255)
                                                                                    aaaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatc                                                            <  2:191763/148‑1 (MQ=255)
                                                                                    aaaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatc                                                            <  5:234330/148‑1 (MQ=255)
                                                                                     aaGAGGCGGCAAGGCAGGGCGCTAAAACGTTTGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatcat                                                          <  1:48715/149‑1 (MQ=255)
                                                                                     aaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACgct                                                                                                            <  7:287726/99‑1 (MQ=255)
                                                                                     aaGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACgct                                                                                                            >  8:287726/1‑99 (MQ=255)
                                                                                      agagGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGcatcat                                                          >  8:6981/1‑148 (MQ=255)
                                                                                       gagGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCATCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATg                                                         >  6:45784/1‑148 (MQ=255)
                                                                                          gcggcAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGGTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAAc                                                     <  3:113874/149‑1 (MQ=255)
                                                                                            ggcAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAAc                                                     >  4:159465/1‑147 (MQ=255)
                                                                                                    aGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTTGTGCAGCATCATGAAACCAGCTATCGt                                           <  4:234951/149‑1 (MQ=255)
                                                                                                    aGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGt                                           >  7:144596/1‑149 (MQ=255)
                                                                                                    aGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAACCTGGTGCAGCATCATGAAACCAGCTATCGt                                           >  2:41663/1‑149 (MQ=255)
                                                                                                         ggtAAAACGTTGGTCGTCAAACGAACAGTGAAGTCACAAATGTACTACAAGCCGCACAGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTACCCTGCAGCCCCAGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTTACAt                                      >  8:84164/3‑149 (MQ=255)
                                                                                                            aaaaCGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCCGCTATCGTGACATTcc                                   >  5:253861/1‑148 (MQ=255)
                                                                                                             aaaCGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAACCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAg                                  <  1:41663/149‑1 (MQ=255)
                                                                                                                 gTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAAt                              >  8:236049/1‑149 (MQ=255)
                                                                                                                 gTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAAt                              >  6:91574/1‑149 (MQ=255)
                                                                                                                     gtcgtcATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTa                           <  3:159465/148‑1 (MQ=255)
                                                                                                                     gtcgtcATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTAt                          <  5:124895/149‑1 (MQ=255)
                                                                                                                     gtcggcATTCCCACTGGGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCGGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTa                           <  6:276467/148‑1 (MQ=255)
                                                                                                                      tcgtcATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACTCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTAtt                         <  5:378010/149‑1 (MQ=255)
                                                                                                                      tcgtcATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTAtt                         <  6:169310/149‑1 (MQ=255)
                                                                                                                       cgtcATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTAtt                         >  1:350552/1‑148 (MQ=255)
                                                                                                                           aTTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTgaga                     >  4:180697/1‑148 (MQ=255)
                                                                                                                            ttCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTCAAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGaaa                   >  2:251966/1‑149 (MQ=255)
                                                                                                                               gCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACAAGCTATCGTGACATTCAGCAATTTATTGAGAAacca                >  7:304850/1‑149 (MQ=255)
                                                                                                                                   tgtgCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGTTGCAGCATCATGACACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCgg            >  4:236672/1‑149 (MQ=255)
                                                                                                                                     tgCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCggg           <  8:21263/148‑1 (MQ=255)
                                                                                                                                      gCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGc          >  3:244571/1‑148 (MQ=255)
                                                                                                                                        cGTCACAAATGTACTACACGCCGAAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGATATCGTGACATTCAGCAATTTATTGAGAAACCACCTGGCAAg       <  2:324653/149‑1 (MQ=255)
                                                                                                                                         gTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGc      >  8:1646/1‑149 (MQ=255)
                                                                                                                                            acaAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCCGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGc   >  6:178828/1‑149 (MQ=255)
                                                                                                                                            acaAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCATTTTATTGAGAAACCACCGGGCAAGCTGc   <  4:86365/149‑1 (MQ=255)
                                                                                                                                            acaAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGc   <  5:45784/149‑1 (MQ=255)
                                                                                                                                            acaAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGc   >  5:238099/1‑149 (MQ=255)
                                                                                                                                            acaAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTACAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGc   <  2:220569/149‑1 (MQ=255)
                                                                                                                                              aaaTGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGCg  >  3:186243/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCGATACCTGGCTTTTATCGTAGCGGAGTGTCGCTGGAAGGCATTAACTACCTGGATGGCGGGATCAGTGATGCGATTCCGGTTAAAGAGGCGGCAAGGCAGGGCGCTAAAACGTTGGTCGTCATTCGCACTGTGCCGTCACAAATGTACTACACGCCGCAGTGGTTCAAACGCATGGAACGCTGGCTGGGTGACAGTAGCCTGCAGCCGCTGGTCAATCTGGTGCAGCATCATGAAACCAGCTATCGTGACATTCAGCAATTTATTGAGAAACCACCGGGCAAGCTGCG  >  NC_000913/4612903‑4613192

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