Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,794,037 A→G K17K (AAA→AAG rfaD → ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,0370AG97.6% 125.0 / NA 41K17K (AAA→AAGrfaDADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/15);  minor base T (0/1);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATATTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGG  >  NC_000913/3793790‑3794281
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gggcaaattccagcacatattacccaaacttataggtcgcgacgagatgacgcgttaacacttctgcaaaattcaggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAACTCGTTAAGGcc                                                                                                                                                                                                                                                    >  2:806‑M2/150‑251 (MQ=255)
   caaattccagcacatattacccaaacttataggtcgcgacgagatgacgcgttaacacttctgcaaaattcaggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGGCTGAGATTGTCTCTGACTCCATAATGCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTg                                                                                                                                                                                                                                                 <  1:290334‑M2/105‑1 (MQ=255)
      atatcagcacatattacccaaacttataggtagcgacgagatgacgcgtcaacacttctgcaaaattcaggataaataacgtccatcaatcttccttttggcttatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGTCTCCATAATTCGAAGGTTTGAGTTATGATCATCGTTTCCGGCAGCGCGGGCTTTATCGGCAGCACCAACGTtaaggcccgtaat                                                                                                                                                                                                                                              <  2:9057‑M2/108‑6 (MQ=255)
         ccagcacagatgacccaaacgtagaagtcgcgacgagaggacgcgttaacacttcggcaaaattcaggataaataacgtccatcaatcgcacttttggcgcatggcggaagccatcgccgtgtagtgatatgatagaagtATTCGGGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTAGGATCATCGTTACCGGCGGCGCGGGCTTTAGCGGCAGCAACACCGTTAAGGCCCTGAATGa                                                                                                                                                                                                                                            <  2:305313‑M2/110‑1 (MQ=255)
               catattacccaaacttataggtcgcgacgagatgacgcgttaacacttctgcaaaattcaggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGc                                                                                                                                                                                                                                     <  2:138689‑M2/117‑1 (MQ=255)
                 tattacccaaacttataggtcgcgacgagatgacgcgttaacacttctgcaaaattcaggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCAt                                                                                                                                                                                                                                   >  1:166820‑M2/133‑251 (MQ=255)
                                     cgcgacgagatgacgcgataacacttctgcaaaattcaggataaataacgtccatcaatctcacgtttggcgaatggcggaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCtggtggtg                                                                                                                                                                                                                <  2:116261‑M2/138‑1 (MQ=255)
                                       cgacgagatgacgcgttaacacttctgcaaaattcaggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGAc                                                                                                                                                                                                             <  1:99722‑M2/141‑1 (MQ=255)
                                                        aacacttctgcaaaattcaggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCCACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCAc                                                                                                                                                                                            >  1:384328‑M2/94‑251 (MQ=255)
                                                              tctgcaaaattcaggataaataacgtccatcaatctcacatttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCATCATCGTTAATGCCCTGAATGATAAAGGCAGCACCGATATTCTGGTGTTGGACAACCTGAAAGACGGCACCAAGt                                                                                                                                                                                       <  2:45653‑M2/163‑1 (MQ=255)
                                                                   aaaattaaggataaataacgtccatcaatctcacttttggcgactggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACCTCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGGCAAACTGAAAAACGGCAACACGTTtgtg                                                                                                                                                                                  >  1:363603‑M2/83‑250 (MQ=255)
                                                                         caggataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCtg                                                                                                                                                                            <  1:189010‑M2/174‑1 (MQ=255)
                                                                             ataaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagaATTCGTGTCTGAGATTGTCTCTGACTCCATCATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCACTGCATGATAAAGGCATCCACGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCtggtgg                                                                                                                                                                        >  1:189585‑M2/73‑250 (MQ=255)
                                                                              taaataacgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGAt                                                                                                                                                                      >  1:73673‑M2/72‑251 (MQ=255)
                                                                                     cgtccatcaatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATAGGCAGCAAAATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGGGGTGGCCAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGAtctggct                                                                                                                                                                >  1:245934‑M2/65‑247 (MQ=255)
                                                                                             aatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGCAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCCTCACCGATATTCTGGTGGGGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGAc                                                                                                                                                       >  1:220194‑M2/57‑251 (MQ=255)
                                                                                                          gcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAgg                                                                                                                                          <  1:355643‑M2/208‑1 (MQ=255)
                                                                                                                 gctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTTTATGGATAAGGAAGACtt                                                                                                                                   >  1:436916‑M2/37‑251 (MQ=255)
                                                                                                                     aagccaacgccgtgaagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCctt                                                                                                                               >  1:138367‑M2/33‑250 (MQ=255)
                                                                                                                         caacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATcc                                                                                                                           >  1:238543‑M2/29‑251 (MQ=255)
                                                                                                                                            gatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAaga                                                                                                        <  1:806‑M2/242‑1 (MQ=255)
                                                                                                                                                 aagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtgg                                                                                                                  <  1:304495‑M2/232‑1 (MQ=255)
                                                                                                                                                 aagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCAACAAGGTTGTGAACCTGGGGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTtcctgatccagattatggggatcggaagggggt                                                                                                   >  2:304495‑M2/5‑236 (MQ=255)
                                                                                                                                                   gtATTCGCGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTCAGGCCCTGAATGCTAAAGGCATCACCGGTATTCTGGTGGGTGACAACCCGAAAGGCGGCACCAAGTTTGTGACCCTGGTGGGTCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCCGAttatggctggcggaaggtttggg                                                                                                 >  2:171138‑M2/3‑240 (MQ=255)
                                                                                                                                                               gagaTTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGGCAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGTTCGGCGCTGTCgcaacg                                                                                     >  2:289460/1‑246 (MQ=255)
                                                                                                                                                                  aTTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGTTCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGTTCGGCGAAGTCGAAGCggtt                                                                                  >  2:4637/1‑248 (MQ=255)
                                                                                                                                                                  aTTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAtt                                                                                  <  2:312104/251‑1 (MQ=255)
                                                                                                                                                                             cTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGGCTATATTGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAgg                                                                      >  1:34230/1‑251 (MQ=255)
                                                                                                                                                                                    aTTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAAGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGc                                                               <  2:73673/251‑1 (MQ=255)
                                                                                                                                                                                         aaGGTTACAGTTATGATCATCGTTACCGGCGGCGTGGGCTTTATCGGTCGCAACATCGTTAAGGCCCTGGATGATAAAGGCATCACCGATATTCTGGTGGGGGACAAACTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGAATATATGGATAAGGAAGACTTCATGATCCAGATTATGGCTGGCGAAAGGTTTGGAGATGTAGGAACGGTTTTCCACGAAgtgggatgctcttt                                                          >  2:20724/1‑238 (MQ=255)
                                                                                                                                                                                            gTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATCAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGACGCGATTTTCCACGAAGGCGCATGCTCTTcccc                                                       >  1:280302/1‑249 (MQ=255)
                                                                                                                                                                                                 aGTTAGGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACAGCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCg                                                  <  2:166820/251‑1 (MQ=255)
                                                                                                                                                                                                   ttATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGACCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATTGCGAAAGGATTTTTCCAGAAGGCGCATGCTCTTccaaccacgcg                                                >  2:353819/1‑243 (MQ=255)
                                                                                                                                                                                                              gTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGGGCATGCTCTTCCACCACCGAGTGGGATGGCa                                      >  1:205347/1‑250 (MQ=255)
                                                                                                                                                                                                                 aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTTAATGATAACGGCCTCCCCGATATTCTGGGGGTGGAACACCTGAAAGAAGGCACCAAGTTTGTGAAACTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGGGATGTCGGAGAGCTTTTTCACCAAAGCGCCTGCCCTTCCACCCCCgcgggggatggccagt                                   >  2:256821/1‑235 (MQ=255)
                                                                                                                                                                                                                 aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAATACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAGGGCAAGta                                  <  2:248359/251‑1 (MQ=255)
                                                                                                                                                                                                                 aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt                                   >  1:103130/1‑250 (MQ=255)
                                                                                                                                                                                                                            gggCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGAt                        >  1:14095/1‑250 (MQ=255)
                                                                                                                                                                                                                                tttATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATaaca                    >  1:313521/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ggCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCGCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCa             <  2:153269/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCATCGTTAAGGCACTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAATGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTc       <  1:345626/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGCGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCaaag  >  1:14012/1‑247 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AATATTGATACTAAAGCACTATTCACATGCAAAACCAACATCCGCCATGAAGGACTAGCTAAAACCCAAACTAGTTTGTTGCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGG  >  NC_000913/3793790‑3794281

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