Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,621 T→G 8.5% Q203P (CAA→CCA)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,6210TG 8.5% 451.8 / 31.1 130Q203P (CAA→CCA) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base T (61/58);  new base G (5/6);  total (66/64)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 130 total aligned reads displayed.

CTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCT  >  NC_000913/4080475‑4080761
                                                                                                                                                  |                                                                                                                                            
cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                             >  1:1039389/1‑148 (MQ=255)
cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                             <  2:1039389/148‑1 (MQ=255)
     aTCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                        <  1:1033112/148‑1 (MQ=255)
        gcgGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGc                                                                                                                                     <  4:306999/148‑1 (MQ=255)
         cgGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCt                                                                                                                                    <  3:933684/148‑1 (MQ=255)
                ttCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACAcc                                                                                                                             >  3:488410/1‑148 (MQ=255)
                       tCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAg                                                                                                                      >  4:919727/1‑148 (MQ=255)
                       tCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAg                                                                                                                      >  3:447630/1‑148 (MQ=255)
                        ccAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                        >  3:931166/1‑129 (MQ=255)
                        ccAGCGACCAGTAATGCAGTTTGCCGCACTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCa                                                                                                                                        <  4:931166/129‑1 (MQ=255)
                         cAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                    <  4:447630/148‑1 (MQ=255)
                          aGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                   >  3:1083421/1‑148 (MQ=255)
                            cGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                 <  3:919727/148‑1 (MQ=255)
                               ccAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACaa                                                                                                              <  4:1083421/148‑1 (MQ=255)
                                  gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCCCACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAAcc                                                                                                           <  3:322395/148‑1 (MQ=255)
                                  gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAAcc                                                                                                           >  4:322395/1‑148 (MQ=255)
                                   tAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCg                                                                                                          <  4:488410/148‑1 (MQ=255)
                                                  ctctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                <  2:344369/143‑1 (MQ=255)
                                                   tctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                 <  1:579617/125‑1 (MQ=255)
                                                                   gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCa                                                                                                                                        <  1:863874/86‑1 (MQ=255)
                                                                   gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCa                                                                                                                                        >  2:863874/1‑86 (MQ=255)
                                                                   gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                    >  3:791426/1‑106 (MQ=255)
                                                                   gcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                    <  4:791426/106‑1 (MQ=255)
                                                                       tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                <  2:1125490/122‑1 (MQ=255)
                                                                       tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                >  3:555038/1‑122 (MQ=255)
                                                                       tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                <  4:555038/122‑1 (MQ=255)
                                                                       tttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCggg                                                                      >  1:1102641/1‑148 (MQ=255)
                                                                         tACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCactca                                                                                                                                    >  3:628882/1‑80 (MQ=39)
                                                                         tACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCactca                                                                                                                                    <  4:628882/84‑5 (MQ=39)
                                                                           cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGc                                                                                                                          <  3:818259/92‑1 (MQ=255)
                                                                           cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                          >  2:441416/1‑124 (MQ=255)
                                                                           cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                          <  1:441416/124‑1 (MQ=255)
                                                                               cACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCg                                                              <  3:1094839/148‑1 (MQ=255)
                                                                                 cTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAat                                                                                                    <  4:1019124/108‑1 (MQ=255)
                                                                                  tACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTAt                                                                                <  4:653404/127‑1 (MQ=255)
                                                                                       gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      <  2:1102641/148‑1 (MQ=255)
                                                                                       gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      >  2:746794/1‑148 (MQ=255)
                                                                                           ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                           <  2:899230/139‑1 (MQ=255)
                                                                                           ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                           >  1:899230/1‑139 (MQ=255)
                                                                                            tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATa                                                                                                                       >  3:1113446/1‑78 (MQ=255)
                                                                                            tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATa                                                                                                                       <  4:1113446/78‑1 (MQ=255)
                                                                                            tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                         >  3:97326/1‑92 (MQ=255)
                                                                                            tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                         <  4:97326/92‑1 (MQ=255)
                                                                                            tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGAGGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTgg                                                 >  3:982713/1‑148 (MQ=255)
                                                                                             cAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTc                                                                              <  1:250293/118‑1 (MQ=255)
                                                                                               gTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACGCCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      <  4:51027/140‑1 (MQ=255)
                                                                                               gTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACGCCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      >  3:51027/1‑140 (MQ=255)
                                                                                                     cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                  >  1:615504/1‑138 (MQ=255)
                                                                                                     cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                        >  1:351780/1‑148 (MQ=255)
                                                                                                     cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                        >  3:1073007/1‑148 (MQ=255)
                                                                                                     cacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATAAGATTg                                                  <  2:615504/138‑1 (MQ=255)
                                                                                                       caGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                          <  3:645512/128‑1 (MQ=255)
                                                                                                       caGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGcc                                      >  3:556558/1‑148 (MQ=255)
                                                                                                         ggTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                         <  3:186389/79‑1 (MQ=255)
                                                                                                         ggTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                         >  4:186389/1‑79 (MQ=255)
                                                                                                           ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCt                                                                                                                                    <  4:824452/50‑1 (MQ=255)
                                                                                                           ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                  >  2:653527/1‑148 (MQ=255)
                                                                                                           ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                  >  3:680249/1‑148 (MQ=255)
                                                                                                           ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                  <  4:680249/148‑1 (MQ=255)
                                                                                                            tgcagtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGt                                                                                            >  1:1085701/1‑89 (MQ=255)
                                                                                                                gtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGc                             <  4:1073007/148‑1 (MQ=255)
                                                                                                                 tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGt                                                                                  >  4:622019/1‑94 (MQ=255)
                                                                                                                 tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                           >  1:988421/1‑117 (MQ=255)
                                                                                                                 tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                           <  2:988421/117‑1 (MQ=255)
                                                                                                                  gcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                          >  1:68492/1‑133 (MQ=255)
                                                                                                                       tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGa                                                                                               <  4:42978/75‑1 (MQ=255)
                                                                                                                       tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGa                                                                                               >  3:42978/1‑75 (MQ=255)
                                                                                                                       tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      <  1:130640/116‑1 (MQ=255)
                                                                                                                       tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      >  2:130640/1‑116 (MQ=255)
                                                                                                                       tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGAc                          <  3:1074852/144‑1 (MQ=255)
                                                                                                                       tAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCAGCGCCCAGTAGAGCGAc                          >  4:1074852/1‑144 (MQ=255)
                                                                                                                            tAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGcc                 >  3:278392/1‑148 (MQ=255)
                                                                                                                             aGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGccc                >  2:450367/1‑148 (MQ=255)
                                                                                                                               ccccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAatat                                                                                                  <  3:184722/64‑1 (MQ=255)
                                                                                                                               ccccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAatat                                                                                                  >  4:184722/1‑64 (MQ=255)
                                                                                                                                 ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                >  3:565244/1‑64 (MQ=255)
                                                                                                                                 ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                <  4:565244/64‑1 (MQ=255)
                                                                                                                                 ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      <  3:634680/106‑1 (MQ=255)
                                                                                                                                 ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAg                                                      >  4:634680/1‑106 (MQ=255)
                                                                                                                                   tGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                        >  1:58359/1‑118 (MQ=255)
                                                                                                                                   tGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                        <  2:58359/118‑1 (MQ=255)
                                                                                                                                   tGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATaa          <  4:556558/148‑1 (MQ=255)
                                                                                                                                      gtggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACg       <  2:351780/148‑1 (MQ=255)
                                                                                                                                      atggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                          >  2:648085/2‑113 (MQ=255)
                                                                                                                                      atggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                          <  1:648085/112‑1 (MQ=255)
                                                                                                                                      atggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                          >  4:310124/2‑113 (MQ=255)
                                                                                                                                      atggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                          <  3:310124/112‑1 (MQ=255)
                                                                                                                                       tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTag                                >  4:9790/1‑122 (MQ=255)
                                                                                                                                       tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTag                                <  3:9790/122‑1 (MQ=255)
                                                                                                                                       tggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAg                        <  3:1003364/130‑1 (MQ=255)
                                                                                                                                       tggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCg                                                              >  2:544661/1‑92 (MQ=255)
                                                                                                                                        ggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAg                                                                            <  1:769762/77‑1 (MQ=255)
                                                                                                                                        ggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAg                                                                            >  2:769762/1‑77 (MQ=255)
                                                                                                                                           gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAt                                                    >  2:198506/1‑98 (MQ=255)
                                                                                                                                           gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAt                                                    <  1:198506/98‑1 (MQ=255)
                                                                                                                                           gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                               >  2:377903/1‑119 (MQ=255)
                                                                                                                                           gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                               <  1:377903/119‑1 (MQ=255)
                                                                                                                                           gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTGGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTg                                          >  4:961541/1‑108 (MQ=255)
                                                                                                                                              aaTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                           >  4:803941/1‑88 (MQ=255)
                                                                                                                                              aaTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                           <  3:803941/88‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCT  >  NC_000913/4080475‑4080761

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: