Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,555 A→T 6.5% intergenic (+44/‑176) metH → / → yjbB homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,5550AT 6.5% 488.8 / 18.3 138intergenic (+44/‑176)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (58/71);  new base T (3/6);  total (61/77)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 140 total aligned reads displayed.

AAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAA  >  NC_000913/4227408‑4227702
                                                                                                                                                   |                                                                                                                                                   
aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                                                                    <  3:304403/149‑1 (MQ=255)
 aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                                                                    >  1:770331/1‑148 (MQ=255)
  aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                                                                                  <  3:356598/149‑1 (MQ=255)
       gcgcGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAgcac                                                                                                                                              <  1:54735/148‑4 (MQ=255)
         gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                                                            >  1:123079/1‑148 (MQ=255)
         gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                                                            <  1:768882/148‑1 (MQ=255)
         gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCt                                                                                                                                           <  3:726361/149‑1 (MQ=255)
             tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                                                                        <  2:153927/148‑1 (MQ=255)
              cAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                                                                                   <  2:349229/136‑1 (MQ=255)
              cAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                                                                      <  4:697431/149‑1 (MQ=255)
                 gTTGAAGATTATGCCCGCCGTAAAGGTATGAGAGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                                                                    <  1:108426/148‑1 (MQ=255)
                   tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGACACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                                                                   <  2:771318/147‑1 (MQ=255)
                    gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGt                                                                                                                                                <  3:501565/133‑1 (MQ=255)
                     aaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                                                                                  >  4:23961/1‑130 (MQ=255)
                     aaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                                                                                  <  3:23961/130‑1 (MQ=255)
                      aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                                                            >  1:745203/1‑135 (MQ=255)
                        aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                                                                             >  1:288827/1‑148 (MQ=255)
                        aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                                                                             <  1:616308/148‑1 (MQ=255)
                         ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatta                                                                                                                               >  1:80938/1‑145 (MQ=255)
                          tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                                                                          <  4:381095/149‑1 (MQ=255)
                                cgccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                     >  1:96261/1‑148 (MQ=255)
                                 gccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACCACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                                                                   <  4:537936/149‑1 (MQ=255)
                                        aGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                                                             <  1:569915/148‑1 (MQ=255)
                                         ggTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                                                                           <  4:153096/149‑1 (MQ=255)
                                            aTGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                                                                     >  2:718657/1‑116 (MQ=255)
                                                 cGTTACCGAAGTTGAGCGCTGGCTGGCACCGCATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                                                                                    >  2:766079/1‑148 (MQ=255)
                                                   ttACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCt                                                                                                  >  1:60506/1‑148 (MQ=255)
                                                        gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACtt                                                                                                                  >  1:672705/1‑127 (MQ=255)
                                                        gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACtt                                                                                                                  <  2:672705/127‑1 (MQ=255)
                                                          aGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                                                                          <  3:257876/149‑1 (MQ=255)
                                                            ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                                                                        <  4:672753/149‑1 (MQ=255)
                                                             tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                                                                        >  1:642235/1‑148 (MQ=255)
                                                              gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                                                                      <  3:807771/149‑1 (MQ=255)
                                                               aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                                                                                   >  2:207700/1‑87 (MQ=255)
                                                               aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                                                                                   <  1:207700/87‑1 (MQ=255)
                                                                gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                                                                     <  2:770331/148‑1 (MQ=255)
                                                                gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTACTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                                                                                    >  3:523738/1‑149 (MQ=255)
                                                                  gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACCGGAAACAGGCTGGAATAATCTTAGCCg                                                                                  <  4:452659/149‑1 (MQ=255)
                                                                     ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                                                                               <  3:201707/149‑1 (MQ=255)
                                                                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                                                                     >  4:797568/1‑142 (MQ=255)
                                                                        tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                                                                            <  4:3015/149‑1 (MQ=255)
                                                                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                                                                        >  4:643208/1‑117 (MQ=255)
                                                                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                                                                        <  3:643208/117‑1 (MQ=255)
                                                                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                                                                          >  4:798718/1‑131 (MQ=255)
                                                                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                                                                          <  3:798718/131‑1 (MQ=255)
                                                                              cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                                                                       >  1:399851/1‑148 (MQ=255)
                                                                               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                                                                      <  1:751890/148‑1 (MQ=255)
                                                                                aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                                                                       <  1:815962/130‑1 (MQ=255)
                                                                                aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                                                                       >  2:815962/1‑130 (MQ=255)
                                                                                   cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                                                                  <  1:251674/148‑1 (MQ=255)
                                                                                       ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                                                             <  4:872389/149‑1 (MQ=255)
                                                                                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                                                             <  1:314716/148‑1 (MQ=255)
                                                                                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                                                             >  2:195339/1‑148 (MQ=255)
                                                                                               gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                                                                       >  3:194117/1‑65 (MQ=25)
                                                                                               gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                                                                       <  4:194117/67‑3 (MQ=25)
                                                                                                 ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                                                                       >  4:248834/1‑129 (MQ=255)
                                                                                                 ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGa                                                   <  4:379396/149‑1 (MQ=255)
                                                                                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAt                                                  <  4:424637/149‑1 (MQ=255)
                                                                                                     tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                                                             <  2:681662/135‑1 (MQ=255)
                                                                                                     tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc                                               >  3:232810/1‑149 (MQ=255)
                                                                                                      gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTTAAAAGGCACAGTAAAGATGcc                                               <  1:423830/148‑1 (MQ=255)
                                                                                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAg                                             >  3:235349/1‑149 (MQ=255)
                                                                                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAg                                             >  2:669883/1‑148 (MQ=255)
                                                                                                         tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTc                                           >  4:677070/1‑149 (MQ=255)
                                                                                                            caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                         >  2:149143/1‑148 (MQ=255)
                                                                                                            caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGc                                        <  3:2724/149‑1 (MQ=255)
                                                                                                              aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCg                                       >  1:744815/1‑148 (MQ=255)
                                                                                                               aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                                                                        >  3:407651/1‑114 (MQ=255)
                                                                                                               aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                                                                        <  4:407651/114‑1 (MQ=255)
                                                                                                                   gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGttt                                  <  2:641234/148‑1 (MQ=255)
                                                                                                                       cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTAGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGt                             >  4:536592/1‑149 (MQ=255)
                                                                                                                          ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                                                                                    <  4:303946/75‑1 (MQ=255)
                                                                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTg                       <  3:756221/149‑1 (MQ=255)
                                                                                                                              ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                      >  4:757703/1‑149 (MQ=255)
                                                                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                                                             <  2:264534/109‑1 (MQ=255)
                                                                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                      >  1:689788/1‑148 (MQ=255)
                                                                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGt                     >  3:403848/1‑149 (MQ=255)
                                                                                                                                 caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTCTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaa                                                      >  1:889489/1‑114 (MQ=255)
                                                                                                                                 caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTCTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaa                                                      <  2:889489/114‑1 (MQ=255)
                                                                                                                                  aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt                  >  3:743343/1‑149 (MQ=255)
                                                                                                                                    caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAtt                 <  2:399851/148‑1 (MQ=255)
                                                                                                                                    caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAtt                 <  2:60506/148‑1 (MQ=255)
                                                                                                                                     aaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCa               >  4:309722/1‑149 (MQ=255)
                                                                                                                                      aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCa               <  1:880154/148‑1 (MQ=255)
                                                                                                                                       acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCATCGTGTTTCCGGTCCTTTGGTCATTCAAt             <  3:777677/149‑1 (MQ=255)
                                                                                                                                       acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCaa              <  1:688393/148‑1 (MQ=255)
                                                                                                                                        caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAt             <  2:893347/148‑1 (MQ=255)
                                                                                                                                        caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTACGGTCCTTTGGTCATTCAAt             <  2:813928/148‑1 (MQ=255)
                                                                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCTAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg           >  1:295112/1‑148 (MQ=255)
                                                                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc          <  4:231564/149‑1 (MQ=255)
                                                                                                                                          gggCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg           >  1:887372/1‑148 (MQ=255)
                                                                                                                                          gggCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg           <  2:887372/148‑1 (MQ=255)
                                                                                                                                           ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCt         >  3:340766/1‑149 (MQ=255)
                                                                                                                                            gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCt         <  1:149143/148‑1 (MQ=255)
                                                                                                                                              gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTc       <  1:669883/148‑1 (MQ=255)
                                                                                                                                              gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTc       <  2:69499/148‑1 (MQ=255)
                                                                                                                                              gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGACTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCa      >  4:712957/1‑149 (MQ=255)
                                                                                                                                                tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACTGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAg     <  1:587439/148‑1 (MQ=255)
                                                                                                                                                tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAgg    >  4:836699/1‑149 (MQ=255)
                                                                                                                                                  aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGaa  <  3:536592/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAGGAA  >  NC_000913/4227408‑4227702

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