Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,997 C→T N230N (AAC→AAT ygdH → UPF0717 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,9970CT95.1% 270.5 / 3.6 82N230N (AAC→AATygdHUPF0717 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (25/53);  total (27/55)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

TGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCGGCCAACAAA  >  NC_000913/2926851‑2927141
                                                                                                                                                  |                                                                                                                                                
tGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAAt                                                                                                                                              >  2:1033309/1‑151 (MQ=255)
  cACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCTCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTg                                                                                                                                            <  1:68866/151‑1 (MQ=255)
   aCCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTg                                                                                                                                            >  2:481366/1‑150 (MQ=255)
   aCCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGAGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCCCTGAACCGCCTAACCCGCTGGTCAATGAATTg                                                                                                                                            >  1:422451/1‑150 (MQ=255)
    ccGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGa                                                                                                                                           <  1:191462/150‑1 (MQ=255)
     cGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGatc                                                                                                                                         <  1:711711/151‑1 (MQ=255)
        cTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGatcatc                                                                                                                                      >  1:884230/1‑151 (MQ=255)
          gtgGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGAtcatca                                                                                                                                     <  2:945560/150‑1 (MQ=255)
          gtgGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATcatcat                                                                                                                                    <  2:549712/151‑1 (MQ=255)
          gtgGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATcatcat                                                                                                                                    <  2:587444/151‑1 (MQ=255)
             gTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGcc                                                                                                                                 <  2:751816/151‑1 (MQ=255)
               ccGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCa                                                                                                                                <  2:1125870/150‑1 (MQ=255)
                cGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGa                                                                                                                              <  1:101002/151‑1 (MQ=255)
                 ggtaGCGATGGAAGCGTCGATAAAAGGTGCTGCGGTCTGACACGCGCAGCAGCGTTTCAAACACACTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAg                                                                                                                               <  1:260912/146‑1 (MQ=255)
                  ggAGCGATGGAAGCGCCGATGAAAGGTGCTTCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGata                                                                                                                            <  2:685984/151‑1 (MQ=255)
                   gAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGatat                                                                                                                           >  2:16873/1‑151 (MQ=255)
                     gCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGATATc                                                                                                                          <  2:882017/150‑1 (MQ=255)
                          ggAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGaaaaa                                                                                                                    <  1:408248/151‑1 (MQ=255)
                          ggAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGaaaaa                                                                                                                    <  1:826310/151‑1 (MQ=255)
                           gAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAACGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAAc                                                                                                                   <  1:876290/151‑1 (MQ=255)
                            aaGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACg                                                                                                                  <  2:1256500/151‑1 (MQ=255)
                            aaGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACg                                                                                                                  <  2:69321/151‑1 (MQ=255)
                               cgcCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCt                                                                                                               >  2:795817/1‑151 (MQ=255)
                                  cGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTgg                                                                                                             >  1:784399/1‑150 (MQ=255)
                                    aTGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCTTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGaa                                                                                                           <  1:315745/150‑1 (MQ=255)
                                    aTGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCGGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGaa                                                                                                           <  1:1103014/150‑1 (MQ=255)
                                    aTGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGa                                                                                                            <  2:110224/149‑1 (MQ=255)
                                      gAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCg                                                                                                        <  2:252114/151‑1 (MQ=255)
                                        aaGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGtt                                                                                                      >  2:1211083/1‑151 (MQ=255)
                                          ggTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGGCAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTg                                                                                                    <  2:314304/151‑1 (MQ=255)
                                          ggTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTg                                                                                                    <  1:197709/151‑1 (MQ=255)
                                              ctgcGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCg                                                                                                >  1:1207905/1‑151 (MQ=255)
                                                 cGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTAt                                                                                             <  2:368953/151‑1 (MQ=255)
                                                          acGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACgg                                                                                    <  2:230870/151‑1 (MQ=255)
                                                                  cagcGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTa                                                                            <  2:602539/151‑1 (MQ=255)
                                                                   agcGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTAt                                                                           <  2:365586/151‑1 (MQ=255)
                                                                    gcGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATc                                                                          >  2:95025/1‑151 (MQ=255)
                                                                      gTTACAAAGACAGTCGTTTTATTGGTATGACTGAGCCGTCGATTAGCGCCGCGGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCt                                                                         <  2:401625/150‑1 (MQ=255)
                                                                      gTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCtt                                                                        >  2:1009675/1‑151 (MQ=255)
                                                                        tACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTc                                                                       <  2:1229136/150‑1 (MQ=255)
                                                                          cAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCt                                                                    <  1:487397/151‑1 (MQ=255)
                                                                           aaaGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTg                                                                   >  2:1099396/1‑151 (MQ=255)
                                                                           aaaGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTg                                                                   <  1:550826/151‑1 (MQ=255)
                                                                            aaGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTgg                                                                  >  2:383903/1‑151 (MQ=255)
                                                                             aGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGc                                                                 <  1:782239/151‑1 (MQ=255)
                                                                             aGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGc                                                                 >  2:235399/1‑151 (MQ=255)
                                                                               aCAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCgg                                                               >  1:665529/1‑151 (MQ=255)
                                                                                cAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCgg                                                               <  1:790692/150‑1 (MQ=255)
                                                                                 aGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGgtg                                                             >  1:39777/1‑151 (MQ=255)
                                                                                  gTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGgtgt                                                            <  1:68738/151‑1 (MQ=255)
                                                                                       tttATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGt                                                        <  2:292725/150‑1 (MQ=255)
                                                                                       tttATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGt                                                        >  2:535944/1‑150 (MQ=255)
                                                                                        ttATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTa                                                       <  2:113555/150‑1 (MQ=255)
                                                                                         tATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTAc                                                      <  1:1080616/150‑1 (MQ=255)
                                                                                         tATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTAc                                                      <  2:1007205/150‑1 (MQ=255)
                                                                                            tGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGc                                                   <  2:587045/150‑1 (MQ=255)
                                                                                            tGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGc                                                   <  1:654857/150‑1 (MQ=255)
                                                                                             ggTATGACAGAGCCGTTGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCa                                                  <  1:60131/150‑1 (MQ=255)
                                                                                              gTATGACAGAGCCGTCGATTATCGCCGTCGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACg                                                     <  1:779108/146‑1 (MQ=255)
                                                                                              gTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGc                                                   >  1:125499/1‑148 (MQ=255)
                                                                                                aTGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCagaag                                              <  2:1042945/151‑1 (MQ=255)
                                                                                                  gccAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCagaaga                                             <  2:327213/148‑1 (MQ=255)
                                                                                                    cAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGtt                                          >  2:387496/1‑151 (MQ=255)
                                                                                                    cAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGtt                                          <  2:791817/151‑1 (MQ=255)
                                                                                                        gCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGct                                       <  1:830876/150‑1 (MQ=255)
                                                                                                          cGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGctct                                     >  1:551946/1‑150 (MQ=255)
                                                                                                          cGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTa                                    <  2:449339/151‑1 (MQ=255)
                                                                                                           gTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTa                                    <  1:686314/150‑1 (MQ=255)
                                                                                                            tCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTAt                                   <  1:1046946/150‑1 (MQ=255)
                                                                                                             cGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTAttt                                 <  1:677780/151‑1 (MQ=255)
                                                                                                                    cgccgcTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTggg                           <  1:1186618/150‑1 (MQ=255)
                                                                                                                     gccgcTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGa                          >  1:226791/1‑150 (MQ=255)
                                                                                                                        gcTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAAttt                      <  2:879222/151‑1 (MQ=255)
                                                                                                                        gcTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAAttt                      >  2:1156398/1‑151 (MQ=255)
                                                                                                                         cTGAACCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAAtttt                     >  1:437563/1‑151 (MQ=255)
                                                                                                                            aaCCGCCTAACCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTaa                   <  1:455302/150‑1 (MQ=255)
                                                                                                                                    aaCCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAAccc            <  1:541824/149‑1 (MQ=255)
                                                                                                                                    aaCCCGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCgg          >  2:869457/1‑151 (MQ=255)
                                                                                                                                        cGCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCGGcc        >  1:661933/1‑149 (MQ=255)
                                                                                                                                         gCTGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCGGCcaac     >  2:759158/1‑151 (MQ=255)
                                                                                                                                           tGGTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCGGCcaaca    >  2:440998/1‑150 (MQ=255)
                                                                                                                                            ggTCAATGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCGGCCAACaaa  <  1:242857/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATCATGCCAGATATCGAAAAACGTCTGGAAGCGTTTGTCCGTATCGCTCACGGTATCATTATCTTCCCTGGCGGTGTGGGTACGGCAGAAGAGTTGCTCTATTTGCTGGGAATTTTAATGAACCCGGCCAACAAA  >  NC_000913/2926851‑2927141

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: