Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,844,146 T→C G61G (GGT→GGC adeD → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,1460TC99.2% 433.0 / NA 122G61G (GGT→GGCadeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (62/59);  minor base A (1/0);  total (63/59)
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

Only 100 of 122 total aligned reads displayed.

AGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGC  >  NC_000913/3844000‑3844293
                                                                                                                                                  |                                                                                                                                                    
aGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGtt                                                                                                                                                   <  2:310421/150‑1 (MQ=255)
aGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGtt                                                                                                                                                   >  2:1103896/1‑150 (MQ=255)
aGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTg                                                                                                                                                  >  1:125/1‑151 (MQ=255)
    gggCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGgcgc                                                                                                                                              >  4:329141/1‑151 (MQ=255)
    gggCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGgcgc                                                                                                                                              >  1:78774/1‑151 (MQ=255)
     ggCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGgcgc                                                                                                                                              >  1:24819/1‑150 (MQ=255)
         gAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAAt                                                                                                                                         <  3:296133/151‑1 (MQ=255)
         gAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAAt                                                                                                                                         >  2:17935/1‑151 (MQ=255)
             aCCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATacac                                                                                                                                     >  4:115788/1‑151 (MQ=255)
              ccAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACt                                                                                                                                    >  3:426549/1‑151 (MQ=255)
               cAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTg                                                                                                                                   <  1:889658/151‑1 (MQ=255)
                  gAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATg                                                                                                                                >  2:901529/1‑151 (MQ=255)
                   aaTTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGc                                                                                                                               >  2:809377/1‑151 (MQ=255)
                      tGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTcc                                                                                                                            >  3:130766/1‑151 (MQ=255)
                       gTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTcca                                                                                                                           <  2:118413/151‑2 (MQ=255)
                        ttAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCgg                                                                                                                          >  4:358640/1‑151 (MQ=255)
                             cGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTg                                                                                                                     >  1:895084/1‑151 (MQ=255)
                              gTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGc                                                                                                                    >  3:381601/1‑151 (MQ=255)
                                 tCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAg                                                                                                                  <  1:1281756/150‑1 (MQ=255)
                                  cccGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGc                                                                                                                 <  3:157245/150‑1 (MQ=255)
                                  cccGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTCACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCg                                                                                                                >  4:245079/1‑151 (MQ=255)
                                      tGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGAtt                                                                                                            >  2:1220317/1‑151 (MQ=255)
                                        gCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGa                                                                                                          <  1:931734/151‑1 (MQ=255)
                                           aCGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGc                                                                                                       <  2:1153500/151‑1 (MQ=255)
                                            cGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCt                                                                                                      <  1:473300/151‑1 (MQ=255)
                                            cGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCt                                                                                                      <  1:1212530/151‑1 (MQ=255)
                                             gCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTc                                                                                                     >  1:246624/1‑151 (MQ=255)
                                              cTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcg                                                                                                    >  2:453696/1‑151 (MQ=255)
                                                gTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTcg                                                                                                    <  1:1028300/149‑1 (MQ=255)
                                                 ttGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCgcgg                                                                                                 >  2:676672/1‑151 (MQ=255)
                                                   gCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCgcggcg                                                                                               >  4:282184/1‑151 (MQ=255)
                                                     cGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGgcgc                                                                                              >  1:1060336/1‑150 (MQ=255)
                                                     cGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCa                                                                                             <  2:659897/151‑1 (MQ=255)
                                                     cGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCa                                                                                             <  1:653341/151‑1 (MQ=255)
                                                         tatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcgg                                                                                         >  2:475100/1‑151 (MQ=255)
                                                          atattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggc                                                                                        >  2:673202/1‑151 (MQ=255)
                                                          atattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggc                                                                                        >  1:398068/1‑151 (MQ=255)
                                                            attattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcgg                                                                                      <  2:434801/151‑1 (MQ=255)
                                                              tattGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAcggcg                                                                                       >  3:297002/1‑148 (MQ=255)
                                                               attGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTg                                                                                    <  4:345742/150‑1 (MQ=255)
                                                                 tGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGcc                                                                                 >  3:246547/1‑150 (MQ=255)
                                                                  gATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCa                                                                                <  1:1262709/150‑1 (MQ=255)
                                                                  gATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAg                                                                               <  2:573611/151‑1 (MQ=255)
                                                                   aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAgg                                                                              <  1:856732/151‑1 (MQ=255)
                                                                   aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAgg                                                                              >  1:9185/1‑151 (MQ=255)
                                                                   aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAgg                                                                              >  3:308450/1‑151 (MQ=255)
                                                                     aaTGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGt                                                                            <  2:1056280/151‑1 (MQ=255)
                                                                       tGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGAGATTAAAGGACGTTACATTGCCGGAGTTGGCGCAGAATACACTGATGCCCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGcgcggg                                                                             >  2:1026624/1‑143 (MQ=255)
                                                                        gTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGttt                                                                          <  2:507780/150‑1 (MQ=255)
                                                                        gTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTa                                                                         <  4:2228/151‑1 (MQ=255)
                                                                          ctctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCTGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTAtt                                                                       >  1:430528/1‑151 (MQ=255)
                                                                          ctctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTAtt                                                                       >  2:512450/1‑151 (MQ=255)
                                                                            ctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTAtt                                                                       >  4:228351/1‑149 (MQ=255)
                                                                              aTTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAtt                                                                       >  1:47675/1‑147 (MQ=255)
                                                                              aTTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTg                                                                      <  1:719287/148‑1 (MQ=255)
                                                                              aTTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                   <  1:1192907/151‑1 (MQ=255)
                                                                              aTTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                   <  2:77925/151‑1 (MQ=255)
                                                                               ttCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGAt                                                                    >  1:393176/1‑149 (MQ=255)
                                                                                tctcGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGAttgt                                                                 >  1:377631/1‑148 (MQ=255)
                                                                                   cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCt                                                                 <  2:680505/148‑1 (MQ=255)
                                                                                    gACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTGGTGGTcccc                                                             >  2:313521/1‑151 (MQ=255)
                                                                                     aCCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCt                                                            <  2:38968/151‑1 (MQ=255)
                                                                                      ccTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGTAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGTTTTTTTATGCCCAcc                                                             >  1:992419/1‑149 (MQ=255)
                                                                                         gATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGGTGCTCACCTGCa                                                         >  1:249778/1‑150 (MQ=255)
                                                                                         gATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCat                                                        <  2:490479/151‑1 (MQ=255)
                                                                                         gATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCat                                                        <  1:290682/151‑1 (MQ=255)
                                                                                           tCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATTCTCACCTGCatat                                                      >  2:1010770/1‑151 (MQ=255)
                                                                                            cAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATAtt                                                     <  3:324844/151‑1 (MQ=255)
                                                                                             aaTGGCGGAGAAATTTCCGGCCCCATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCCGGGTTTTTTgtgtccaccctc                                                           >  3:299654/1‑134 (MQ=255)
                                                                                                ggCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGaa                                                  <  2:754237/150‑1 (MQ=255)
                                                                                                ggCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGaa                                                  <  2:1265661/150‑1 (MQ=255)
                                                                                                ggCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAAt                                                 <  2:199716/151‑1 (MQ=255)
                                                                                                  cGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATcc                                               <  2:898458/151‑1 (MQ=255)
                                                                                                   ggAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCa                                              <  2:5117/151‑1 (MQ=255)
                                                                                                    gagaAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGGCAGGGTTTATTGATGCTCACCTGCATATTGAATcc                                               >  3:369267/1‑149 (MQ=255)
                                                                                                       aaaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGTTTTTTGATGCCCACCTGCATATTGAATCCCGc                                            >  2:479205/1‑149 (MQ=255)
                                                                                                       aaaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGc                                            <  1:711316/149‑1 (MQ=255)
                                                                                                       aaaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCat                                          <  1:125991/151‑1 (MQ=255)
                                                                                                        aaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatg                                         <  1:957692/151‑1 (MQ=255)
                                                                                                        aaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTTATTGTCACCTTCCTATTGAATCCa                                              >  3:425178/1‑146 (MQ=255)
                                                                                                        aaTTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCat                                          >  1:82056/1‑150 (MQ=255)
                                                                                                          tttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTTATGCTCACCCGGATATTGAATTCAGCatg                                         >  2:277982/1‑149 (MQ=255)
                                                                                                             ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCAtgatga                                     <  2:872850/150‑1 (MQ=255)
                                                                                                             ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGAc                                    <  2:243231/151‑1 (MQ=255)
                                                                                                             ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGAc                                    >  1:328983/1‑151 (MQ=255)
                                                                                                             ccGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGAc                                    <  2:930101/151‑1 (MQ=255)
                                                                                                                gCCCAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGcc                                 <  2:1279945/151‑1 (MQ=255)
                                                                                                                  ccAATTGTGATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCgg                               <  1:901704/151‑1 (MQ=255)
                                                                                                                          gATTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACtttt                       <  2:389989/151‑1 (MQ=255)
                                                                                                                           aTTAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTg                      <  1:1258161/151‑1 (MQ=255)
                                                                                                                             tAAAGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGa                     <  1:5733/150‑1 (MQ=255)
                                                                                                                               aaGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaaa                   >  1:584378/1‑150 (MQ=255)
                                                                                                                                aGGACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGaaaac                 >  1:1105023/1‑149 (MQ=255)
                                                                                                                                 ggACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCg                <  1:762459/151‑1 (MQ=255)
                                                                                                                                  gACGTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCg                <  2:987216/150‑1 (MQ=255)
                                                                                                                                     gTTACATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTTAAACCGCTAc            >  2:367803/1‑151 (MQ=255)
                                                                                                                                         cATTGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTAccc          >  3:225608/1‑149 (MQ=255)
                                                                                                                                            tGCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGcc      >  2:644793/1‑150 (MQ=255)
                                                                                                                                             gCCGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcg     <  4:251108/150‑1 (MQ=255)
                                                                                                                                               cGGCGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATGCTCCCCCGCATATTGCATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgcg   >  3:103485/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGC  >  NC_000913/3844000‑3844293

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