Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 8.8% intergenic (‑82/+36) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8550CA 8.8% 653.8 / 37.9 193intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (88/87);  new base A (7/10);  total (96/97)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 201 total aligned reads displayed.

GCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTA  >  NC_000913/4366589‑4367088
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
                       aaTGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCCACCCCAATCCGACTGCAAAACCAATACCTATTACCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTAAAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTGGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTg                                                                                                                                                                                                             <  4:431084/273‑1 (MQ=255)
                                                                                                                      cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                                                                                                           >  7:852335/1‑149 (MQ=255)
                                                                                                                          tCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGAc                                                                                                                                                                                                                                       <  5:81966/149‑1 (MQ=255)
                                                                                                                             ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACg                                                                                                                                                                                                                                      <  4:153792/147‑1 (MQ=255)
                                                                                                                             ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                                                                                                                    >  7:724420/1‑149 (MQ=255)
                                                                                                                              cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                                                                                                                   >  7:913132/1‑149 (MQ=255)
                                                                                                                              cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                                                                                                                   >  6:1380662/1‑149 (MQ=255)
                                                                                                                              cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGc                                                                                                                                                                                                                                   <  6:1265252/149‑1 (MQ=255)
                                                                                                                               ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGcc                                                                                                                                                                                                                                  <  6:588464/149‑1 (MQ=255)
                                                                                                                                cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCt                                                                                                                                                                                                                                 >  5:637051/1‑149 (MQ=255)
                                                                                                                                cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCt                                                                                                                                                                                                                                 <  6:791602/149‑1 (MQ=255)
                                                                                                                                  aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCttt                                                                                                                                                                                                                               >  8:1688776/1‑149 (MQ=255)
                                                                                                                                  aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                                                                                                        <  2:25375/188‑1 (MQ=255)
                                                                                                                                         gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                                                                                                        <  5:916919/149‑1 (MQ=255)
                                                                                                                                           cgctgtagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGt                                                                                                                                                                                                                      <  8:352632/143‑1 (MQ=255)
                                                                                                                                            cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                     >  7:540323/1‑149 (MQ=255)
                                                                                                                                            cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                     <  5:1281141/149‑1 (MQ=255)
                                                                                                                                             aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                                                                                                                    >  5:1024808/1‑149 (MQ=255)
                                                                                                                                              agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACg                                                                                     <  2:565031/275‑1 (MQ=255)
                                                                                                                                                   aTCGCCAGCAAAACTATGATGAATTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTGTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGt                                                                            <  1:166288/279‑1 (MQ=255)
                                                                                                                                                      gCCAACAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                                                                                                           <  7:585561/149‑1 (MQ=255)
                                                                                                                                                        cAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACg                                                                                                                                                                                                         >  5:480309/1‑149 (MQ=255)
                                                                                                                                                          gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                                                                                                          <  6:907244/130‑1 (MQ=255)
                                                                                                                                                          gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                                                                                                       <  7:407746/149‑1 (MQ=255)
                                                                                                                                                                    gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                                                                                                             >  7:229283/1‑149 (MQ=255)
                                                                                                                                                                    gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGACTTTTTTATTTGTACTACCCTGTACGTTAACTGTTCGCTGTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                                                                                                              >  2:475716/1‑196 (MQ=255)
                                                                                                                                                                     atgaGTTCTACAACTAGCATCTTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                                                                                                               <  8:1174250/146‑1 (MQ=255)
                                                                                                                                                                     atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTGAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                                                                                                               >  7:1174250/1‑146 (MQ=255)
                                                                                                                                                                       tttgtttACAACTAGTTTATTGTCTTTCATTGTNTTTTAAACATTTGTTATTAGATTTTNGTGTTTAATTTACTGCTCAAAAGAAAAAAGGTACGTCATCTGACGTGCCTGTTTTATTTGTACTACCCGGTACGATGTCTNTTCGCTTTCATCAGTATAGCGTTTTGCTTTGGTAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGc                                                        <  3:174203/273‑1 (MQ=255)
                                                                                                                                                                          ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                                                                                                         >  2:764085/1‑115 (MQ=255)
                                                                                                                                                                          ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAccc                                                                                                                                                                                                               <  3:37032/125‑1 (MQ=255)
                                                                                                                                                                          ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                                                                                                       >  5:1311811/1‑149 (MQ=255)
                                                                                                                                                                            cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                                                                                                     <  5:1380662/149‑1 (MQ=255)
                                                                                                                                                                            cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                                                                                                     <  8:724420/149‑1 (MQ=255)
                                                                                                                                                                             tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                                                                                                                    >  8:1542129/1‑149 (MQ=255)
                                                                                                                                                                             tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                                                                                                                    <  5:531507/149‑1 (MQ=255)
                                                                                                                                                                             tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                                                                                                                    >  7:144845/1‑149 (MQ=255)
                                                                                                                                                                              aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                     <  4:723046/115‑1 (MQ=255)
                                                                                                                                                                              aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGAc                                                                                   <  2:202782/245‑1 (MQ=255)
                                                                                                                                                                                aaCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtat                                                                                                                                                                                 >  6:1190574/1‑149 (MQ=255)
                                                                                                                                                                                    aGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                                                                                                         >  2:740532/1‑105 (MQ=255)
                                                                                                                                                                                     gCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGa                                                                                                                                                                                                                                        <  6:775437/89‑1 (MQ=255)
                                                                                                                                                                                     gCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGt                                                                                                                                                                            >  6:698498/1‑149 (MQ=255)
                                                                                                                                                                                     gCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGt                                                                                                                                                                            >  5:1429077/1‑149 (MQ=255)
                                                                                                                                                                                       atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGttt                                                                                                                                                                          <  5:222367/149‑1 (MQ=255)
                                                                                                                                                                                        tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGtttt                                                                                                                                                                         >  5:1687044/1‑149 (MQ=255)
                                                                                                                                                                                          ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                                                                                                                       >  5:791389/1‑149 (MQ=255)
                                                                                                                                                                                          ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                                                                                                                       >  7:730429/1‑149 (MQ=255)
                                                                                                                                                                                           tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                                                                                                                      <  6:722375/149‑1 (MQ=255)
                                                                                                                                                                                             gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                                                                                                                    <  6:1311811/149‑1 (MQ=255)
                                                                                                                                                                                             gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                                                                                                                    <  5:698498/149‑1 (MQ=255)
                                                                                                                                                                                             gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                                                                                                                    >  8:1216601/1‑149 (MQ=255)
                                                                                                                                                                                                    ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                                                                                                                    >  6:493574/1‑94 (MQ=255)
                                                                                                                                                                                                     tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCg                                                                                                                                                            >  6:886750/1‑149 (MQ=255)
                                                                                                                                                                                                       ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                                                                                                                         <  1:256614/166‑1 (MQ=255)
                                                                                                                                                                                                        tttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGt                                                                                                                                                         <  6:1554816/149‑1 (MQ=255)
                                                                                                                                                                                                            tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                                                                                                           >  5:15388/1‑63 (MQ=255)
                                                                                                                                                                                                            tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                                                                                                           <  6:15388/63‑1 (MQ=255)
                                                                                                                                                                                                              aCAAGTTGAGATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTGATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTAGTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCa                                                                                                                                                   >  8:489206/1‑149 (MQ=255)
                                                                                                                                                                                                                aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                                                                                                       <  2:127082/111‑1 (MQ=255)
                                                                                                                                                                                                                aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGa                                                                                                                                                 <  8:523695/149‑1 (MQ=255)
                                                                                                                                                                                                                  gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                                                                                                     <  7:845876/143‑1 (MQ=255)
                                                                                                                                                                                                                  gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACATCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                                                                                                     >  8:845876/1‑143 (MQ=255)
                                                                                                                                                                                                                   ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGTTGCATCAGATTc                                                                                                                                              <  8:852335/149‑1 (MQ=255)
                                                                                                                                                                                                                   ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                                                                                                              >  7:1310478/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                                                                                                              <  8:540323/149‑1 (MQ=255)
                                                                                                                                                                                                                     tgtatttgtgtGTAATTATTAAGATAGTACTCTGAAGAAAAAAGGTACGGCATCGGATGGTGCTGTTGTATTTGGACAACACAGGAGGATTACGGGTCGCGTTAATCTGTATATGTTATTTGTGTGTCAGCCGGGTGAATTGGACGCATTACAGATAGTCTAACGTCCAATTCCGCTTCAGTCAACAGACCGCGTTCCAGCACGACTTCACGTTCAACCTTACGGGATGCGGCAGAGAGTTTACCCTCGATGCCACCGTTCAGGTGACGGATGag              <  4:34825/265‑2 (MQ=255)
                                                                                                                                                                                                                     gATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                                                                                                          <  4:758630/167‑1 (MQ=255)
                                                                                                                                                                                                                      atatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                                                                                                           <  8:1310478/149‑1 (MQ=255)
                                                                                                                                                                                                                       tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                                                                                                      >  5:937566/1‑105 (MQ=255)
                                                                                                                                                                                                                       tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                                                                                                      <  6:937566/105‑1 (MQ=255)
                                                                                                                                                                                                                       tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                                                                                                          <  6:827281/149‑1 (MQ=255)
                                                                                                                                                                                                                         ttAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                                                                                                     <  6:286673/104‑1 (MQ=255)
                                                                                                                                                                                                                         ttAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACg                                                                                                                                        >  7:1505975/1‑149 (MQ=255)
                                                                                                                                                                                                                            gATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACg                                                                                                                                        <  8:1512152/146‑1 (MQ=255)
                                                                                                                                                                                                                             aTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgaga                                                                                                                                    >  5:751181/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTCCTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTTTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATACCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTTAACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTTTGTTGACCGATGCAGGGACTCAGgta  >  4:496814/1‑275 (MQ=255)
                                                                                                                                                                                                                                 ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                                                                                                                >  6:664502/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                                                                                                                >  8:663881/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aTTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAatat                                                                                                                              <  5:886750/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ttttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGa                  <  4:230455/255‑1 (MQ=255)
                                                                                                                                                                                                                                       ttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                                                                                                          <  7:478969/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aaGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCaa                                                                                                                        <  8:913132/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ttACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                                                                                                     >  5:588468/1‑149 (MQ=255)
                                                                                                                                                                                                                                            ttACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                                                                                                     >  7:1052273/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTcc                                                                                                                   >  5:365463/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTcc                                                                                                                   <  5:1190574/149‑1 (MQ=255)
                                                                                                                                                                                                                                                tGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                                                                                                                 <  6:1661926/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     acaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                                                                                                             <  7:1630601/148‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                                                                                                           >  6:609440/5‑149 (MQ=255)
                                                                                                                                                                                                                                                      aaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                                                                                                           >  7:690904/5‑149 (MQ=255)
                                                                                                                                                                                                                                                       aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                                                                          >  6:1450371/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCaa                                                                                                        >  5:1130184/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                                                                                                    <  7:1542129/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                                                                                                                  >  7:1067175/1‑135 (MQ=255)
                                                                                                                                                                                                                                                               aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                                                                  >  8:1118590/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aGGCACGTCAGATGACGTGCCTTTTTTATCTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                                                                          >  7:943828/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                 gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGGGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc                                                                                                >  7:1244516/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                                                                         <  5:494448/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa                                                                                       <  8:144845/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
GCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTA  >  NC_000913/4366589‑4367088

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