Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,975,737 G→T 6.2% A198S (GCG→TCG)  wecE → dTDP‑4‑dehydro‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,975,7370GT 6.2% 323.7 / 7.7 96A198S (GCG→TCG) wecEdTDP‑4‑dehydro‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base T (4/2);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAAC  >  NC_000913/3975589‑3975884
                                                                                                                                                    |                                                                                                                                                   
tgtgGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGc                                                                                                                                                    <  1:1179412/150‑1 (MQ=255)
   ggTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGtt                                                                                                                                                 >  2:240377/1‑150 (MQ=255)
    gTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTa                                                                                                                                                >  2:465373/1‑150 (MQ=255)
      agaagaTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAAt                                                                                                                                              <  1:110042/150‑1 (MQ=255)
         agaTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGa                                                                                                                                           <  1:1188501/150‑1 (MQ=255)
           aTGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAAc                                                                                                                                         <  2:452354/150‑1 (MQ=255)
            tGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACg                                                                                                                                        >  2:1135536/1‑150 (MQ=255)
            tGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACg                                                                                                                                        <  2:146193/150‑1 (MQ=255)
                gCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGcc                                                                                                                                    >  2:1002133/1‑150 (MQ=255)
                        cGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatcat                                                                                                                            >  2:1518297/1‑150 (MQ=255)
                        cGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatcat                                                                                                                            >  2:382825/1‑150 (MQ=255)
                           gATGTCCACTTACAAAGGGCGGGCACTGGGAACCAGTGGGCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGGGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCg                                                                                                                         <  2:1071402/150‑1 (MQ=255)
                             tGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTg                                                                                                                       <  2:1087009/150‑1 (MQ=255)
                             tGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTg                                                                                                                       <  1:39474/150‑1 (MQ=255)
                              gTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGa                                                                                                                      >  1:486999/1‑150 (MQ=255)
                               tCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTAATc                                                                                                                                             <  2:776099/126‑1 (MQ=255)
                               tCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTAATc                                                                                                                                             >  1:776099/1‑126 (MQ=255)
                               tCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGaa                                                                                                                     >  1:901729/1‑150 (MQ=255)
                                ccACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGaaa                                                                                                                    >  1:1285737/1‑150 (MQ=255)
                                 cACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGaaaa                                                                                                                   <  1:1420082/150‑1 (MQ=255)
                                   cTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAgg                                                                                                                 <  2:1036342/150‑1 (MQ=255)
                                     tACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGc                                                                                                               >  2:407293/1‑150 (MQ=255)
                                     tACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGc                                                                                                               <  1:1135536/150‑1 (MQ=255)
                                       cAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCAc                                                                                                             >  2:614394/1‑150 (MQ=255)
                                       cAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCAc                                                                                                             <  1:917693/150‑1 (MQ=255)
                                       cAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCAc                                                                                                             >  2:626505/1‑150 (MQ=255)
                                         aaGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTa                                                                                                           >  2:994228/1‑150 (MQ=255)
                                            ggCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAAcc                                                                                                        <  1:482945/150‑1 (MQ=255)
                                             gCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatca                                                                                                                             >  2:1012003/1‑128 (MQ=255)
                                             gCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGatca                                                                                                                             <  1:1012003/128‑1 (MQ=255)
                                                tGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAg                                                                                                    >  2:1381379/1‑150 (MQ=255)
                                                 gCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGc                                                                                                   <  1:1493678/150‑1 (MQ=255)
                                                    cTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGAGCACTAACCGCAGCCAg                                                                                                >  2:643798/1‑150 (MQ=255)
                                                     tGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGt                                                                                               <  1:1168736/150‑1 (MQ=255)
                                                      gggAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTaa                                                                                                                                               <  2:895080/101‑1 (MQ=255)
                                                      gggAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGCTTTATCGTTaa                                                                                                                                               >  1:895080/1‑101 (MQ=255)
                                                       ggAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGttc                                                                                             >  1:1094563/1‑150 (MQ=255)
                                                            cATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCg                                                                                        <  2:1020885/150‑1 (MQ=255)
                                                              ttGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTg                                                                                      <  1:323413/150‑1 (MQ=255)
                                                               tGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTgg                                                                                     <  2:1130546/150‑1 (MQ=255)
                                                                ggTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGCGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGt                                                                                    >  1:1043396/1‑150 (MQ=255)
                                                                ggTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGt                                                                                    >  1:1216297/1‑150 (MQ=255)
                                                                  tCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCa                                                                                  <  2:1094563/150‑1 (MQ=255)
                                                                  tCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCa                                                                                  <  2:212501/150‑1 (MQ=255)
                                                                   cATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAg                                                                                 >  2:1174174/1‑150 (MQ=255)
                                                                    atatTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAgg                                                                                >  2:973948/1‑150 (MQ=255)
                                                                     tatTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGTGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGt                                                                               >  1:818297/1‑150 (MQ=255)
                                                                         ggCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGAt                                                                           >  2:199689/1‑150 (MQ=255)
                                                                         ggCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGAt                                                                           >  2:908416/1‑150 (MQ=255)
                                                                          gctgctTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATa                                                                          <  2:1508159/150‑1 (MQ=255)
                                                                          gctgctTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATa                                                                          <  1:382825/150‑1 (MQ=255)
                                                                           ctgctTTAGCTTCCATGAATCCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATaa                                                                         <  1:871873/150‑1 (MQ=255)
                                                                            tgctTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATaaa                                                                        <  1:1301135/150‑1 (MQ=255)
                                                                             gctTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAat                                                                       <  1:614394/150‑1 (MQ=255)
                                                                                 tAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATAc                                                                   >  2:742890/1‑150 (MQ=255)
                                                                                 tAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATAc                                                                   >  2:1037092/1‑150 (MQ=255)
                                                                                 tAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGATAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATAc                                                                   >  2:604654/1‑150 (MQ=255)
                                                                                       ccATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGgcg                                                             >  2:1471805/1‑150 (MQ=255)
                                                                                        cATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGgcgc                                                            >  2:542387/1‑150 (MQ=255)
                                                                                          tGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGa                                                          <  2:56941/150‑1 (MQ=255)
                                                                                           gAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGat                                                         >  1:534345/1‑150 (MQ=255)
                                                                                              aCCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATAtt                                                      >  1:1502858/1‑150 (MQ=255)
                                                                                               ccAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTg                                                     >  2:1397356/1‑150 (MQ=255)
                                                                                               ccAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTg                                                     <  2:369010/150‑1 (MQ=255)
                                                                                                cAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTgg                                                    <  1:742890/150‑1 (MQ=255)
                                                                                                    aaCTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTcc                                                <  1:931428/150‑1 (MQ=255)
                                                                                                         cacGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAat                                                                       >  2:1309836/1‑122 (MQ=255)
                                                                                                         cacGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAat                                                                       <  1:1309836/122‑1 (MQ=255)
                                                                                                             gcgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGTGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTg                                       <  1:643798/150‑1 (MQ=255)
                                                                                                              cgggcggTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGa                                      <  2:780831/150‑1 (MQ=255)
                                                                                                                    gTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCg                                >  2:1129690/1‑150 (MQ=255)
                                                                                                                    gTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCg                                >  1:1503701/1‑150 (MQ=255)
                                                                                                                        aGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCt                            <  1:973948/150‑1 (MQ=255)
                                                                                                                        aGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCt                            >  1:1081790/1‑150 (MQ=255)
                                                                                                                        aGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCt                            <  1:1172519/150‑1 (MQ=255)
                                                                                                                         ggcggcGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTg                           >  2:1503583/1‑150 (MQ=255)
                                                                                                                           cggcgCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCa                         <  2:459210/150‑1 (MQ=255)
                                                                                                                            ggcgCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCaa                        >  1:555281/1‑150 (MQ=255)
                                                                                                                             gcgcGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAg                       <  2:119004/150‑1 (MQ=255)
                                                                                                                              cgcgACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGc                      <  1:994228/150‑1 (MQ=255)
                                                                                                                                cgACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTg                    >  2:895586/1‑150 (MQ=255)
                                                                                                                                cgACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTg                    <  2:1530112/150‑1 (MQ=255)
                                                                                                                                  cgttcgATTAACGCTTTATCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCa                  >  1:934880/6‑150 (MQ=255)
                                                                                                                                   cGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCAt                 <  1:465373/150‑1 (MQ=255)
                                                                                                                                     cTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATAc               >  2:791236/1‑150 (MQ=255)
                                                                                                                                     cTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATAc               <  1:981562/150‑1 (MQ=255)
                                                                                                                                       gATTAACGCTTTATCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACAGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCt             >  1:9085/1‑150 (MQ=255)
                                                                                                                                       gATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCACGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCt             >  1:819261/1‑150 (MQ=255)
                                                                                                                                         ttAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCtgt           >  1:17616/1‑150 (MQ=255)
                                                                                                                                            aCGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTggg        >  1:28820/1‑150 (MQ=255)
                                                                                                                                            aCGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTggg        <  2:901729/150‑1 (MQ=255)
                                                                                                                                              gATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGgcg      <  1:518723/150‑1 (MQ=255)
                                                                                                                                               aTAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGgcgc     >  1:1218352/1‑150 (MQ=255)
                                                                                                                                                 aaaGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATaaa                                                                        >  2:214836/1‑81 (MQ=255)
                                                                                                                                                 aaaGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATaaa                                                                        <  1:214836/81‑1 (MQ=255)
                                                                                                                                                  aaGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAAc  >  2:407395/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTACACGGCGGGCGGTGAAGGCGGCGCGACGCTGATTAACGATAAAGCGTTAATCGAACGAGCCGAGATCATCCGTGAAAAGGGCACTAACCGCAGCCAGTTCTTCCGTGGTCAGGTCGATAAATATACCTGGCGCGATATTGGCTCCAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAAC  >  NC_000913/3975589‑3975884

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