Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,479,936 A→C N267H (AAC→CAC)  ynbC → putative esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,479,9360AC98.7% 269.9 / ‑6.4 76N267H (AAC→CAC) ynbCputative esterase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (36/39);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGAACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGA  >  NC_000913/1479703‑1480167
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       
aTCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGAAc                                                                                                                                                                                                                                       >  1:206589/1‑236 (MQ=255)
 tCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGcaca                                                                                                                                                                                                                                      <  1:639565/236‑1 (MQ=255)
   cGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAg                                                                                                                                                                                                                                     >  2:45043/1‑235 (MQ=255)
     cTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAg                                                                                                                                                                                                                                     >  1:597791/1‑233 (MQ=255)
      tGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGgcg                                                                                                                                                                                                                                 <  2:73787/236‑1 (MQ=255)
             acGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCAt                                                                                                                                                                                                                          >  2:5560/1‑236 (MQ=255)
               gggCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCAttt                                                                                                                                                                                                                        >  1:48676/1‑236 (MQ=255)
                          gTTAATCTCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             <  2:388631/236‑1 (MQ=255)
                          gTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             <  1:464855/236‑1 (MQ=255)
                            tAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCaa                                                                                                                                                                                                            <  1:45043/235‑1 (MQ=255)
                                  cTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTtta                                                                                                                                                                                                     >  1:97527/1‑236 (MQ=255)
                                   ttGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTttat                                                                                                                                                                                                    >  1:109262/1‑236 (MQ=255)
                                     gCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             <  1:681831/225‑1 (MQ=255)
                                     gCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             >  2:681831/1‑225 (MQ=255)
                                      cTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAg                                                                                                                                                                                                 <  2:454383/236‑1 (MQ=255)
                                           tctcTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGtt                                                                                                                                                                                            <  2:48676/236‑1 (MQ=255)
                                               tACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTtata                                                                                                                                                                                        <  1:764570/236‑1 (MQ=255)
                                                aCAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGCCGTTTatat                                                                                                                                                                                       >  1:579847/1‑236 (MQ=255)
                                                    aaCGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGTTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTa                                                                                                                                                                                   >  1:436410/1‑236 (MQ=255)
                                                    aaCGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTa                                                                                                                                                                                   >  1:485292/1‑236 (MQ=255)
                                                         cTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGt                                                                                                                                                                                                <  2:256595/218‑1 (MQ=255)
                                                         cTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGt                                                                                                                                                                                                >  1:256595/1‑218 (MQ=255)
                                                                      taGAGTTGAGGGCGCGATTACGTTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTg                                                                                                                                                                 <  2:394918/236‑1 (MQ=255)
                                                                          gaTGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTa                                                                                                                                                             >  1:407892/1‑236 (MQ=255)
                                                                           aTGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACCAATCGCAACAATTTGATTa                                                                                                                                                             >  1:439817/1‑235 (MQ=255)
                                                                           aTGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTa                                                                                                                                                             >  1:265329/1‑235 (MQ=255)
                                                                             gcggcggcgATTACGCTCCCCACGCATCTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATc                                                                                                                                                           <  1:390637/235‑1 (MQ=255)
                                                                             gcggcggcgATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             >  2:591680/1‑185 (MQ=255)
                                                                             gcggcggcgATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             <  1:591680/185‑1 (MQ=255)
                                                                             gcggcggcgATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATc                                                                                                                                                           <  2:206589/235‑1 (MQ=255)
                                                                               ggcggcgATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGc                                                                                                                                                        >  2:766322/1‑236 (MQ=255)
                                                                                  ggcgATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCAt                                                                                                                                                      <  1:602807/235‑1 (MQ=255)
                                                                                       ttACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGa                                                                                                                                                 <  1:5560/235‑1 (MQ=255)
                                                                                       ttACGATCCCCACGAAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGa                                                                                                                                                 <  1:386552/235‑1 (MQ=255)
                                                                                          cGCTCCCCACGCATCTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTg                                                                                                                                              <  2:254422/235‑1 (MQ=255)
                                                                                             tCCCAACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGAATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGcaca                                                                                                                                          <  2:453669/236‑1 (MQ=255)
                                                                                               cccACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             >  1:229742/1‑167 (MQ=255)
                                                                                               cccACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCa                                                                                                                                                                                                             <  2:229742/167‑1 (MQ=255)
                                                                                                              atatCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCg                                                                                                                         <  1:766322/236‑1 (MQ=255)
                                                                                                                        aTGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATggcg                                                                                                                   >  1:726385/1‑232 (MQ=255)
                                                                                                                        aTGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTa                                                                                                               >  2:513182/1‑236 (MQ=255)
                                                                                                                               tgtgGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGca                                                                                                                                            <  2:724226/200‑1 (MQ=255)
                                                                                                                               tgtgGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGca                                                                                                                                            >  1:724226/1‑200 (MQ=255)
                                                                                                                                    tgCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGcaca                                                                                                                                          >  2:627584/1‑197 (MQ=255)
                                                                                                                                    tgCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGcaca                                                                                                                                          <  1:627584/197‑1 (MQ=255)
                                                                                                                                     gCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGAc                                                                                                  <  2:726385/236‑1 (MQ=255)
                                                                                                                                      cATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATggcgg                                                                                                                  >  1:469877/1‑219 (MQ=255)
                                                                                                                                      cATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATggcgg                                                                                                                  <  2:469877/219‑1 (MQ=255)
                                                                                                                                         cGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCg                                                                                               <  2:109262/235‑1 (MQ=255)
                                                                                                                                         cGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGt                                                                                              >  2:275766/1‑236 (MQ=255)
                                                                                                                                              aCCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCa                                                                                          >  2:265074/1‑235 (MQ=255)
                                                                                                                                                   agtttgattTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCg                                                                                    <  2:97527/233‑1 (MQ=255)
                                                                                                                                                   aGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCg                                                                                    <  1:234338/236‑1 (MQ=255)
                                                                                                                                                     attgattTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCCGCTGCCAGGCTTTTATCACGACCCGTTGGGTGAAGAGCACAGGGCGCTTGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAATAATTTGATTATCATCATGAAGACTGCACAGGACCATCAACGGATCGATGTCGGCTACGTTCTCGCGGACACCTGCCCTTTACGcc                                                                                  >  2:571874/1‑236 (MQ=255)
                                                                                                                                                               aTCCGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAGAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTa                                                                        <  2:513295/236‑1 (MQ=255)
                                                                                                                                                               aTCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTa                                                                        >  2:45256/1‑236 (MQ=255)
                                                                                                                                                                        tACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTa                                                                   <  1:265074/232‑1 (MQ=255)
                                                                                                                                                                          cGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTa                                                                        <  2:452276/225‑1 (MQ=255)
                                                                                                                                                                          cGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTa                                                                        >  1:452276/1‑225 (MQ=255)
                                                                                                                                                                               cccTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTACCAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCACCGAATCGATGGCGGTCACTTTCTGGTGGACCCGTCCCATTATCGCCGGTTGATTAACCGCATCCCTTTAttc                                                        >  2:73363/1‑234 (MQ=255)
                                                                                                                                                                                ccTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCg                                                       <  2:597791/236‑1 (MQ=255)
                                                                                                                                                                                ccTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCg                                                       >  2:56669/1‑236 (MQ=255)
                                                                                                                                                                                    tGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGATCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCgta                                                     >  2:118463/1‑232 (MQ=255)
                                                                                                                                                                                    tGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAg                                                   >  1:633605/1‑236 (MQ=255)
                                                                                                                                                                                            cTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGaa                                           >  1:125062/1‑236 (MQ=255)
                                                                                                                                                                                            cTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGaa                                           >  1:556023/1‑236 (MQ=255)
                                                                                                                                                                                              gCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAAt                                         >  2:265184/1‑236 (MQ=255)
                                                                                                                                                                                                 tCTGCTTCCAGGCTTTTATCACGTATCGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAGTCTTTATTAGTCGGTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTTCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGt                                      <  2:584738/236‑1 (MQ=255)
                                                                                                                                                                                                   tgctgcCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTc                                    <  1:513182/236‑1 (MQ=255)
                                                                                                                                                                                                       gcCAGGCTTTTATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCggg                                 >  1:413424/1‑235 (MQ=255)
                                                                                                                                                                                                                ttATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAACTGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCa                            <  2:413424/231‑1 (MQ=255)
                                                                                                                                                                                                                ttATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACg                              <  1:159903/229‑1 (MQ=255)
                                                                                                                                                                                                                ttATCACGACACGTTGGGTGAAGAGCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCa                            <  1:265184/231‑1 (MQ=255)
                                                                                                                                                                                                                ttAGCACGACACGTTGGTTGAAGAGCACAGGGCGCTGGCAATTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACg                              <  2:265329/229‑1 (MQ=255)
                                                                                                                                                                                                                                 gTGAAGTGCACAGGGCCTTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTAATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGa      <  2:439817/236‑1 (MQ=255)
                                                                                                                                                                                                                                     agagCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGcac                                                                                                                                           <  1:671834/99‑1 (MQ=255)
                                                                                                                                                                                                                                     agagCACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGcac                                                                                                                                           >  2:671834/1‑99 (MQ=255)
                                                                                                                                                                                                                                      gagcACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGa  >  1:13754/1‑235 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       
ATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGTCTGAACGTATTATTAGAGATGCGGCGGCGATTACGCTCCCCACGCAACTTCTGATATCAGGCGATGACTATGTGGTGCATCGCCAACCGCAGATTGATTTTTATCAGAGATTACGTAGCCCTCTGAAAGAGCTGCATCTGCTGCCAGGCTTTTATCACGACACGTTGGGTGAAGAGAACAGGGCGCTGGCATTTGAAAAAATGCAAAGCTTTATTAGTCGTTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGA  >  NC_000913/1479703‑1480167

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