Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,150,316 G→A E219K (GAG→AAG)  yghA → putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,150,3160GA95.0% 1225.2 / 35.7 380E219K (GAG→AAG) yghAputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base A (165/196);  total (175/205)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 385 total aligned reads displayed.

GGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTG  >  NC_000913/3150168‑3150463
                                                                                                                                                    |                                                                                                                                                   
  aaGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCgaagaa                                                                                                                                                 <  7:671726/151‑1 (MQ=255)
    gCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAg                                                                                                                                               <  2:808241/151‑1 (MQ=255)
      gATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGt                                                                                                                                             <  3:508301/151‑1 (MQ=255)
      gATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGt                                                                                                                                             <  3:1048192/151‑1 (MQ=255)
         cccccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTAtt                                                                                                                                          <  7:126102/151‑1 (MQ=255)
           cccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCg                                                                                                                                        <  4:1211183/151‑1 (MQ=255)
           cccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCg                                                                                                                                        >  1:1641373/1‑151 (MQ=255)
           cccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCg                                                                                                                                        >  7:1519213/1‑151 (MQ=255)
             gtGCTTCCGATAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAATTCCGCATTTACTGGTCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTTGCAAAACAGGTCGCGATTAAAGGTATTCggg                                                                                                                                      <  6:905282/150‑1 (MQ=255)
              tGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCACATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCgggg                                                                                                                                     >  6:19420/1‑150 (MQ=255)
               gCTACAGAAAGGTGAAAGTATCATCACCAATTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGt                                                                                                                                     <  2:1252432/150‑1 (MQ=255)
                  aCCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGaa                                                                                                                                  <  1:30322/150‑1 (MQ=255)
                    cGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAatat                                                                                                                               <  1:1118812/151‑1 (MQ=255)
                      aaaTGTGCAAGTATCATCACCACTTCGTCTATCCAGGGATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCGGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTg                                                                                                                             <  8:635711/151‑1 (MQ=255)
                      aaaGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTg                                                                                                                             <  1:1394201/151‑1 (MQ=255)
                       aaGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGt                                                                                                                            >  8:1642966/1‑151 (MQ=255)
                       aaGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGAGAAGAAAGGTATTCGGGTGAATATTGt                                                                                                                            >  1:1052089/1‑151 (MQ=255)
                         ggTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCATGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTcg                                                                                                                          <  7:925839/151‑1 (MQ=255)
                             cAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGcc                                                                                                                      >  4:432618/1‑151 (MQ=255)
                              aaGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCa                                                                                                                     <  5:1375996/151‑1 (MQ=255)
                              aaGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCa                                                                                                                     <  2:219358/151‑1 (MQ=255)
                              aaGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCa                                                                                                                     <  5:1599604/151‑1 (MQ=255)
                                     atcaCCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGa                                                                                                              <  4:1355970/151‑1 (MQ=255)
                                       caccacTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATc                                                                                                            <  6:800468/151‑1 (MQ=255)
                                        accacTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCt                                                                                                           <  5:475924/151‑1 (MQ=255)
                                             ttCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGAcc                                                                                                      >  4:145782/1‑150 (MQ=255)
                                              tCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACa                                                                                                      <  5:1546001/150‑1 (MQ=255)
                                                  cAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACt                                                                                                 >  2:989216/1‑151 (MQ=255)
                                                   aaTCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTg                                                                                                >  1:1030798/1‑151 (MQ=255)
                                                    aTCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGATGGCGGGGATTCTGAACTACAGCCGTGGCTTGGCTAAGCAGGTGGGGATGAAGGGGGTTGGGGTGATGATTGTCGGGGCAGGGCCGCTGTgggggggggggg                                                                                               >  2:1388744/1‑141 (MQ=255)
                                                        aGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAAGAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGCCAGCACTGCAAAt                                                                                           >  5:230742/1‑151 (MQ=255)
                                                          gCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAAttt                                                                                         >  4:269367/1‑151 (MQ=255)
                                                          gCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAAttt                                                                                         >  7:390479/1‑151 (MQ=255)
                                                          gCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCACGCCCGATCTGGACAGCACAGCAAAttt                                                                                         >  2:754806/1‑151 (MQ=255)
                                                           cATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTc                                                                                        >  8:1077801/1‑151 (MQ=255)
                                                           cATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTc                                                                                        >  1:1590296/1‑151 (MQ=255)
                                                              aCCAGCCAAGTCCTCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAATCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcg                                                                                      <  5:1235161/150‑1 (MQ=255)
                                                              aCCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATGCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcg                                                                                      <  1:1316165/150‑1 (MQ=255)
                                                               ccAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcggc                                                                                    <  5:977572/151‑1 (MQ=255)
                                                                    cAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATACTCAACAACAGACGTGGCTTGGCAAAACAGGTCGCGAAGAAACGAATTCGGGTGAATATTGTCGCGCCACGCCCGATCTGGACAGCACTCCAAATTTCCGGCCGAc                                                                                >  4:992061/1‑150 (MQ=255)
                                                                     aaGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACa                                                                               >  7:1170250/1‑150 (MQ=255)
                                                                       gTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAAc                                                                            <  6:613063/151‑1 (MQ=255)
                                                                       gTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAAc                                                                            <  3:1274297/151‑1 (MQ=255)
                                                                         ccGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACg                                                                           <  5:130452/150‑1 (MQ=255)
                                                                         ccGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGc                                                                          >  5:72976/1‑151 (MQ=255)
                                                                            cATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAgg                                                                       >  3:1651595/1‑151 (MQ=255)
                                                                             aTTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAgg                                                                       >  1:211738/1‑150 (MQ=255)
                                                                             aTTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAgg                                                                       >  8:1504726/1‑150 (MQ=255)
                                                                             aTTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATAGTCGCGCCAGGCCCTATCTGGACAGCACTGGAAAATTCAGTCGTACAAACGCacg                                                                       >  8:1389671/1‑148 (MQ=255)
                                                                               ttACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGGCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGAt                                                                     <  6:985123/150‑1 (MQ=255)
                                                                               ttACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGAt                                                                     >  2:357569/1‑150 (MQ=255)
                                                                                tACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATa                                                                    <  5:813659/150‑1 (MQ=255)
                                                                                  cTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATa                                                                    >  6:1212395/1‑148 (MQ=255)
                                                                                      aCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATccc                                                             <  5:56669/151‑1 (MQ=255)
                                                                                        tATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGc                                                           >  5:422217/1‑151 (MQ=255)
                                                                                         aTGCGGCTACGAAGGCGGCGATTCAGAACTACAGCCGTGGCTTGGCAAAACAGCTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTACCGGCGGACAAACGCAGGATAAGATCCCg                                                            <  2:594334/149‑1 (MQ=255)
                                                                                          tGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGc                                                           >  2:495088/1‑149 (MQ=255)
                                                                                            cGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGtt                                                       <  7:459320/151‑1 (MQ=255)
                                                                                             ggCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGttt                                                      >  8:732161/1‑151 (MQ=255)
                                                                                             ggCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGttt                                                      >  4:1603737/1‑151 (MQ=255)
                                                                                              gCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTg                                                     <  7:1426245/151‑1 (MQ=255)
                                                                                                 aCGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGt                                                   <  1:931652/150‑1 (MQ=255)
                                                                                                 aCGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                  >  8:495293/1‑151 (MQ=255)
                                                                                                 aCGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTgg                                                    <  3:751642/149‑1 (MQ=255)
                                                                                                 aCGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                  <  6:696017/151‑1 (MQ=255)
                                                                                                 aCGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                  <  2:1137483/151‑1 (MQ=255)
                                                                                                 aCGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                  >  2:1520417/1‑151 (MQ=255)
                                                                                                  cGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGAGAGCTCTGCATATTTCGGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCa                                                 <  7:270115/151‑1 (MQ=255)
                                                                                                   gAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                  >  4:757839/1‑149 (MQ=255)
                                                                                                             aTTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGa                                      <  6:1535633/151‑1 (MQ=255)
                                                                                                              ttCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGAt                                     >  3:2013/1‑151 (MQ=255)
                                                                                                               tCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCTGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGAt                                     <  3:86622/150‑1 (MQ=255)
                                                                                                                 tGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGa                                   >  5:887509/1‑150 (MQ=255)
                                                                                                                       aCAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGc                            >  4:1020977/1‑151 (MQ=255)
                                                                                                                        cAGCCGTGGCCTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCg                           <  6:1068500/151‑1 (MQ=255)
                                                                                                                         cgCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCgg                          <  1:758266/150‑1 (MQ=255)
                                                                                                                         aGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCgg                          <  3:43800/151‑1 (MQ=255)
                                                                                                                           ccGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCgggg                        <  2:270995/151‑1 (MQ=255)
                                                                                                                             gTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGc                       >  8:37123/1‑150 (MQ=255)
                                                                                                                               ggCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTTAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCaa                     <  1:647288/150‑1 (MQ=255)
                                                                                                                               ggCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCaa                     <  8:864897/150‑1 (MQ=255)
                                                                                                                                   tGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACcg                  <  1:592037/149‑1 (MQ=255)
                                                                                                                                   tGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACcg                  <  3:1173762/149‑1 (MQ=255)
                                                                                                                                    ggCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACcggcg               >  7:7204/1‑151 (MQ=255)
                                                                                                                                    ggCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACcggcg               <  8:1254042/151‑1 (MQ=255)
                                                                                                                                      cAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCATCAAACGCCGATGAAACGTGCGGGGCAACcggcgg              >  1:412749/1‑150 (MQ=255)
                                                                                                                                      cAAAACAGCTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGa             <  4:671183/151‑1 (MQ=255)
                                                                                                                                       aaaaCAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGaa            <  7:213871/151‑1 (MQ=255)
                                                                                                                                        aaaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAAc           >  6:39015/1‑151 (MQ=255)
                                                                                                                                        aaaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAAc           >  4:920213/1‑151 (MQ=255)
                                                                                                                                         aaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACt          >  3:1493901/1‑151 (MQ=255)
                                                                                                                                         aaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACt          <  3:992212/151‑1 (MQ=255)
                                                                                                                                         aaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACt          <  6:141268/151‑1 (MQ=255)
                                                                                                                                           cAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTgg        >  7:79693/1‑151 (MQ=255)
                                                                                                                                            aGGTCGCGGAGAAAGGTATTCGGGGGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGc       <  5:1118880/151‑1 (MQ=255)
                                                                                                                                            aGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGc       <  1:499560/151‑1 (MQ=255)
                                                                                                                                            aGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGc       <  1:116747/151‑1 (MQ=255)
                                                                                                                                                cgcgAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCt   <  3:605939/151‑1 (MQ=255)
                                                                                                                                                cgcgAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCt   <  6:1333006/151‑1 (MQ=255)
                                                                                                                                                cgcgAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGAAACCGGCGGAACTGGCCCCt   <  2:421806/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTG  >  NC_000913/3150168‑3150463

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