Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,974 C→T 10.1% intergenic (+448/+92) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9740CT10.1% 122.0 / 5.6 55intergenic (+448/+92)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (2/3);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

GCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGC  >  NC_000913/1286736‑1287207
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                                                                                                                           <  1:632623/239‑1 (MQ=255)
    cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAGGGATTACTCATCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCCGCGCTCTCTCGCTGACTCTCGAGccag                                                                                                                                                                                                                                        >  1:1054183/1‑236 (MQ=255)
    cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                                                                                                                        >  2:95395/1‑238 (MQ=255)
                          tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTATTCGTCGCCTCGCGCCTCACCCCTCGGGCCGTTGCCTGTTGCAACGTTCTCTCGCTGTCGCCCGAGTCGAACCTTAGTCGAAGGTTCTCac                                                                                                                                                                                                                   >  2:307817/1‑236 (MQ=14)
                                   cAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTAGTCGAAGCTTCGCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                                                                                                                         <  2:152390/238‑1 (MQ=255)
                                       gatcgctctctgccgatcagaccGACCCTGGGTGGAAGTTTAGGATCGTTCCCCGGTGGGGCAGAGTACTTGATTGCGGATTAGTTTGAGATTTCGGGGATTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                                                                                                                    <  1:144260/217‑1 (MQ=11)
                                         tgtctcGTTGGCGGTGAAGTCGATGCGTGGTGATAGTTTATCATCGTTCCCCGTTTGGGCAGAATATTTGATTGGGGAGTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                                  <  1:989915/237‑1 (MQ=14)
                                         tctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                                  <  2:440462/239‑1 (MQ=255)
                                         tctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                                  <  2:676393/239‑1 (MQ=255)
                                             ccGCTGGCGCACAATCCGAGCCTGTGTCGACAGTTGTCAGCCTTCCCCGCTTGGGCAGAATGTTTGATTGCGGATTCGTTGGAAAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCCCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCgatgagtgcaaa                                                                                                                                                                                               <  1:307817/237‑12 (MQ=2)
                                                         aaGTCGAACCTTGGTCGAAGCTTCTCATCCATCCTCGATTGGGCAGAATGTCAGATGGGTAATTGGGTTAAGAATTCCGTGATTTAGGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                                                                                                                                  <  1:787716/239‑1 (MQ=14)
                                                                           agtggctcaccgtgcaccGTTTGGGCAGAATATTGGATTGCGGCTGAGTGTGAGGATTGCGGGGCTTGTGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                                                                                                                 <  1:668402/222‑1 (MQ=17)
                                                                            gCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATGACTCTTCGCTGCGTGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCTGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                                                                               >  1:168169/1‑239 (MQ=255)
                                                                                ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCGCGGCGCGCTTCGCCGTTCGGGTCGTTGCTTGCGGCTGCCCTCTCTAGCTGCCGCTCGCGCCCGACCTTAGTCCTTGCTTCTCCCCCTTCCCCCCCTCGGCGCCACTTTTCCTTGCGGATTCCCTTGTCACTTCAgggc                                                                                                                                                             >  1:282892/1‑236 (MQ=9)
                                                                                   aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATGCGTTTGAGAATTCCGGGGTTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                                        <  2:771011/239‑1 (MQ=255)
                                                                                      cgtcaacGCTGGGGCAGAAGATTTGATTGGGGAGTCGTTTGAGAATTCGGGGGCTTTAGAAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                                                                                                     <  1:1023349/233‑1 (MQ=14)
                                                                                        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGTTTCGCTTGAATGTTCAGGGACTTTTGAAAg                                                                                                                                                   >  2:601896/1‑239 (MQ=255)
                                                                                        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCATCGCTCTCTCGCTGCCCCTCGCGCCGAACCTTCGTAGAAGCATCTCATCCTTCCCCTGCTGGGCAGGGTATTTTTTTTCGGATTCGTTTGCCTGTTCAGGGGCttttttgaa                                                                                                                                                   >  1:980660/1‑234 (MQ=11)
                                                                                          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTCCGGATTAGCTGTAGAGTTAAGGGACTTTTGGAAGTg                                                                                                                                                 >  2:992446/1‑239 (MQ=255)
                                                                                           ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                                                                                                <  2:1054183/239‑1 (MQ=255)
                                                                                                       agtatTTGATTGCGGCGTCGTTTGAGAATGCCGGGGTTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                    <  2:711828/237‑1 (MQ=255)
                                                                                                       aaTATTTGATTGCGGATTCTTTTGAGAATTCCGGGCGTTTTTAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCACCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                                                                                                <  1:876034/227‑1 (MQ=255)
                                                                                                       aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                    >  2:761317/1‑239 (MQ=255)
                                                                                                        tgggttGATTGCGGAGTAGTGACAGATTTCGGGGGCGTTTGAAGGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                    <  1:1123584/234‑1 (MQ=14)
                                                                                                        gtgttggATTGGGGATGGGTATGAGTATTGGGGGGTTTTAGGAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGTGCTTCGTCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATATGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                    <  1:1104171/232‑1 (MQ=14)
                                                                                                        atgtttGATTGGGGATTCGTTGCAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                                                                                                                      <  2:443350/233‑1 (MQ=255)
                                                                                                        atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                    <  2:43725/238‑1 (MQ=255)
                                                                                                          atTTGACTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                                                                                                                 <  2:1207552/239‑1 (MQ=255)
                                                                                                            tttaTTGCCGCTTGGTTAGAGGTTTCCGGTCGTTATGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGTTTCGCCCTTCGGGTCGTGGGCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                               <  1:129140/236‑1 (MQ=14)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCTTTGTCCAAGCTTCTCATCCTTCCCCGCATGCGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGTTGGTGGGGGAAGGATTAc                                                                                                                           >  1:1144992/1‑238 (MQ=255)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                                                                                                                      >  2:876034/1‑227 (MQ=255)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCACGCATGGGCAGAATATTTGATTGCGGATTGGCTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGCAGGATTAc                                                                                                                           >  1:726463/1‑238 (MQ=255)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCGGGATATTTTATTGCGGTTTCGTTTTACAGTTCAGGGACTTTTGAAGGTGATGGTGGTGGGGGAAGGAttgc                                                                                                                           >  1:629946/1‑236 (MQ=14)
                                                                                                                   cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGGGGGGGAAGGATTACTCa                                                                                                                        >  2:220113/1‑239 (MQ=255)
                                                                                                                                        gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGCTTCGCTTGCGAGTTCAGGGACTTTTGAAAGTGATGGAGGTGGGGGACGGATTTCTCAGCGCGTAGCGCTTCGCCCTTc                                                                                                   >  1:918250/1‑239 (MQ=14)
                                                                                                                                             tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                                                                                                                                >  2:1099238/1‑141 (MQ=37)
                                                                                                                                             tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                                                                                                                                <  1:1099238/141‑1 (MQ=37)
                                                                                                                                                        tggtgggggGGGGAGGAGTCCGCAGCGCTGGGCGCTTCGACCATCGGGGCGTTGGGTGCGGCATCGCTCTCGCGCTGGCGCTCAAGTCGACCCTTAGTCGTAGCTTATCATCCGTCCCCGCATGGGCAGAGTATCTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                    <  1:796376/238‑1 (MQ=11)
                                                                                                                                                               gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGTTCGTCGCCTGCGGCTACgg                                                                             >  2:218478/1‑237 (MQ=255)
                                                                                                                                                                       gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcga                                                                                                                          >  2:852719/1‑182 (MQ=255)
                                                                                                                                                                       gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcga                                                                                                                          <  1:852719/185‑4 (MQ=255)
                                                                                                                                                                           aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCATTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                 >  2:158751/1‑238 (MQ=255)
                                                                                                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                              <  1:95395/239‑1 (MQ=255)
                                                                                                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                              <  2:1027984/239‑1 (MQ=255)
                                                                                                                                                                              cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGTGTGTTTTATCAACGCTGGGCGCTTCCCCCGTCGGGCCTTGGGCTGGCGCAATCCTTCCTCGGTGGCGCt                                                              >  1:516663/1‑238 (MQ=14)
                                                                                                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTCGGGCATCGCTCTCTCGCTGgcccccccgtgga                                                      >  1:895741/1‑227 (MQ=255)
                                                                                                                                                                                         gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAcc                                                   <  2:597001/238‑1 (MQ=255)
                                                                                                                                                                                                         ccttGCCTGCGGCCACGCTGTCTCGCCGGCGGTGTAACGGTATCTTTGACGGACCGTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                   <  1:992446/236‑1 (MQ=14)
                                                                                                                                                                                                            tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCCCGAGTCGGACCTTTGTCGAAGCTTGTCAt                                 >  2:705375/1‑237 (MQ=255)
                                                                                                                                                                                                             gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGTCGAACCTTGGTCGTAGCTTCTCATctt                              >  2:447836/1‑237 (MQ=255)
                                                                                                                                                                                                                     cAACGCTCTCCCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc                       <  2:1194105/238‑1 (MQ=255)
                                                                                                                                                                                                                         gcgttctcGCCGACCTTCGAGACGAACCCTAGTCGAATCTTCTCAAGCGTCCGATGTGGAGAAGAGAAAGTGATTGGGGCTGCGCTTGAGACTTCGGGTACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                  <  1:403820/235‑1 (MQ=11)
                                                                                                                                                                                                                         gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                  <  2:594520/239‑1 (MQ=255)
                                                                                                                                                                                                                           tttctcGCTGACGGTCGAGGCGAACCTTAGTCGAAGCTTCTGATCCTTCCCCGGATGGGCAGAATATTTGATTGCGGATGCGTTTGAGAGTTCGGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                <  1:989003/237‑1 (MQ=255)
                                                                                                                                                                                                                           gtttcggcggcggcgcGAGCGGAATCTTAGAGGTAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                <  1:601896/224‑1 (MQ=14)
                                                                                                                                                                                                                                       gCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt    >  2:744713/1‑239 (MQ=255)
                                                                                                                                                                                                                                         tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTCTCGCGGGCGCTCGCAGCGAACCTTGGTCGAATCTGCTCATCCTTCCCCTCATGGCCAGAGTAGTTACtttc  >  1:1168618/1‑237 (MQ=14)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
GCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGC  >  NC_000913/1286736‑1287207

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